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A Symposium of Large Scale Digital 
Calculating Machinery 


Under the joint sponsorship of Harvard University and the Bureau of 
Ordnance, United States Navy, 336 representatives of university, industrial 
and government laboratories and research groups met in Cambridge, Mass. 
from January 7—10, 1947 to participate in a four-day Symposium on Large 
Scale Digital Calculating Machinery. 

The meetings were held in Harvard’s new Computation Laboratory and, 
in part, the Symposium celebrated the formal opening of this new facility. 

In addition to the extensive group of technical papers, the program 
included a demonstration of the Automatic Sequence Controlled Calculator 
(usually referred to as the Mark I Calculator), and a preview of the Mark II 
Calculator which is now being assembled at Harvard for the Bureau of 
Ordnance. The latter machine will later be reassembled at the Dahlgren 
Proving Ground of the Navy. 

The program of the Symposium was as follows: 


I. Tuesday, 7 January, 10:00 a.m. 
Opening Addresses by 
Edward Reynolds, Administrative Vice-President of Harvard University 
Rear Admiral C. T. Joy, USN, Naval Proving Ground, Dahlgren, Virginia 
Prof. H. H. Aiken, Technical Director of the Computation Laboratory, Harvard 
University 
Inspection of the Computation Laboratory and of Mark II Calculator 


II. Tuesday, 7 January, 2:00 p.m. 
Existing Calculating Machines, Prof. W. E. Bleick, Chairman 
. “The Work of Charles Babbage”’ by Mr. R. H. Babbage? 
. “Mark I Calculator” by Mr. R. M. Bloch 
. “Brief description and operating characteristics of the ENIAC” by Dr. L. P. Tabor 
. “Bell Telephone Laboratories relay computing systems” by Mr. S. B. Williams 
. “Mark II Calculator’ by Mr. R. V. D. Campbell 


III. Wednesday, 8 January, 9:30 a.m. 
The Logic of Large Scale Calculating Machines, Prof. W. H. Furry, Chairman 

1. ‘Problems of mathematical analysis involved in machine computations” by Dr. A. W. 
Wundheiler 

2. “The organization of large scale calculating machinery” by Dr. G. R. Stibitz 


IV. Wednesday, 8 January, 2:00 p.m. 
Storage Devices, Dr. J. H. Curtiss, Chairman 
. “Mercury delay lines as a memory unit’’ by Dr. T. K. Sharpless 
. “Slow electromagnetic waves” by Prof. L. N. Brillouin 
. “High Speed electrostatic storage” by Mr. J. W. Forrester 
. “Magnetic and phosphor coated disks” by Dr. B. L. Moore 
. “The selectron—a tube for selective electrostatic storage” by Dr. Jan Rajchman 
. “Optical and photographic storage techniques” by Dr. A. W. Tyler 


V. Thursday, 9 January, 9:30 a.m. 


Numerical Methods and Suggested Problems for Solution, Dr. Mina S. Rees, Chairman 
. ‘Method of finite differences for the solution of partial differential equations” by Prof. 
Richard Courant 
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“On computational techniques for certain problems in fluid dynamics” by Dr. R. J. Seeger 


. “Computational problems arising in connection with economic analysis of interindustrial 


relationships” by Prof. W. W. Leontief 


. “On the accumulation of errors in processes of integration on high speed calculating 


machines” by Prof. H. A. Rademacher 


. “Fluid mechanics computations” by Prof. H. W. Emmons 
. “Firing tables” by Dr. L. S. Déderick 


VI. Thursday, 9 January 2:00 p.m. 
Sequencing, Coding, and Problem Preparation, Prof. J. A. Stratton, Chairman 


. “Coding for large scale calculating machinery’ by Dr. H. H. Goldstine 
. “Preparation of problems for EDVAC-type machines’’ by Dr. J. W. Mauchly 
. “The preparation of problems for large scale calculating machinery” by Mr. J. O. 


Harrison, Jr. 


VII. Friday, 10 January, 9:30 a.m. 
Input and Output Devices, Prof. E. L. Chaffee, Chairman 


. “Application of printing telegraph techniques to large scale calculating machinery” by 


Mr. F. G. Miller 


. “Some physical aspects of magnetic recording’ by Mr. Otto Kornei 
. “The numeroscope” by Mr. H. W. Fuller 
. “Input and Output Devices for electronic digital calculating machinery” by Mr. S. N. 


Alexander 


. “An input device using multiple gates’’ by Dr. Morris Rubinoff 
. “Photographic methods of handling input and output data” by Dr. K. G. Macleish 
. “Transfer between external and internal memory” by Mr. C. B. Sheppard 


VIII. Friday, 10 January, 2:00 p.m. 
Conclusions and Open Discussion, Prof. C. C. Bramble, Chairman 


. “Publication, classification, and patents” by Prof. S. H. Caldwell 
. “Le domaine du calcul mécanique” by Dr. Louis Couffignal (paper not read, but to be 


included in the published report). 


. “New Vistas in mathematics” by Dr. A. T. Waterman 


All of the remarks and papers of the 34 speakers, as well as the related 


discussions, are to be published in a later volume of the Annals of the 
Computation Laboratory. 


Both Harvard University and the Bureau of Ordnance are to be con- 


gratulated for their enterprise in making possible this significant meeting. 
The extensive program reflected every phase of the sweeping progress being 
made in the field of large-scale calculators, and the large attendance offered 
a unique opportunity for the numerous smail group meetings which do much 
to facilitate the exchange of helpful information. 


The following members registered for the Symposium: 


Milton Abramowitz, mathem., NBSMTP 

Erik Ackerlind, group leader, Northrop Aircraft Corp., Hawthorne Field, Cal. 

L. V. Ahlfors, prof. math., Harvard Univ. 

H. H. Aiken, prof. appl. math., Harvard Univ. 

L. W. Alberts, grad. stu., Harvard Law Sch. 

S. N. Alexander, group leader, computer components group, National Bureau of Standards 


(NBS) 
Alterman, engin., Reeves Instrument Corp., New York 


F. J. 
E. G. Andrews, engin., Bell Tel. Lab., New York 
R. C 


. Archibald, prof. math., Brown Univ. 
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George Arfken, Jr., grad. stu., Yale Univ. 

K. J. Arnold, assist. prof. math., Univ. Wisconsin 

R. H. Babbage, editor, 241 St. James St. West, Montreal, Canada 

Rear Admiral W. D. Baker, commandant, U. S. Naval Base, Boston 

M. D. Baller, group leader, computers, Army Air Force (AAF), Malden, Mass. 
Colin Barnes, prof. physics, Univ. Toronto 

Dr. J. L. Barnes, consultant, Raytheon Mfg. Co., Waltham, Mass. 

R. M. Barrett, electronic engin., Watson Lab., AAF AMC, Cambridge 

Phyllis P. Barrows, res. assist., Harvard Univ. 


C. E. Barthel, Jr., electronic computing and control group, Naval Ordnance Lab., Wash. 


Walter Bartky, prof. math., Univ. Chicago. 

Dr. A. E. Benfield, visiting lect. appl. physics, Harvard Univ. 

William Bentinck-Smith, ed. Harvard Alumni Bull., Harvard Univ. 

R. J. Bergemann, electronic engin., Office Naval Res., Boston 

Dr. Stefan Bergman, res. lect. aeronautical engin., Harvard Univ. 

E. C. Berkeley, methods analyst, Prudential Life Ins. Co., Newark, N. J. 

A. A. Berlinsky, lab. electrician, Bureau of Census, Wash. 

C. H. Berry, prof. mech. engin., Harvard Univ. 

J. H. Bigelow, chief engin., Institute Adv. Studies, Princeton, N. J. 

Garrett Birkhoff, prof. math., Harvard Univ. 

W. E. Bleick, assoc. prof. math., post grad. sch., U. S. Naval Acad., Annapolis, Md. 
Daniel Blitz, res. engin., Raytheon Mfg. Co., Waltham, Mass. 

R. M. Bloch, res. assist., Harvard Univ. 

J. H. Blythe, nautical scientist, Hydrographic Office, Navy Dept., Wash. 
Richard Blythe, res. assoc. electr. engin., Mass Inst. Techn. 

J. M. Boermeester, actuarial assoc., John Hancock Mutual Life Ins. Co., Boston 
B. A. Boley, instr. aeron. engin., Polytechnic Inst., Brooklyn, N. Y. 

Capt. C. T. Bonney, prof. naval sci. and tactics, Harvard Univ. 

H. R. Boyd, engin., servomechanisms lab., Mass. Inst. Techn. 

C. C. Bramble, prof. math., post-grad. sch., U. S. Naval Acad. 

Dr. A. E. Brandt, statist., Naval Ordnance Lab., Wash. 

Mark Breiter, mathem., office chief ordnance, Pentagon Bldg., Wash. 

Ruth A. Brendel, instr. math., Univ. Buffalo 

P. W. Bridgman, prof. math. and nat. phil., Harvard Univ. 

L. N. Brillouin, prof. appl. math., Harvard Univ. 

D. R. Brown, electronic engin., servomechanisms lab., Mass. Inst. Techn. 

O. E. Brown, assoc. prof. math., Case Sch. Appl. Sci., Cleveland, O. 

R. L. Brown, research assoc., Mass. Inst. Techn. 

T. H. Brown, prof. business statistics, Harvard Univ. 

Capt. C. A. Buchanan, res. and develop. div., Navy Dept., Bur. Ord., Wash. 


Dr. R. S. Burington, chief mathem., evaluation and analysis group, Navy Dept., Bur. 


Ord., Wash. 
S. H. Caldwell, prof. electr. engin., Mass. Inst. Techn. 
H. Campaigne, mathem., Office Chief Naval Operations, Wash. 
Commander A. A. A. Campbell, Recruiting Station, First Naval District, Boston 
R. V. D. Campbell, res. assoc., Harvard Univ. 
W. L. Campbell, prof. food techn., Mass. Inst. Techn. 
Dr. E. W. Cannon, chief, machinery devel. sect., NBS 
Achille Capecelatro, physicist, Underwater Sound Lab., New London, Conn. 
F. E. Carey, fellow in journalism, Harvard Univ. 


Commander J. H. Carmichael, res. and development div., Navy Dept., Bur. Ord., Wash. 


Arthur Casagrande, assoc. prof. civil engin., Harvard Univ. 

E. L. Chaffee, prof. physics and commun. engin., Harvard Univ. 
Harold Chestnut, engin., Gen. Electr. Co., Schenectady, N. Y. 
B, F, Cheydleur, analysis sect., Naval Ord. Lab., Wash, 
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Chuan Chu, res. assoc., Univ. Penn. 

W. H. Claflin, Jr., treasurer Harvard Coll., 24 Milk St., Boston 

G. M. Clemence, director, Naut. Alm. Off., Naval Obs., Wash. 

Dr. R. F. Clippinger, mathem., Ballistic Res. Lab., Aberdeen Proving Gr., Md 
Admiral W. T. Cluverius, president, Worcester Polyt. Inst. 

Charles Concordia, engin., Gen. Electr. Co., Schenectady, N. Y. 
Patrick Conley, teaching fellow appl. physics, Harvard Univ. 

J. M. Coombs, engin., Engineering Res. Assoc. Inc., St. Paul, Minn. 
Lester Corrsin, grad. stu., Mass. Inst. Techn. 

Mary Coulan, admin. assist., Harvard Law School 

Richard Courant, prof. math., New York Univ. 


Perry Crawford, Jr., director, Spec. Devices Center, Office Naval Res., Port Washington, 


L. I. 
Dr. J. H. Curtiss, assist. to the director, NBS 
H. N. Cutler, stu., Harvard Coll. 
J. F. Daly, statist., Bureau of Census, Wash. 
Dr. G. B. Dantzig, techn. adv. (math.), AAF Office Air Comptr., Pentagon Bldg., Wash. 
D. J. Davis, stu., Harvard Coll. 
Dr. L. S. Dederick, assoc. director, Ballistic Res. Lab., Aberdeen Proving Gr., Md. 
R. L. DeVeer, sales engin., C. P. Clare Co., Boston 
Rear Admiral M. L. Deyo, commandant, First Naval Distr., Boston 
A. H. Dickinson, engin., Watson Lab., IBM Corp., New York 
S. H. Dodd, Jr., res. engin., Mass. Inst. Techn. 
Dr. F. W. Dresch, principal mathem., Naval Proving Gr., Va. 
B. M. Durfee, engin., IBM Corp., Endicott, N. Y. 
Dr. J. J. Eachus, Office Chief Naval Oper., Navy Dept., Wash. 
Dr. W. J. Eckert, director pure sci., Watson Lab., IBM Corp., New York 
N. E. Edlefsen, engin., Raytheon Mfg. Co., Waltham, Mass. 
Murray Ellis, res. assist., Harvard Univ. 
G. V. Eltgroth, divisional patent attorney, Bendix Aviation Corp., Riatnnnin Md. 
H. W. Emmons, assoc. prof. mech. engin., Harvard Univ. 
Bernard Epstein, res. assoc. aeronaut. engin., Harvard Univ. 
R. R. Everett, res. engin., Mass. Inst. Techn. 
Harris Fahnestock, res. engin., Mass. Inst. Techn. 
Dr. J. H. Faull, Jr., head sci. sect., Office Naval Res., Boston 
Samuel Favor, mech. engin., Computation Lab., Harvard Univ. 
H. Feistel, physicist, Watson Lab., AAF AMC, Cambridge, Mass. 
L. A. Fine, instr. math., Harvard Univ. 
H. W. Flagg, sales repr., Marchant Calc. Machine Co., Boston 
D. A. Flanders, assoc. prof. math., New York Univ. 
Margaret I. Florencourt, res. engin., servomechanisms lab., Mass. Inst. Techn. 
N. C. Ford, assist. actuary, Mass. Mutual Life Ins. Co., Springfield, Mass. 
J. W. Forrester, assoc. dir., servomechanisms lab., Mass. Inst. Techn. 
Donald Foster, physicist, Watson Lab., AAF AMC, Cambridge 
L. V. Foster, optical engin., Bausch & Lomb Co., Rochester, N. Y. 
P. G. Frank, lect. physics and math., Harvard Univ. 
Stanley Frankel, assist. prof., Univ. Chicago 
L. F. Freedman, mech. engin., Reeves Instr. Corp., New York 
A. E. Frost, assist. to automatics res. engin., Western Union Tel. Co., New York 
H. W. Fuller, res. assist., Harvard Univ. 
W. H. Furry, assoc. prof. physics, Harvard Univ. 
I. R. Gale, expediter purchaser, Raytheon Mfg. Co., Waltham, Mass. 
Herbert Galman, physicist, Exper. Br., Frankford Arsenal, Philadelphia, Pa. 


Dr. H. L. Garabedian, prin. physicist, Clinton Labs., Monsanto Chem. Co, Knoxville, Tenn. 


C. S. Gardner, stu., Harvard Univ. 
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E. B. Gardner, operations evaluations group, Mass. Inst. Techn. 
Dr. H. F. Gingerich, mathem., Naval Communi. Annex, Wash. 


’ Capt. J. A. Glick, insp. naval material, 495 Summer St., Boston 


R. H. Goddard, develop. engin., Sylvania Electric Products, Inc., Boston 

J. A. Goff, dean, Towne Sci. Sch., Univ. Penn. 

Dr. H. E. Goheen, mathem., Office Naval Res:, Boston 

Dr. H. H. Goldstine, res. fellow, Institute Adv. Studies, Princeton, N. J. 
Harold Grad, grad. stu., New York Univ. 

R. B. Graham, chief engin. SPDE, Bendix Aviation Corp., Teterboro, N. J. 
J. V. Granger, grad. stu., Harvard Univ. 

E. F. Grant, radio engin., Watson Lab., AAF AMC, Cambridge 

Serge Gratch, instr. mech. engin., Univ. Penn. 

H. V. Greenough, Jr., director, Technichord Records, Cambridge 

R. H. Griest, res. engin., Hughes Aircraft Corp., Culver City, Cal. 

W. C. Griffith, teaching fellow engin. sci., Harvard Univ. 

S. E. Gudmundson, civil engin., Aeronaut. Res. Inst. of Sweden, Ulvsunde, Sweden 
G. Hagg, prof., Univ. Upsala, Sweden 

B. Hamilton, superv. IBM, Mass. Gen. Hosp., Boston 

F. E. Hamilton, engin., IBM Corp., Endicott, N. Y. 

Dr. R. W. Hamming, engin., Bell Tel. Labs., Murray Hill, N. J. 

E. B. Hammond, Jr., res. -engin., Sperry Gyroscope Co., Great Neck, L. I. 

M. H. Hansen, statist. assist. to director, Bureau of Census, Wash. 

J. O. Harrison, Jr., res. assoc., Harvard Univ. 

R. W. Hart, electr. engin., Office Naval Res., Boston 

Bernhard Haurwitz, assoc. prof. meteorology, Mass. Inst. Techn. 

Capt. E. C. Hawk, off. superv. insp. naval material, First Naval District, Boston 
R. L. Hawkins, chief technician, Computation Lab., Harvard Univ. 

H. L. Hazen, prof. electr. engin., Mass. Inst. Techn. 

Dr. R. W. Hickman, lect. appl. physics, Harvard Univ. 

G. W. Hobbs, electronics engin., Gen. Electric Co., Schenectady, N. Y. 

J. V. Holberton, Ballistics Res. Lab., Aberdeen Proving Gr., Md. 

L. M. Hollingsworth, chief components lab., Watson Lab., AAF AMC, Cambridge 
R. E. Hopkins, assoc. prof., Univ. Rochester 

Dr. Grace M. Hopper, res. fellow, Harvard Univ. 

A. B. Horton, Jr., res. assist. electr. engin., Mass. Inst. Techn. 

Dr. A. S. Householder, prin. physicist, Clinton Laboratories, Oak Ridge, Tenn. 
J. H. Howard, director lab., Engineering Res. Associates Inc., St. Paul, Minn. 
B. C. Howland, reporter, Boston Post 

E. S. Hughes, Jr., designer, IBM Corp., Endicott, N. Y. 

F. B. Hunt, jun. engin., N. E. Power Service Co., Boston 

F. V. Hunt, assoc. prof. appl. physics, Harvard Univ. 

Dr. J. M. Ide, techn. director, Underwater Sound Lab., New London, Conn. 
Fritz John, assoc. prof. appl. math., New York Univ. 

Rear Admiral C. T. Joy, commanding officer, Naval Proving Gr., Dahlgren, Va. 
H. R. Kaiser, res. engin., Univ. Calif., Los Angeles 

G. N. Kamm, grad. stu., Harvard Univ. 

Martin Katzin, consultant, Naval Res. Lab., Wash. 

R. Kelner, stu., Harvard Univ. 

Constance Kepke, res. assist., Harvard Univ. 

Capt. C. H. J. Keppler, counsellor foreign students, Harvard Univ. 

Dr. R. B. Kershner, appl. physics lab., The Johns Hopkins Univ. 

Dr. G. W. King, res. chemist, Arthur D. Little Inc., Cambridge 

R. W. P. King, assoc. prof. appl. physics, Harvard Univ. 

Bertram Klein, res. assoc., Polytechnic Inst., Brooklyn, N. Y. 

Ruth A. Knowlton, secretary, Harvard Univ. 
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Dr. Zdenek Kopal, res. assoc. electr. engin., Mass. Inst. Techn. 

Otto Kornei, res. engin., The Brush Development Co., Cleveland, O. 

Capt. W. E. Kraft, superv. inspec. naval material, First Naval District, Boston 

George Kramer, reporter, The Harvard Crimson 

J. A. Krumhansl, assist. prof. physics, Brown Univ. 

Jack Laderman, mathem., NBSMTP 

Emory Lakatos, engin., Bell Tel. Labs., Murray Hill, N. J. 

C. D. Lake, development engin., IBM Corp., Endicott, N. Y. 

Dr. O. E. Lancaster, mathem., Bureau Aeronautics, Navy Dept., Wash. 

Dr. Philippe LeCorbeiller, lect. appl. physics, Harvard Univ. 

Dr. Joseph Lehner, mathem., Hydrocarbon Research Inc., New York 

Alfred Leitner, res. assist., Yale Univ. 

W. W. Leontief, prof. economics, Harvard Univ. 

Dr. A. A. Levin, physicist, Bausch & Lomb Optical Co., Rochester, N. Y. 

Dr. J. H. Levin, chief machines sect. computing lab., Ballistic Res. Lab., Aberdeen Proving 
Gd., Md. 

M. L. Libman, patent attorney, Office Chief Ord., War Dept., Wash. 

Dr. D. M. Little, secy. to univ., Harvard Univ. 

Barbara Livingston, secy. to dean, Harvard Law School 

W. S. Loud, res. engin., Mass. Inst. Techn. 

Samuel Lubkin, engin., Ballistic Res. Lab., Aberdeen Proving Gd., Md. 

Herman Lukoff, res. assist., Univ. Penn. 

J. R. Macdonald, res. engin., Mass. Institute Techn. 

Capt. Louis McGlone, director training, First Naval District, Boston 

G. W. Mackey, assist. prof. math., Harvard Univ. 

Saunders MacLane, prof. math., Harvard Univ. 

Dr. K. G, Macleish, physicist, Eastman Kodak Co., Rochester, N. Y. 

A. B. Macnee, res. assist. electr. engin., Mass. Inst. Techn. 

Kathleen McNulty, Ballistic Res. Lab., Moore Sch., Univ. Penn. 

J. L. McPherson, statistician, Bureau of Census, Wash. 

C. Malik, electronic engin., Watson Lab., AFC AMC, Cambridge 

R. A. Mangini, manager planning dept., John Hancock Mutual Life Ins. Co., Boston 

Charles Manneback, prof. theoretical physics, Univ. Louvain 

L. C. Mansur, radar engin., Watson Lab., AAF AMC, Cambridge 

T. M. Marcus, NRC fellow, Mass. Inst. Techn. 

L. J. Marsee, planning engin., Mass. Mutual Life Ins. Co., Springfield, Mass. 

R. C. Martin, Jr., stu., Harvard Univ. 

A. H. Mathewson, purchasing agent, Mass. Gen. Hospital, Boston 

Dr. J. W. Mauchly, Electronic Control Co., 1215 Walnut St., Philadelphia, Penn. 

Margaret W. Mayall, res. assist., Harvard Coll. Observatory, Cambridge 

D. H. Menzel, prof. astrophysics, Harvard Coll. Observatory, Cambridge 

Nicholas Metropolis, assist. prof., Univ. Chicago 

F. A. Metz, electronic engin., Watson Lab., AAF AMC, Cambridge 

J. W. Miles, physicist, Univ. Calif., Los Angeles 

F. G. Miller, electrical engin., Computation Lab., Harvard Univ. and Dahlgren Proving 
Gr., Dahlgren, Va. 

H. R. Mimno, prof. appl. physics, Harvard Univ. 

Richard von Mises, prof. aerodynamics and appl. math., Harvard Univ. 

Herbert Mitchell, grad. stu., Computation Lab., Harvard Univ. 

E. B. Mode, prof. math., Boston Univ., Boston 

R. J. Monroe, assist. director, Inst. Statistics, Univ. N. Carolina, Raleigh, N. C. 

C. N. Mooers, stu. reporter, Mass. Inst. Techn. 

Dr. B. L. Moore, res. fellow. Harvard Univ. 

Commander H. M. Mott-Smith, res. and develop. div., Navy Dept., Bur. Ord., Wash. 
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. C. F. Muckenhoupt, assist. chief sci. center, Office Naval Res., Boston 

G. Muller, engin., Bell Tel. Labs., Murray Hill, N. J. 

E. Nash, physicist, Underwater Sound Lab., New London, Conn. 

. B. Neal, project engin., Sperry Gyroscope Co., Great Neck., L. I. 

. H. Newell, prof. electr. engin., Worcester Polyt. Inst. 

. B. Nichols, science ed., Christian Sci. Monitor, Boston 

Natalie N. Nicholson, libr., Grad. Sch. Engin., Harvard Univ. 

r. E. N. Nilson, analytical engin., United Aircraft Corp., East Hartford, Conn. 
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~) 


Capt. K. H. Noble, assist. chief res. and development, Navy Dept., Bur. Ord., Wash. 


W. J. Nolan, Jr., res. assoc. electr. engin., Mass. Inst. Techn. 

Dr. Harry Nyquist, engin., Bell Tel. Lab., New York 

J. A. O’Brien, res. assist. electr. engin., Mass. Inst. Techn. 

J. W. O'Connor, office of vice-pres., John Hancock Mutual Life Ins. Co., Boston 
Otto Oldenberg, prof. physics, Harvard Univ. 

B. M. Oliver, engin., Bell Tel. Lab., New York 

R. D. O'Neal, physicist, Eastman Kodak Co., Rochester, N. Y. 

I. B. Parkhurst, assist. business manager, Harvard Coll. 

Dr. C. L. Pekeris, IBM Corp., Watson Lab., Columbia Univ. 

Lt. Comdr. J. T. Pendergrass, Office Chief Naval Operations, Navy Dept., Wash. 
C. G. Peterson, agent, Marchant Calculating Machine Co., Boston 

J. H. Phillips, British Supply Office, 1785 Mass. Av., N. W., Wash. 

Dr. E. R. Piore, electronic sec., Office Naval Res., War Dept., Wash. 

Dr. L. A. Pipes, res. physicist, Hughes Aircraft Co., Culver City, Cal. 

G. P. Ploetz, electronic engin., Watson Lab., AAF AMC, Cambridge 

W. A. Porter, constr. supervisor, Harvard Univ. 

H. A. Rademacher, prof. math., Univ. Penn. 

Dr. Jan Rajchman, res. physicist, Radio Corp. Amer. Labs., Princeton, N. J. 

E. A. Rasor, actuarial mathem., Social Security Admin., Wash. 

I. G. Recant, aeronautical engin., Nat. Adv. Comm. Aeronautics, Langley Field, Va. 
Dr. Mina S. Rees, chief, math. sect., Office Naval Res., Wash. 

D. K. Reynolds, electronic res. lab., Cruft Lab., Harvard Univ. 

Edward Reynolds, administrative vice-pres., Harvard Univ. 

George Reynolds, engin. aide, Watson Lab., AAF AMC, Cambridge 

C. H. Richards, electr. designer, Harvard Univ. 

Dr. F. D. Rigby, mathem., Navy Dept., Office Naval Res., Wash. 

R. F. Rinehart, prof. math., Case Sch. Appl. Sci., Cleveland, O. 

J. W. Roche, technician, Computation Lab., Harvard Univ. 

E. Rodriguez, stu., Harvard Coll. 

Commander J. B. Roman, material office for personnel, First Naval District, Boston 
J. H. Rosenbloom, physicist, Frankford Arsenal, Philadelphia 

Dr. Morris Rubinoff, instr. physics, Harvard Univ. 

Reinhold Rudenberg, prof. electr. engin., Harvard Univ. 

J. B. Russell, assoc. prof. electr. engin., Columbia Univ. 

E. A. Ryavec, mathem., Navy Dept., Office Naval Res., Wash. 

R. C. Sandborn, stu., Harvard Coll. 

B. L. Sarahan, mathem., Naval Res. Lab., Anacostia, Wash. 

Lt. Comdr. W. L. Savidge, distr. intelligence officer, First Naval District, Boston 
J. B. Scarborough, prof. math., U. S. Naval Acad., Annapolis, Md. 

A. C. Schaeffer, assoc. prof. math., Stanford Univ. 

F. F. Schaller, jun. engin., N. E. Power Service Co., Boston 

T. E. Schuette, project engin., Barber Colman Co., Rockford, III. 

R. E. Scott, res. assist. electr. engin., Mass. Inst. Techn. 

H. E. Seaton, electr. engin., The Autocall Co., Shelby, O. 

R. R. Seeber, Jr., res. assist., Watson Lab., IBM Corp., New York 
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Dr. R. J. Seeger, physicist, chief. mech. div., Naval Ordnance Lab., White Oaks, Silver 
Springs, Md. 

Dr. Harner Selvidge, director special products dev., Bendix Aviation Corp., Teterboro, N. J. 

Dr. C. E. Shannon, engin., Bell Tel. Lab., New York 

Dr. T. K. Sharpless, project super. EDVAC, Univ. Penn. 

C. B. Sheppard, Electronic Control Co., 1215 Walnut St., Philadelphia 

Martin Shwartz, mathem., Office Naval Res., Boston 

L. B. Slichter, prof. geophysics, Univ. Wisconsin 

A. E. Smith, mathem., Navy Dept., Office Naval Res., Washington 

Dr. C. L. V. Smith, senior engin., Raytheon Mfg. Co., Waltham, Mass. 

V. G. Smith, prof: electr. engin., Univ. Toronto 

Frances E. Snyder, mathem., ENIAC, Ballistic Res. Lab., Aberdeen Proving Gr., Md. 

Andrew Sobczyk, mathem., Watson Lab., AAF AMC, Cambridge 

R. A. Spence, engin., office bldgs. and grounds, Harvard Univ. 

Dr. R. C. Spencer, chief, antenna lab., Watson Lab., AAF AMC, Cambridge 

L. R. Steinhardt, director sp. projects, Engineering Res. Associates Inc., St. Paul, Minn. 

Dr. T. E. Sterne, chief ballistician, Ballistic Res. Lab., Aberdeen Proving Gr., Md. 

Dr. G. R. Stibitz, consultant appl. math., 393 South Prospect St., Burlington, Vt. 

Dr. W. W. Stifler, Jr., physicist, Underwater Sound Lab., New London, Conn. 

J. A. Stratton, director, Res. Lab. Electronics, Mass. Inst. Techn. 

Dr. L. P. Tabor, chief electronics sect., Moore Sch., Univ. Penn. 

G. R. Tatum, assoc. director res. appl. physics, The Johns Hopkins Univ. 

Margaret O. Taylor, mathem., Gulf Res. and Development Co., Pittsburgh, Pa. 

Richard Taylor, assist. prof. electr. engin., Mass. Inst. Techn. 

L. H. Thomas, physicist, IBM Corp., New York 

Lt. P. D. Thompson, AAF, res. assist., Institute Adv. Studies, Princeton, N. J. 

Capt. F. P. Tibbits, ordnance officer, Charlestown Naval Base, Boston 

P. D. Tilton, res. assist. electr. engin., Mass. Inst. Techn. 

Dr. C. B. Tompkins, engin., Engineering Res. Associates Inc., St. Paul, Minn. 

Dr. A. M. Turing, research fellow, National Res. Lab., Teddington, England 

Dr. A. W. Tyler, physicist, Eastman Kodak Co., Rochester, N. Y. 

C. F. Van Aken, chief machine tabul. div., Bureau of Census, Washington 

Capt. A. M. Van Eaton, commanding officer, Office Naval Res., Boston 

F. M. Verzuh, res. assoc. electr. engin., Mass. Inst. Techn. 

J. H. Wakelin, res. supervisor, Engineering Res. Associates Inc., St. Paul, Minn. 

Henry Wallman, assoc. prof. math., Mass. Inst. Techn. 

John Walsh, electr. engin., Watson Lab., AAF AMC, Cambridge 

J. H. Warden, president, The Autocall Co., Shelby, O. 

K. L. Warthman, engin., Bell Tel. Lab., New York 

Dr. A. T. Waterman, chief sci. planning div., Navy Dept., Office Naval Res., Wash. 

J. R. Weiner, engin., Raytheon Mfg. Co., Waltham, Mass. 

Albert Wertheimer, head engin., Navy Dept., Bur. Ord., Washington 

Dr. F. J. Weyl, mathem., Navy Dept., Bur. Ord., Wash. 

D. P. Wheatland, Cruft Lab., Harvard Univ. 

Stephen White, science news writer, New York Herald-Tribune 

Lt. Comdr. A. W. Whitney, aide to commandant, First Naval District, Boston 

Hassler Whitney, prof. math., Harvard Univ. 

D. V. Widder, prof. math., Harvard Univ. 

Maria Wilhelm, science sect., Newsweek Mag., New York 

R. E. Wilkins, mech. engin., Harvard Univ. 

S. B. Williams, cons. electr. engin., 366 Clermont Ave., Brooklyn, N. Y. 

Capt. J. H. Willingham, Jr., commandant, Naval Ammunition Depot, Hingham, Mass. 


M. K. Wood, techn. assist., progr. monitoring div., headquarters, Army Air Force, Wash. 


L. F. Woodruff, techn. chief sci. branch, W.D.G.S. 2, Pentagon Bldg., Wash. 
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R. L. Wooley, electronic engin., Watson Lab., AAF AMC, Cambridge 

Dr. A. W. Wundheiler, mathem., Navy Dept., Bur. Ord., Wash. 

Jeff Wylie, bureau head, Boston Office, Time Mag. 

Patrick Youtz, res. assoc. eleetr. engin., Mass. Inst. Techn. 

Dr. V. K. Zworykin, director electronic res., Radio Corp. Amer., Princeton, N. J. 


2 <<. 


1 Richard Henry Babbage (1883- ) is a son of Henry Whitmore B. (1855-1911), son 
of Henry Prevost B. (1824-1918), son of Charles Babbage (1792-1871). He is an assistant 
editor of Canada’s national farm magazine, Family Herald and Weekly Star. 


APPENDIX 


A. Existing and Practically completed Digital Computing Machines de- 
veloped by or presently on loan or contract to the Government. 


I. Department: War 
Branch or Bureau: Ordnance Dept. 
Permanent Location of Machine: Ballistic Research Leb. Aberdeen Proving Ground, Md. 
Type: A. Relay (small) manufactured by IBM Corp. 
B. Relay (small) manufactured by IBM Corp. 
C. Electronic (ENIAC) manufactured by Moore School of Electrical Engineer- 
ing, University of Pennsylvania. See MTAC, v. 2, p. 97f. 
D. Relay (large) manufactured by Bell Telephone Labs. 
Il. Department: War 
Branch or Bureau: Army Ground Forces 
Permanent Location of Machine: Army Ground Forces, no. 4, Fort Bliss, Texas. 
Type: Relay (small) manufactured by Bell Telephone Labs. 


III. Department: Navy 
Branch or Bureau: Office of Naval Research 
Permanent Location of Machine: Naval Research Lab., Bellevue, Md. 
Type: Relay (small) manufactured by Bell Telephone Labs. 
IV. Department: Navy 
Branch or Bureau: Bureau of Ordnance 
Permanent Location of Machine: Harvard Computation Lab., Cambridge, Mass. 
Type: Electromechanical (IBM Sequence Controlled Calculator = ‘‘ Mark I Calculator’’). 
See MTAC, v. 2, p. 185f. 
V. Department: Navy 
Branch or Bureau: Bureau of Ordnance 
Permanent Location of Machine: Navy Proving Ground, Dahlgren, Va. [In probable 
operation by July 1947]. 
Type: A. Relay (small) manufactured by IBM Corp. 
B. Relay (large) ‘Mark II Calculator,’”’ manufactured by the Computation Lab., 
Harvard Univ. [Epiroriat Note: It is about 12 times faster than Mark I, 
and cost about $400 000. ] 


VI. Department: National Advisory Committee for Aeronautics 


Permanent Location of Machine: Langley Memorial Lab., Langley Field, Va. 
Type: Relay, manufactured by Bell Telephone Labs. (under contract expiring in 1948). 


B. Present Automatic Digital Computing Machine Construction and 
Development Projects supported entirely or in part by the Government. 
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I. Department: War 
Branch or Bureau: Ordnance Department 
Description of activity: 
A. Moore School of Electrical Engineering, Construction of an electronic digital 
machine (““EDVAC’”’). 
B. Institute for Advanced Study and RCA Labs., Princeton, N. J. Construction 
of an electronic digital machine (financed only in part by federal funds). 
C. National Bureau of Standards, Washington, D. C., Long-range component 
development program. 
II. Department: Navy 
Branch or Bureau: Office of Naval Research 
Description of Activity: 
A. National Bureau .of Standards, Construction of an electronic digital machine. 
B. Servomechanisms Lab., Mass. Institute of Technology, Construction of an 
electronic digital machine to be used in a large guided-missile flight simulator. 
III. Department: Navy 
Branch or Bureau: Bureau of Ordnance 
Description of Activity: 
A. Harvard Computation Lab., Research and preparation of specifications for an 
electronic digital machine. 
B. Naval Ordnance Lab., White Oaks, Md., Construction of an electronic digital 
machine, temporarily abandoned (Dec. 1946). 
IV. Department: Commerce 
Branch or Bureau: Bureau of the Census 
Description of Activity: National Bureau of Standards, Construction of an electronic 
digital machine. 
J. H. Curtiss 
National Bureau of Standards 


EpiroriaL Norte: Other projected digital machines, unaided in their construction 
by the U. S. Government, are being built by: 
(i) The Eastman Kodak Co., Rochester, N. Y. 

(ii) The University of California, Berkeley. 

(iii) The National Physical Laboratory, Teddington, England, under the direction of Dr. 
ALAN M. Turinc. The planned Automatic Computing Engine will work at the speed 
of the ENIAC or possibly somewhat higher, and will take advantage of new technical 
developments, making possible both a greater memory capacity and a higher degree of 
complexity in the instructions. 


P. G. Scheutz, Publicist, Author, Scientific 
Mechanician, and Edvard Scheutz, En- 
gineer,— Biography and Bibliography 


PeuR GeEorG SCHEUTZ (1785-1873), son of a tavern-keeper, after passing 
examinations in law, practised in different localities before settling in Stock- 
holm in 1812. In 1817 as owner of a printing establishment founded by the 
well-known writer F. CEDERBORGH, Scheutz devoted himself from that time 
principally to literary interests. He soon became part owner and co-editor 
of a newspaper which, with changed name Argus (1820-1836), became 
Sweden’s most important political newspaper in the 1820’s. During 1826- 
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1842 he published five other journals and newspapers devoted to manu- 
factures and management; art, sloyd, and kindred sciences; engineering; 
industry and trade. In the field of literature he worked as a translator of 
Shakespeare, Werner, Walter Scott, Boccaccio, and others; his editions of 
Merchant of Venice (1820) and of Julius Caesar (second ed., 1831), the 
first translations of Shakespeare into Swedish, are in the Library of Harvard 
University. (His unpublished translation of King Lear was performed at 
Lindeberg’s Theatre). He also published numerous handbooks, and books 
of instruction. In 1842 he gave up his printing establishment and became 
one of the regular staff of the Aftonbladet (Evening News) Stockholm, with 
special interest in technical and economic subjects, and continued this 
interest until a few years before his death (Nos. 24, 25). 

But it was neither as a literary man nor as a publicist that Scheutz 
achieved his greatest commendation, but rather as inventor of a calculating 
machine. The Difference Engine of CHARLES BABBAGE (1792-1871) had 
from time to time been cursorily noticed in several periodicals, when a 
circumstantial and elaborate disquisition on its merits and construction 
appeared in the Edinburgh Review for July 1834. It was from perusal of this 
article that Scheutz, at that time the editor of a technological journal in 
Stockholm, derived the first conception of constructing a machine for effect- 
ing the same purpose as that of Babbage, namely of calculating and simul- 
taneously printing numerical tables. But after he had satisfied himself of 
the practicability of the scheme, by constructing various models, composed 
of wood, pasteboard, and wire, he postponed to a future period the further 
prosecution of the design (No. 10). Three years afterwards, in the summer 
of 1837, his son EpvarpD (1821-1881), then a student at the Royal Techno- 
logical Institute, Stockholm, took up the problem of constructing a working 
model in metal. Working with his father, many improvements in the original 
conception were introduced, and numerous alterations effected. By 1840 the 
apparatus was so far completed that it correctly calculated the value of 
series to 5D and one difference also of five figures. By April 1842 the model 
was extended so as to calculate similar series with two and three orders of 
differences. In 1843 the printing apparatus and all other parts of the model 
were in readiness for the inspection of the Royal Swedish Academy of 
Sciences. After several trials a certificate of the machine’s performance was 
given in Sept. 1843. Making use of this certificate as a recommendation the 
inventors sought for orders in various countries; but meeting with no success, 
the model was shut up in its case during the ensuing seven years (No. 10). 

In 1850 another inspection was made by a Committee of the Royal 
Academy of Sweden, and in 1851 Georg Scheutz made an application to the 
Government for the means to construct a large and still more improved 
machine. Such a grant was finally made. The new machine, Difference 
Engine no. 1, was completed in October 1853, being manufactured by C. W. 
Bergstrom of Stockholm. From the first it was found to work perfectly. 
During the latter part of 1854 the inventors visited England and France. 
The machine aroused considerable interest among a number of men of 
science. Upon being placed in the Great Exhibition at Paris, in 1855, the jury 
awarded a gold medal to the inventors. In February 1856 Georg Scheutz was 
made a member of the Swedish Academy of Sciences and in the same year 
created a Knight of the Order of Vasa and St. Anna. In 1858 he was also 
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created a Knight of the Order of the North Star. The Academy voted him 
an annuity of 1200 riksdollar' in 1860; and in 1872 the Academy awarded 
the Carl John Prize, to “Georg Scheutz, the first who successfully clothed 
Shakespeare in a Swedish costume, and for whom literature, even though 
as author, an occupation, comprehensive and pursued to the evening of a 
long life, has in addition for a long period had a connection which has not 
been effaced by the fact that the man of letters has also acquired a respected 
name for himself in a field which lies beyond the boundaries of belles lettres.”’ 
(No. 25). 

Through the interest of the astronomer B. A. GouLp, then Director of 
the Dudley Observatory at Albany, N. Y., Joun F. RATHBONE,? a manu- 
facturer of that city, purchased Difference Engine no. 1 for $5000, and 
presented it to the Observatory in 1856. There it remained unused until sold 
in 1924 to the Felt & Tarrant Manufacturing Co., Chicago, IIl., where it 
may still be seen, in a closed display case at the company’s museum at 
1735 North Paulina St., Chicago 22. 

Scheutz’s Difference Engine no. 2, with improvements, was in 1858 con- 
structed for the office of the Register General at Somerset House, by Messrs. 
Bryan Donkin & Co., from drawings of Edvard Scheutz. It contained some 
modifications of D. E. no. 1, and was us.d in the next few years for compu- 
tations indicated in Nos. 17 and 21 below. In 1914 the machine was pre- 
sented to the Science Museum, South Kensington where it is still exhibited 
and occasionally operated for visiting experts. 

The descriptions and illustrations in No. 31 are interesting, and the 
‘ capabilities of the machine are somewhat elaborately set forth in No. 18. 
Charles Babbage and his son Henry were enthusiastic in their praise of the 
achievements of the Swedish inventors (Nos. 4, 5, 9). In No. 9 Charles 
Babbage, who urged upon the Royal Society that their Engine was “highly 
deserving of a Medal,”’ added 

“The principle of Calculation by Differences is common to Mr. Scheutz’s 
engine and to my own, and is so obviously the only principle, at once 
extensive in its grasp and simple in its mechanical application, that I have 
little doubt it will be found to have been suggested by more than one 
antecedent writer. 

“Mr. Scheutz’s engine consists of two parts,—the Caiculating and the 
Printing; the former being again divided into two,—the Adding and the 
Carrying parts. 

“With respect to the Adding, its structure is entirely different from my 
own, nor does it even resemble any one of those in my drawings. 

“The very ingenious mechanism for carrying the tens is also quite 
different from my own. 

“The Printing part will, on inspection, be pronounced altogether unlike 
that represented in my drawings; which, it must also be remembered, were 
entirely unknown to Mr. Scheutz.” 

Among Edvard Scheutz’s own inventions was a rotary steam engine 
which was used on some steamboats. In 1860 the Swedish Academy awarded 
him a prize for this. (Nos. 24, 30; see also Poggendorff’s Biographisch- 
Literarisches Handwérterbuch, v. 32, Leipzig, 1898.) 


? A riksdollar is worth a little more than an American dollar. 


2 There is a brief biography of le FINLEY RATHBONE in A ppleton’s Cyclopedia of 
American Biography, New York, v. 5, 1888. 
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13. F. N. M. MorGno, Cosmos, v. 13, 1858, p. 78-84. 

[An account, with a full-page illustration, of the Scheutz machine, and reports, with 


comments, of discussions by Babinet, Leverrier, and others, at a meeting of the Académie 
des Sciences, Paris, 12 July 1858.] 
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°14. Three articles in Aftonbladet [Evening News ], 1859, nos. 14, 17, 253. 
[Listed in No. 25.] 

°15. Report on the Calculating Machine recently constructed by Mr. Donkin. 
London, 1859. 4to. 4 p. [Listed in No. 25. An account of Difference Engine 
no. 2. | 

16. G. ScHeuTz, ‘“Scheutz’s Raknemaskin,” Svenska Vetenskapsakad., 
Stockholm, Ofversigt af . . . Férhandlingar, v. 16, 1859, p. 391. 

°17. Mountain Barometer Tables: Calculated and Stereoglyphed by Messrs. 
Scheutz’s Calculating Machine no. 2, and printed by machinery, London, 
1859. 19 p. 

[A small booklet, 2 X 6} ins., “for private circulation,” printed for Wm. Gravatt 
F.R.S., foreign member of the Royal Academy of Sweden—the first publication of results 
from Difference Engine no. 2. There has been a copy in the Science Museum, South Kensing- 
ton, since 1914. Listed in Nos. 25 (with date 1859) and 31 (with date 1862). The title here 
seems to have been influenced by the criticism in No. 10 (ii).] 


18. H. MEIDLINGER, ‘‘Die Scheutz’sche Rechenmaschine,”’ Dinglers Poly- 
technisches Jn., Stuttgart and Augsburg, s. 4, v. 6, 1860, p. 241-256, 321-336. 

19. A. DeEMorGan, “Table,’’ English Cyclopaedia, Arts and Science 
Section, London, v. 7, 1861, col. 1007, 9 lines. 

[The date of publication of No. 17 is here given as 1859.] 


20. CHARLES BABBAGE, Passages from the Life of a Philosopher, London, 
1864, p. 47-48. 

21. English Life Tables. Tables of Lifetimes, Annuities, and Premiums. 
With an Introduction by William Farr, M.D., F.R.S., D.C.L. Published by 
Authority of the Registrar General of Births, Deaths, and Marriages in England. 
London, printed for Her Majesty’s Stationery Office, 1864. clvi, 606 p. 
“English Life Table no. 3’’ p. 6-605. 

[Appendix, ‘‘Scheutz’s calculating machine and its use in the construction of the 
English Life Table no. 3,” p. cxxxix—cxliv. “‘This volume is the result; and thus—if I may 
use the expression—the soul of the machine is exhibited in a series of Tables which are 
submitted to the criticism of the consummate judges of this kind of work in England and 
in the world,” p. cxl. The tables calculated and printed by Difference Engine no. 2 are on 
p- 6-11, 42-47, 76-81, 114-115, 129-133, 142-145, 175-245 (mostly), 355-425 (mostly), 
481-551 (mostly). There are copies of these tables at Harvard University, and in the 
South Kensington Museum. ] 


22. ““The Swedish tabulating machine of G. & E. Scheutz,” Annals of 
the Dudley Observatory, v. 1, Albany, 1866, p. 116-126 +- large folding plate 
of the machine on a stand. 

23. J. W. L. GviatsHer, “Report of the Committee on mathematical 
tables,” reprint from B.A.A.S., Report 1873, p. 59, 159. 6 lines about No. 10. 

24. HERMAN HoOFBERG, Svensk Biografiskt Handlexicon, Stockholm, v. 2, 
1876, p. 242-243. 

[Sketches of P. G. and E. Scheutz. Copy at Library of Congress. ] 

New edition, Stockholm, v. 2, 1906, p. 428-429: 


[Biographical sketches, with portraits, of Per Georg Scheutz and Edvard Scheutz. Copy 
at Harvard University. This, No. 25, No. 2, and No. 10, were the main bases for the compila- 
tion of our sketch. ] 


25. C. F. Bercstept, ‘Georg Scheutz,’’ Svenska Vetenskapsakad., 
Lefnadsteckningar, v. 2, part 1, Stockholm, 1878, p. 155-179. 








244 P. G. SCHEUTZ AND EDVARD SCHEUTZ 


[The most extensive biography of Scheutz, with a complete bibliography of his pub- 
lications. ] 

26. MauRIcE D’OcaGNnE, Le Calcul Simplifié par les Procédés Mécaniques 
et Graphiques, Paris, 1894, p. 37-38. 

Second ed., Paris, 1905, p. 83-86, 213-218. 

Third ed., Paris, 1928, p. 74-77, 182-187 (‘‘Note sur la machine a 
différences Scheutz’”’ by Colonel BERTRAND). 

[In the third edition, there is a picture of the machine on p. 75. On p. 76 it is incorrectly 
stated that after Difference Engine no. 1 was at Albany it was utilized for the calculation 
of tables of logarithms, sines, and logarithmic sines (see No. 9). The statement concerning 
Difference Engine no. 2, ‘‘Cette machine a calculé et imprimé 605 tables (grand in-4°) qui 
constituent le fondement du calcul des rentes viagéres servies par les caisses d’épargne 
postale anglaises”’ is also highly inaccurate and misleading; see No. 21.] 


27. V.G. von Boo., Pribory i mashiny dita mekhanicheskago Proizvodstva 
arifmeticheskikh Deistvit [Apparatuses and machines for mechanical execu- 
tion of arithmetical operations ]. Moscow, 1896, p. 184-189. 


[There is a copy of this work in Library of Congress, by VLADIMIR GEORGIEVICH VON 
Boot (1835-1889). A picture of the Scheutz Difference Engine no. 1 fills p. 1893. 


28. R. MEHMKE & M. p’OcaGneE, ‘“‘Calculs numériques,” Encycl. d. Sci. 
Math., I, 4. 2. Paris, 1908, p. 268-269. 


[Mostly accurate and with good bibliography; there are a number of errors in the 
reference to no. 27.] 


29. L. Jacos, Le Calcul Mécanique, Paris, 1911, p. 115-117. [Seems to be 
an abridgment of No. 26, second ed. ] 

30. Kungl. Svenska Vetenskapsak., Personforteckningar, 1739-1915, ed. 
E. W. DABLGREN, Stockholm, 1915. 

31. D. BAXANDALL, Mathematics. I. Calculating Machines and Instru- 
ments (Catalogue of the Collections in The Science Museum, South Kensington 
with Descriptive and Historical Notes and Illustrations.) London, 1926, p. 32, 
34-36 plate IX (i. Scheutz’s Difference Engine no. 2; ii. Details of wheelwork; 
iii. Portion of a calculated and printed table). 

[The entries of the Catalogue are: no. 86, portrait of Edvard Scheutz during his student 
days at the Royal Technological Institute, Stockholm. No. 87, no. 7. No. 88, photograph 
of the Scheutz Difference Engine no. 1 at the Dudley Observatory. No. 89, No. 10. No. 90, 
Scheutz Difference Engine no. 2, presented by the General Register Office, Somerset 
House, 1914. No. 91, No. 17, but with date 1862. No. 92, No. 21. No. 93, moulds and 
stereotypes produced by the Scheutz Difference Engine no. 2. No. 94, Stereotype block 
from a matrix produced by the Scheutz Difference Engine no. 2. No. 95, Lead mould, 
papier-maché moulds (3) and stereotype produced by the Scheutz Difference Engine no. 2.] 


32. Svensk Uppslagsbok, v. 24, Malmé, 1937. 
[Sketch and portrait of P. G. Scheutz, col. 119]. 


33. H. H. AIKEN, etc., A Manual of Operations for the Automatic Sequence 
Controlled Calculator (Annals of the Computation Laboratory of Harvard 
University, v. 1). Cambridge, Mass., 1946, p. 6-7. 

(Quotation: “In 1834 George Scheutz, a printer in Stockholm, built a less ambitious 
difference engine with the aid of a grant from the Swedish government. This machine was 
completed in 1853 and used for the computation and printing of tables of logarithms, sines, 
and logarithms of sines.” In view of the details given above the first sentence clearly conveys 
more than one erroneous impression. The latter sentence is equally misleading since only 
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the tables described in No. 10 are in question. Because the reviewer in MTAC, v. 2, p. 186, 
‘ followed the Harvard Annals v. 1, he listed ‘““Scneutz (1834)”, rather than Sceutz (1853).] 


R.C.A. 


A New Approximation to z 


A, EpitoriaL Notes: In MTAC, v. 2, p. 143-145 we noted various 
formulae which had been used for calculating + to many places of decimals. 
These included that of MAcuINn (1706) 


(1) = = 4tan} — tan gts, 


which was used by WILLIAM SHANKS (1812-1882) to compute z to 707D. 
The accuracy of this computation to 500D was verified by an independent 
| calculation completed and published in 1854. No one appears to have 
| checked the later figures until 1945, when Mr. D. F. FERGUSON, now con- 
| nected with the Department of Mathematics of the University of Man- 
chester, undertook the task. As we have already noted he used the formula 


(2) = = 3tan} + tan gy + tan” ress, 


| given in LonNeEy’s Plane Trigonometry (1893). We published his new compu- 
tation for + from the 526th through the 620th decimal place (p. 145). Mr. 
Ferguson gave an account of his work in Mathematical Gasette, v. 30, May 
1946, p. 89-90, and recorded there his figures for x from the 521st to the 
540th decimal place. Mr. Ferguson found that Shanks’ approximation to + 
| was incorrect beyond 527D. By November 1946 he had carried on his 
calculations of the value of x to 700D, and by January 1947 to 710D. 

In December 1945 we suggested to Dr. Jonn W. WRENCH, JR. that he 
might take up the wholly independent computation of + by means of 
Machin’s formula (1). In April 1946 he reported that he was in communica- 
tion with Mr. Levi B. Smita of Talbotton, Georgia, who began his work on 
computing tan-! g$s in November 1940 and had by February 1944 com- 
| pleted the work to 820D, through the term [173-239!%]~. Then Dr. W. 
took up actively the computation of tan~' } so that his results might be 
combined with those of Mr. S. as in Machin’s formula. He found the errors 
in work of Shanks, earlier pointed out by Mr. F., and others described below. 
| Early in January 1947 Dr. W. sent to us his new approximation to z to 





808D given below, as a companion to the value of e to 808D (MTAC, v. 2, 
April 1946, p. 69). The value found by Mr. F. to 710D agrees with this. 


mw = 3.14159 26535 89793 23846 26433 83279 50288 41971 69399 37510 
58209 74944 59230 78164 06286 20899 86280 34825 34211 70679 
82148 08651 32823 06647 09384 46095 50582 23172 53594 08128 
48111 74502 84102 70193 85211 05559 64462 29489 54930 38196 
44288 10975 66593 34461 28475 64823 37867 83165 27120 19091 
45648 56692 34603 48610 45432 66482 13393 60726 02491 41273 
72458 70066 06315 58817 48815 20920 96282 92540 91715 36436 
78925 90360 01133 05305 48820 46652 13841 46951 94151 16094 
33057 27036 57595 91953 09218 61173 81932 61179 31051 18548 
07446 23799 62749 56735 18857 52724 89122 79381 83011 94912 
98336 73362 44065 66430 86021 39494 63952 24737 19070 21798 
60943 70277 05392 17176 29317 67523 84674 81846 76694 05132 
00056 81271 45263 56082 77857 71342 75778 96091 73637 17872 
14684 40901 22495 34301 46549 58537 10507 92279 68925 89235 
42019 95611 21290 21960 86355 44191 19716 02977 46113 ner 
1 
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B. REPORT OF Mr. SmitH & Dr. WRENCH, January 1947 


As stated in the Editorial Notes formula (1) was employed by the 
present writers in their: joint calculation of x. The calculation of tan- 4 
consisted essentially of checking and extending to 850D the values of the 
individual terms of the corresponding series as published to 530D by Shanks.! 
On the other hand, the computation of tan- x} was carried out to 820D 
entirely independently of earlier calculations of that number, and the re- 
sulting values of the terms of the series were checked against the correspond- 
ing data of Shanks only when the investigation was nearly completed. 

An important preliminary step in the new evaluation of tan—! } consisted 
of the formation of a definitive table of powers of 2. This original table 
contains the exact values of 2" for m = 1(2)1207, and has been collated with 
an unpublished table of non-consecutive powers to 2°! computed by Pro- 
fessor H. S. Uhler® and also with a table of 2", m = 13(12)721, given by 
Shanks.! No discrepancies were found. In addition to this comparison with 
previous tables of high powers of 2, every entry beyond 2! in the new table 
was checked by the Fermat-Euler theorem. 

The quotient arising from the application of this congruential check on 
the accuracy of the tabular value of 2?"-' comprised, together with the 
appropriate number of antecedent zeros, the sequence of digits occupying 
the first 2n — 1 decimal places of the approximation to the mth term of the 
series for tan 3. If r denotes the residue determined by the preceding check, 
then the decimal evaluation of r/(2n — 1) is also required, corresponding to 
all integral m between 1 and 425. Thus it was found desirable to calculate 
de novo a table of the complete periods of the reciprocals of all prime-powers 
(p',k = 1, p ¥ 2,5) less than 800. In each of the many cases where the 
period consisted of an even number of digits an effective check involved the 
juxtaposition and subsequent addition of the two halves of the period so as 
to yield an unbroken sequence of 9’s.* All remaining cases of evaluation of 
r/(2n — 1) were checked by duplicate machine calculation, as was the 
summation of the terms of the series. 

The calculation of the successive terms of the series for tan-! s$y was 
performed by the recurrence formula 


Un41 (2n 7. 1) U,/(2n + 1)239", 


where U, and U,+1 denote respectively the mth term and its successor. The 
computations were carried to at least 820D and were checked modulo 10” + 1 
every hundred decimal places. The same checking procedure was applied 
to the respective sums of the positive and negative terms of the series. 

The final values of tan 3 and tan~ 335 were compared with the corre- 
sponding data to 709D of Shanks, and several errata in the latter were dis- 
covered. In addition to the two errors (described in C) which were inde- 
pendently discovered by Mr. Ferguson in Shanks’ value of tan-! }, there 
exist in the same number a unit error in the 533rd place and a residual error 
of approximately 5.2193762669 < 10-®. Shanks’ approximation to tan 345 
is also erroneous, for it exceeds the present estimate of that number by 
nearly 4.77447473 10-52. 

Appended to this report are values of (a) tan-'}, and (b) tan— gz, 
both curtailed to 811D from more extended approximations appearing on 
the work sheets. 
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(a) 
0.19739 
21015 
12628 
30327 
92806 
11323 
81553 
08560 
02324 


(b) 

0.00418 
19508 
54975 
09280 
60154 
58881 
83097 
14512 
51034 
07086 
14237 
34496 
02023 
27821 
78092 
57415 
28586 


40760 
27019 
45066 
08462 
54814 
27339 
70020 
32220 
35262 
31401 
67708 
44014 
49475 
13089 
05326 
88258 
38373 


49880 
94024 
36913 
10747 
68061 
21796 
29057 
73686 
49211 
96139 
24383 
87633 
46026 
78720 
41703 
66648 


02074 
61288 
64448 
71962 
80599 
78533 
64652 
94463 
82445 
93869 
36515 
24187 
97790 
84258 
86491 
29094 
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72386 
71817 
63756 
33077 
24498 
46518 
82763 
68627 
12693 
51950 
69182 
24950 
97469 
30836 
62667 
05850 


45382 
78341 
05245 
33759 
30211 
04574 
46532 
66155 
55670 
37058 
52702 
92283 
50111 
77190 
40271 
90382 


65194 
24378 
98867 


59311 
10525 
19285 
08123 
78815 
25572 
83775 
32030 
07922 


13623 
79954 


14959 
42289 
83947 
47634 
46039 
25481 
96910 
20831 
49968 
64107 
00229 
86239 
88854 
91008 
68444 
07331 


28545 
32737 
89311 
60331 
12539 
35867 
48183 
67964 
44452 
70855 
30895 
54553 
76673 
39098 
24477 
75908 


34475 
97828 
35574 
85526 
61329 


64919 
15800 
90799 
13846 
87783 
51006 
67995 
37431 
58081 
43132 


27410 
82605 
86505 
84734 
49527 
46447 
86543 
26465 
47904 
85540 
44950 
33565 
97957 
51655 
31579 
43199 


48065 
81362 
89221 
14570 
60779 
51979 
56078 
74655 
33177 
45223 
04358 
11719 
31537 
10419 
64520 
77843 


78785 


95216 
73186 
93397 


94017 
22514 
02007 
13529 
52464 
31257 


67431 
13759 
83892 


30763 
29094 
28833 
33198 
68815 
10232 
91959 
11032 
28393 
55388 
54409 
73747 
09303 
22416 
27549 
27604 


1W. SHANKs, Contributions to Mathematics comprising chiefly the Rectification of the 


Circle to 607 Places of Decimals, London, 1853. 
2M 224, N66. 


AC, v. 2, p. 
3H. RADEMACHER & O. ToEpLitz, Von Zahlen und Figuren, second ed. Berlin, 1933, 


p. 133-135. 


C. Report oF Mr. FERGUSON, January 1947 


In calculating my value of x to 710D I made use of formula (2) and 
gp, and (c) tan sys: 


found the following results for (a) tan! 3, (b) tan 


(a) 

0.24497 
40671 
63834 
20373 
87873 
35238 
02760 
13351 
37795 
50590 
10855 
23900 
39341 
42486 
81885 


86631 
27375 
80449 
$2012 
06750 
13507 
62335 
40849 
44423 
15074 
05491 
98621 
51460 
09165 
50930 


26864 
91466 
00371 
77087 


15417 
73551 
18374 
38758 
56233 
55568 
36107 


20824 
19587 
29548 
16557 


81211 
64209 
54209 
21586 
63956 
58302 
78574 
46394 
78861 
59686 
80499 
50318 
95391 
62989 


27581 
65745 
95059 
71598 
36978 
66262 
13846 
26115 
19777 
03414 
01695 
46160 
27283 
80883 


09141 
34157 
97695 
26385 
08521 
33025 
85151 
26552 
28818 
25179 
01187 
29411 
64199 
85199 


44098 
66870 
89869 
50632 
59107 
99157 
60692 
60206 
69907 
53473 
90130 
55558 
28418 
09565 


38118 
19913 
60614 
05220 
32458 
53281 
64028 
13960 
30107 
22290 
24132 
05517 
54722 
57220 
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0.04995 83957 21942 76141 00062 87034 84488 14912 77080 42350 
71744 10853 45482 99835 95476 71033 50612 64888 70485 01265 
49675 88718 56799 74803 45043 78235 17343 64195 86075 35558 
34705 50031 66812 64425 55070 35889 99864 21844 62020 22011 
35398 44491 94479 55125 91884 70605 15358 82203 57911 15507 
66709 70265 20884 04697 53559 08904 34425 33211 75071 00898 
99983 99369 89611 53196 70717 40134 40774 24235 31335 37603 
73612 47259 31779 72222 52596 59464 82850 02739 09656 29682 
61838 28530 42311 66214 89812 84597 81323 80425 73403 65277 
35640 82643 15372 91283 73850 56089 49548 56557 20164 84879 
81610 72192 83012 94406 89240 40051 11637 64820 56557 65999 
47240 24101 35373 55511 99718 11544 33853 54021 44594 33781 
36222 03768 16540 61055 38956 20032 50668 29159 05403 66710 
93525 58744 35937 48968 67734 87127 37233 28015 36651 95672 
97735 23477 (67) 


0.00050 37782 94913 08568 94071 15151 20340 68155 82974 27671 
70794 50754 94924 49051 47331 91562 48721 41344 62457 31663 
54356 72097 25292 46735 43656 70658 85121 98509 57714 34631 
16201 87555 72950 86848 41560 75739 76456 56371 24816 21875 
12054 38001 54144 09788 49786 28731 58173 88022 96546 91890 
13988 03384 98768 15748 37461 32808 20136 07891 78079 24576 
34848 81139 81141 28185 10421 04373 41755 93445 09515 54421 
06064 77781 30180 52296 70643 13015 42264 54341 08263 07459 
83039 69957 29909 17547 54316 44112 04827 48412 99637 24037 
14450 28084 07818 72581 81602 03033 62489 37749 65168 46978 
10408 29672 64677 37415 73398 53325 49265 37800 02819 17053 
46292 72603 22836 58266 41037 79381 23834 28205 57905 00949 
45767 62167 06957 55242 02247 36629 36425 52266 02749 98589 
82687 23984 76364 21164 11631 63537 47742 14457 26882 79973 
42113 22633 (35) 


(6) + (c) = tan—(5/99), which was independently computed and fur- 
nished complete agreement to 710D+. 

My procedure for calculating tan—! } was as follows: 

(i). calculated (1/4)?"* by dividing (1/4)?"—! by 16; 
(ii). multiplied the result by 4096 and compared with (1/4)?"-5; 
(iii). divided (1/4)?! by (2n + 1); 
(iv). multiplied this by 16(2” + 1) and compared with (1/4)?"—'; 
(v). after copying (1/4)"*"/(2m + 1) for the purpose of the series, checked 
by multiplying the copied figures by 2” + 1. 
At all steps of the work ample margins of overlap were allowed. 

In the course of my work I discovered two errors in the results of Shanks. 
The explanation of the first of these which vitiates his final result beyond 
527D was noted in January 1946 and is as follows: It was a question of an 
omission in the evaluation of the term [497-5“7]— in the 531st decimal 
place. I found the value to be (through 547D) 


00804 82897 38430 58350, 
while Shanks, carelessly omitting a zero, used 
00848 28973 84305 83501 


The second Shanks’ error was the omission from 569D+ of the term 
5-*/29 which comes in the series for tan~ }. 
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RECENT MATHEMATICAL TABLES 


359[A].—A. L. CrELLE, Rechentafeln welche alles Multiplizieren und Dividieren 
mit Zahlen unter Tausend ganz ersparen, bei grésseren Zahlen aber die 
Rechnung erleichtern und sicherer machen. Neue Ausgabe besorgt von O. 
Seeliger [1907]. Neudruck, mit Tafeln der Quadrat- und Kubikzahlen von 
11-1000. Berlin, Gruyter, 1944. viii, 501 p. 24.8 K 36.7 cm. See also 
MTAC, v. 2, p. 179. 


A. L. CRELLE (1780-1855) was the founder (1826) of a notable mathematical research 
journal still in existence, and author of various volumes, including even one on music (1823). 
His calculating tables, of which there have been many editions in English, French, and 
German, was first published in two small volumes over 125 years ago, Berlin, 1820. But in 
such editions, since 1857 at least, the folio format (each page containing four pages of the 
first edition) has been in use up to the present. As the title indicates the war edition under 
review contains two extra pages with the squares and cubes of numbers N, N = 1(1)999. 
A one-volume Japanese edition of Crelle’s tables by TsuNETA YANO, Tokyo, 1913, has been 
referred to, MTAC, v. 2, p. 18; our incidental note, v. 1, p. 436, that this was published by 
an insurance company, is incorrect. This fact was learned after Mr. Epwin G. BEAL, Jr., 
Chief of the Japanese Division in Library of Congress had kindly made a study of the volume 
for us. He also reported that he could find no record of a 1927 Japanese edition of this work. 


RC A. 


360[A, B, F].—J. Ser, La Numération et le Calcul des Nombres, Paris, 
Gauthier-Villars, 1944, 194 p., 25 K 16 cm. 


The author of this work appears to be something of an individualist. He gives but one 
reference to the work of others and this is to an article on the foundations of mathematics, 
a subject not covered in this book. The reader will find many new points of view, “‘méthodes 
personelles,” and unfamiliar nomenclature. Much of this has to do with pencil and paper 
calculation, a discipline all but unknown in this age of mechanized computation. 

The more extensive tables in this work may be described as follows: 


(i) An arithmetical table (p. 60-95) giving for the first 1000 integers N the following 
functions: 

1/N to 6D, the first nine multiples of N (these are used to facilitate multiplication 
and division), N?, VN to 4D, V10N to 4D, N®, N? to 4D, (10)! to 3D, and (100) to 3D. 

(ii) A small table, p. 96, of powers n* for nm = 1(1)9, k = 1(1)16, the first nine multiples 
of 10, 10!, 100% (all but five values have last-digit errata) and the fourth roots of m-10*, 
m = 1(1)9, k = 1, 2,3 to 4D. 

(iii) A table of the roots x of the linear congruence 


Rx = D (mod B) 


where R and B range over the first 50 integers and D is the greatest common divisor of R 
and B (p. 100-101). The reader will find that values of R are given as column headings. 

(iv) A table of the residues with respect to each of the moduli 2, 3, 5, and 7 of the first 
210 integers arranged in two ways (p. 122-123). 

(v) A factor table for each of the first 1000 integers except multiples of 10. To save 
space composite numbers are usually broken into only two factors, one of which is often a 
composite number (p. 144-145). 

(vi) A factor table for those numbers between 1000 and 10000 which are prime to 30. 
Beyond 4020 only the least prime factor is given if the number is composite (p. 146-149). 

(vii) A condensed factor table for numbers under 210000, not divisible by 2, 3, 5, or 7. 
This 12-page table (p. 150-161) is reminiscent of an unfinished project of E. LEBon,' and 
is too complicated to describe here in detail. The reader is warned to study directions before 
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attempting to use the table. For composite numbers whose least factor exceeds 210, the 
table yields the two factors at one “coup d’oeil,’”’ that is, after a little hunting. For other 
numbers some mental calculation involving the factoring of 10 or sometimes 25 three-digit 
numbers is necessary. 
DB &.. 
1E. Lepon, Table de Caractéristiques de Base 30 030 donnant, en un seul Coup d’Oeil, 


les Facteurs Premiers des Nombres Premiers avec 30 030 et Inférieurs & 901 800 900, v. 1, 
pt. 1, Paris, 1920. 


361[B].—ALBERT GLODEN, Table des Bicarrés X* pour 1000 < X < 3000, 
Luxembourg, author, rue Jean Jaurés 11, 1946. Offset printing on one 
side of each of 17 leaves, with paper cover. 20.5 X 29.7 cm. 


There is no text. The values of X* for X = 1001(1)1099 check with the values given 
in BAASMTC, Mathematical Tables, v. 9, 1940, p. 122-123. Some of the printing is unclear 
so that a “3,” one case tested, X = 1003, might easily be mistaken for a ‘‘5.”” This table 
was made preparatory to writing the paper reviewed in RMT 348. 

R. C. A. 


362[D].—Josrer KRovAk, Natiirliche Zahlen der Funktion Cotangens fiir 
Winkel in Zentesimalteilung von 0° bis 100°". Prague, Landesvermessungs- 
amt Béhmen und Mahren, second ed., 1943. iii, 396 p. 15.5 X 21.5 cm. 


The one-page explanation (dated Prague, 1943) of these tables says that they had 
already been announced in the Sechsstellige Tafeln der natiirlichen Werte der Funktionen 
Sinus und Cosinus fiir Winkel in Zentesimalteilung that the Finance Minister had published 
in the previous year. 

The tables were brought into being for use with trigonometrical survey calculations 
with twin calculating machines where the accuracy of measurement is of the order of 
2” (0“.6) or slightly better. In other words, it provides six significant figure values for 
working to 09.0001 or 1“ or one centesimal second or about one third of a sexagesimal second. 
They are a further outcome of Hitler’s decree that German surveyors were to use the 
centesimal division of the quadrant. 

The principal survey problem that is facilitated by cotangents and twin machines is 
that of intersection, i.e. the determination of the co-ordinates of a point whose bearings 
from two known points have been measured. 

The lay-out is shown by the following table: 


Pages From To Interval Diff. for 1‘* 
1- 20 Oe 19 ~ None 
21-160 1 15 2 63.6 to 1.5 
161-230 15 50 10 2.9to 0.3 
231-316 50 93 10 3.2 to 1.5 
317-376 93 99 2 1.5 to 15.9 
377-396 99 100 1 None 


Where the interval is 2“‘, mean differences are given for 1‘; where the interval is 10‘, com- 
plete proportional parts for centesimal seconds are given. 

The tables have been typed and reproduced from photographic plates. The Bremiker 
division of the lines has been used throughout. No description of the source of the values 
is given. 

It seems hard to justify the use of six significant figures throughout if the accuracy of 
measurement is limited to about 1’. For small angles, the last one, two or even three figures 
are meaningless. The same is true of angles near 100°, where the number of decimals in- 
creases steadily to 11. These extra decimals cannot be of any use in survey work, and are 
only likely to be a source of confusion. 

‘It is a blemish on the arrangement of the part of the table that is at interval 10“ (156 
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pages) that the centesimal minutes 0 to 50 are on a right-hand page and 50 to 100 on the 
following left-hand page; they should, of course, have been printed at a single opening. 
L. J. 


363[D, E, L].—MIk és Imre HETENYI, Beams on Elastic Foundation. Theory 
with Applications in the Fields of Civil and Mechanical Engineering. 
(Univ. Michigan Studies. Scientific Series, v. 16). Ann Arbor, Univ. of 
Michigan Press, 1946, p. 217-255. 17 K 25.2 cm. $4.50. 


The tables include (p. 217-239) graphs and 4D tables A, = e~* (cosx + sin x), 
B, =e*sinx, C, = e*(cosx —sinx), D, = e*cosx, for x = 0(.001).02(.01)4(.1)8, 
$x(4m)$r. There are also 5—-7S tables with graphs (p. 241-243) of Ei = }e*(sinh x + sin x), 
F, = }e*(cosh x + cosx), Ei, = $e7(sinh x — sinx), Fy = }e7(cosh x — cos x), for 
= 0(.05)3(.1)5. 
There are also graphs and 4D tables (p. 245-255) of Z:(x) = ber x, Z2(x) = — bei x, 
Zi’ (x), Ze’ (x), Z3(x), Za(x), Z3’(x), Za'(x) for x = 0(.01)6, where 
Zs = 4 ber x — (2/)[Ri — bei x(y + In }x)], 
Zs = — 4 beix + (2/x)[R2 + ber x(y + In 3x)], and 
“2 - 


= (fx)? — (3x)* + => (4x)" — 
2 4 6 ; 
R2'= “ (3x)* — _*e (3x)* + ae “ (jx)#—---, on) =) 1/k, 
I 
y = .577216--- 


364[D, P].—Istiruto GreoGraFico Miuttare, Florence, Tavole per Calco- 
lare le Differenze di Livello nelle Levate Topografiche e per Calcolare le 
Distanze ridotte all’ Orizzonte. (Collezione di Testi Tecnici). Florence, 1943, 
viii, 195 p. 19.3 X 23.6 cm. Full cloth. The Preface is signed by Prof. 
GIOVANNI BoaGa, geodetic chief. 


L is the distance AB of an object, a the angle of its elevation or depression with reference 
to the horizontal plane, D = AC the projection of L on this plane, L’ = AE the projection 
of AC = D, on L, and BC’= h. 


Table I, p. 3: D = Lcosa to 4D, for a = 1°(1°)30°, Z = 1(1)9. 

Table II, p. 14-41: L’ = Dcosa = Lcos*a, to 4 or 5S, for a = 0(5’)5°(2’)11°(1’)19°59’ 
and L = 1(1)9; also for L = 1, a = 20°(1’)45°. 

Table III, p. 43-135: h = D tana, to 5D, for D = 1(1)9, a = 0(15%)15°(30’)20°(1’)45°20’. 

Table IV, p. 191: corrections due to sphericity and refraction, differences of level in meters 
1000(100)25900, coefficient of refraction = .06733. 

Table V, p. 195: tan?a, for a = 0(50‘)309 50‘, Tables for correction of sphericity and 
refraction. 


The previous edition of this work appeared in 1915 (15.7 X 22.4 cm., 53 p.) and con- 
tained four tables. The first, and last two tables are practically equivalent to the first, fourth, 
and fifth tables of the 1943 edition. T. II (1915) is of the same plan as T. III (1943), but 
for the range [0(1’)45°; 4D]. This 1915 edition was an enlarged and corrected edition of a 
previously revised and corrected edition, which appeared in 1896. 

Rc. A, 


365[D, S].—Louts CouFFIGNAL, Tables de Produits de Lignes Trigono- 
métriques. Paris, Gauthier-Villars, 1943. iii p. + 24 thick paper leaves, 
printed on only one side. 31.3 X 23.5 cm. Boards, 210 francs. 


This volume was prepared under the’direction of Dr. CourFiGNaL, the director of 
the laboratory of mechanical calculation in the Centre National de la Recherche Scientifique. 
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The author’s volume, Les Machines @ Calculer. Leur Principes. Leur Evolutions. (Paris, 
1933, ix, 86 p.) is well known. His doctoral dissertation at Paris was entitled Sur Il’ Analyse 
Mécanique. Application aux Machines a Calculer et aux Calculs de la Mécanique Céleste 
(Paris, 1938, 132, 3 p. 4to). f 

On the back of the title-page of the present volume is a brief preface in French, German 
and English, and on the opposite page, again in three languages, are ‘‘Directions for use of 
the Tables.” We are told that ‘‘The establishment of crystal structure from X-ray diagrams 
demands extensive calculations, where products of two or three cosines occur continually. 
On the request of several French crystallographers the French National Office of Scientific 
Research has undertaken to publish tables which might facilitate this kind of work. Besides, 
such Tables may be useful in a great many cases of harmonic analysis.” 

The tables give the products P = f(X)-g(Y)-h(Z), where f, g, # are either sine or 
cosine functions. The arguments X, Y, Z are at interval one hundredth of a circumference, 
that is, 49 = 3°.6. From the table one may read off at once the value of P for any X, Y, Z 
in 49 units up to 100. For example, to evaluate P = sin 41-cos 65-cos7 first turn to 
cos Z = cos 7, p. 7 (in the upper right-hand corner of the page). Then on that page columns 
sin X = sin 41, cos Y = cos 65, indicate that P = — .2850. The results are all to 4D. 
The author states that the error in any P is less than 5-10. 


R. C. A. 


366[F].—A. GLopEN, “‘Compléments aux tables de factorisations de CuN- 
NINGHAM,” Mathesis, v. 55, 1946, p. 254-256. 16.2 X 25 cm. 


The tables referred to are those in which CUNNINGHAM gives! (with many incomplete 
entries) the factors of numbers of the form x* + 1 for x < 1000. The results quoted in this 
note serve to complete all but 51 entries in this table. Previous addenda by KRattcHiK? 
and BEEGER? are given and have been verified. The new. results are by-products of tables 
of the solutions of the congruence 


x* = — 1 (mod p) 


for p < 500000 by GLopEN and DELFELD.‘ Those values of x for which (x* + 1)/d is a 
prime between 10! and 25-10" are listed for d = 1, 2, 17, 34, 41, and 82. For some reason 
the author has failed to list 565 for d = 1 and 640, 648 for d - 2. Six other factorizations 
are given for x = 595, 598, 685, 714, 844, 880. The author has recently given a similar table® 
to Cunningham's for 1000 < x < 3000. The present note closes with a table of the factors 
of x® + 1 for x = 37, 41, 50, 52, 63, 82, 85, and 87. 
3 D. H. L. 

1A, J. C. CunNiNGHAM, Binomial Factorisations, v. 1, London, 1923, p. 113-119. 

2M. Krarrcuik, Recherches sur la Théorie des Nombres, v. 2, Paris, 1929, p. 116-117. 

3N. G. W. H. BEEGER, Additions and Corrections to Binomial Factorisations by Cunning- 
ham. Amsterdam, 1933, 1945. 

‘See MTAC, v. 1, p. 6; v. 2, p. 71-2, 210-211. 

5See MTAC, v. 2, p. 211. 


367[F].— MIKHAIL BoriIsovicH OsTROGRADSKII (1801-1861) Polnoe Sobranie 
Sochinenit Akademika M. B. Ostrogradskogo [Complete collected works 
of Academician M. B. Ostrogradskit], v. 2: Lekisit Algebraicheskogo i 
Transfsendentnogo Analiza [Lectures on algebraic and transcendental 
analysis], Moscow-Leningrad, Academy of Sciences, 1940, 464 p. 17 X 25 
cm. Bound, 19 roubles. 


This volume of Ostrogradskii’s works contains (p. 433-462) his tables of indices and 
powers of a primitive root modulo , for all primes under 200. This set of tables first appeared 
in Akad. Nauk, S.S.S.R., Leningrad, Mémoires, . . . Sci. Math. Phys. et Nat.s.6,v.3 = Sci. 
Math. Phys., s. 6, v. 1, “livraison 4,” 1836, p. 359-385, and apparently was the first of its 
kind to be published. These tables were reproduced and extended by JAcosr in 1839 to 
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p < 1000 to form his famous Canon Arithmeticus (Compare MTAC, v. 1, p. 440). Twelve 
errata were discovered in Ostrogradskii’s tables by Jacobi after the latter had had them 
set in type: 


? Table Arg. For Read 
71(439) I 16 15 22 
71(439) I 26 22 15 
83(440) I 25 8 80 

127(447) N 105 107 108 
127(447) N 116 31 71 
137(449) N 108 88 87 
167(455) I 57 128 28 
173(456) I 57 72 92 
181(458) I 16 165 172 
181(458) I 26 172 165 
181(458) N 78 94 64 
193(460) N 155 173 174 


More than a century later the tables are now reproduced with the same old errata. 
It should be noted that Ostrogradskii’s tables give also all the primitive roots of each 
prime, information not presented in Jacobi’s Canon. 


D. H. L. 


EpiroriaL Note: Ostrogradskii’s portrait is on a plate opposite p. 64 of Les Mathéma- 
tiques dans les Publications de l’ Académie des Sciences 1728-1935. Répertoire Bibliographique, 
oscow, Academy of Sciences, 1936. In this v., p. 108, the date of publication of Ostro- 
gradskii’s ‘“‘Tables des racines . . .” is given incorrectly as 1838. The blue cover of “‘livrai- 
son 4” in the Harvard University library copy is dated 1836. 


368[F].—WILHELM Patz, Tafel der regelmdssigen Kettenbriiche fiir die Quadrat- 


wurzeln aus den natirlichen Zahlen von 1-10000. Leipzig, Akademische 
Verlagsgesellschaft, 1941. Lithoprinted by Edwards Bros., Ann Arbor, 
Michigan, 1946, xvi, 282 p. 15 X 22.9 cm. $6.50. Published and dis- 
tributed in the public interest by authority of the Alien Property Cus- 
todian under license number A-412. 


The regular continued fraction representing the square root of a positive non-square 
integer has been the subject of much experimental work and theoretical investigation since 
the time of EuLer. This table will serve as a useful tool in the further work along these lines. 
For each positive non-square D < 10002 are given the periodic partial quotients i, bs, - -- 
b, in the expansion 


1 1 1 1 


* 
vwD=b = a. aa 
aoe Ores ee ee 


Here bp = 2bo = 2[VD}. Since b, = bp_x(k > 0), it suffices to give by for k < 4p and by; 
and this is done except when p < 6. In a majority of cases p is even (p = 2g) and by is 
printed with an asterisk. Thus for D = 178 and 209 the entries are 

178|13(2, 1, 12*, 1, 2, 26), 209|14(2, 5, 3, 2*, ---, 28). 
In case. p is odd a diamond is printed before the expansion. Those values of D which are 
primes are followed by a small p in the argument column. 


The usual method of expanding VD is explained on p. xi-xiii. Writing the mth complete 
quotient in the form 


Xn = (VD + Px)/Qn = bn + Xm 
we have the four formulae of recurrence 
bn = (bo + Pn)Qn™, Pasi = bnQn — Pr, Qn41 = bn(Pni — Pn) + Qn, 
Quai = (D — Pii41)/Qn 
The last of these was used as an “automatic check”’ on the exactness of the calculation. 
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It is perhaps worth noting that the first two formulae may be replaced to advantage by the 
following pair whenever the value of 5, is not at once obvious 

bo + Pa = bnQn + fr, Psi = bo — fn, 
where 0 <r ¢ Q. 

The usefulness (and also the number of pages) of this volume would have been more 
than doubled had the author included the denominators Q,, as given in the tables of DEGEN,! 
CayLey? and WHITFORD.’ These numbers are important in the application of continued 
fractions, especially to the diophantine equation of Lagrange 


x? — Dy = N 


Besides this application the importance of this table lies in the wealth of statistical 
information it gives about the expansion of square roots of integers. The table throws some 
light on the unsolved questions of whether the period p is even or odd, whether, for a given 
D, the central partial quotient has the value do or 69 — 1 or not, whether p < cvD; and 
so on. Inspection of the latter part of the table reveals quite a large number of very long 
expansions. There are 31 values of D for which the period exceeds 2VD. These range from 
1726 to 9949 with periods of 88 and 217 respectively. The ratio p/ VD reaches a maximum 
of 2.2245 at D = 7606. Thus there is still room for the conjecture that, for all D, p < V5D. 
These long expansions appear to have more than their share of unit values among their b’s. 
In fact more than 43.62 percent of their partial quotients are equal to unity. The average 
for all real numbers is only logs (4/3) = .41503. 

There are three errata listed on p. xiv: D = 2872, for 1, 2, 2, 4, read 1, 2,4; D = 4170, 
for 2, 1, 3, 3, read 2, 1,4, 3; D = 4966, for 1, 4, 1, 2, read 1, 4, 2, 2. The first still occurs in 
the 1946 edition while the last two have been corrected. Nevertheless the above list is given. 
This confusing bit of editing led the reviewer to recalculate the expansions for D = 4170 
and 4966. BEEGER has pointed out (/TAC, v. 2, p. 88) that'in the 1941 edition the diamond 
sign is printed one line too low at D = 6938, 6949, 6953, and 9698. The first three of these 
misprints occur also in the present edition but the diamond is two lines too low at 9697 and 
is missing at 9698. The author has compared his table with those of DEGEN (D < 1000), 
CayLey (1001 < D < 1500), Wuitrorp (1501 < D < 2012) and THreLMANN® (about 140 
isolated D’s under 10*). No errata in these tables are quoted although the first and third 
are known to contain 2 and 4 erroneous continued fraction expansions respectively. The 
table of Roserts® for all primes D = 4n + 1 < 10‘ was not available to the author. A 
comparison of these two tables would give a very good idea of the reliability of the one 
under review. 


D. H. L. 


1C. F. DEGEN, Canon Pellianus . . ., Copenhagen, 1817. 

2A. CayLey, “Report of a committee appointed for the purpose of carrying on the 
tables connected with the Pellian equation from the point where the work was left by Degen 
in 1817,’”” BAAS, Report, 1893, p. 73-120; also Collected Mathematical Papers, v. 13, 1897, 
p. 430-467. [These tables were computed by C. E. Bickmore. ] 

3E. E. WuitForp, The Pell Equation, New York, 1912, p. 164-190. 

4M. von THIELMANN, “Zur Pellschen Gleichung,”’ Math. Annalen, v. 95, 1926, p. 
635-640. 

5D. H. LEHMER, Guide to Tables in the Theory of Numbers, 1941, p. 138, 171. 

°C. A. Roserts, ‘‘Table of the square roots of the prime numbers of the form 4m + 1 
—— 10000 expanded as periodic continued fractions,” Math. Magazine, v. 2, p. 105-120, 
1892. 


369[F].—HEInrIcH TIETZE, ‘‘Einige Tabellen zur Verteilung der Primzahlen 
auf Untergruppen der Gruppe der teilerfremden Restklassen nach 
gegebenem Modul,” Akad. d. Wiss., Munich, Abh., Math. Nat. Abdt., 
n.s., no. 55, 1944. 31 p. 22.3 X 28.5 cm. 


This paper contains 26 short tables giving information about the distribution of primes 
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in certain sets of arithmetical progressions the last term of which is denoted by L. The actual 
forms considered are km + 1;(j = 0, 1, ---) for the 26 following values of m and 1;: 


Table m o(m) "% L 
1 10 + 1,9 571 
2 10 + 1 571 
3 8 4 1 449 
4 8 4 1,3 457 
5 8 4 1,5 449 
6 8 + 1,7 449 
7 9 6 1,4,7 487 
8 9 6 1,8 487 
9 30 8 1, 29 3511 
10 30 8 1, 11 3511 
11 30 8 1, 17, 19, 23 1913 
12 30 8 1, 7, 13, 19 1879 
13 30 8 1, 11, 19, 29 1901 
14 26 12 1, 3, 9, 17, 23, 25 1091 
15 26 12 1,3,9 1069 
16 26 . 12 1, 25 1091 
17 26 12 1 1093 
18 262 130 Quadratic residues 3931 
19 262 130 5th power residues 3929 
20 262 130 10th power residues 3911 
21 262 130 13th power residues 3929 
22 262 130 26th power residues 3467 
23 262 130 65th power residues 3929 
24 262 130 1 298943 
25 262 | 130 259 298153 
26 262 130 17 297911 


Under multiplication modulo m, the set of r’s in each case forms a group I, in fact a subgroup 
of the group H of the g(m) numbers < m and prime to m. 

Let I(x) and p(x) denote the number of primes < x belonging respectively to H and 
I modulo m. If there are h elements of I and if ¢(m) = hi, then, according to the prime 
number theorem (generalized), g(x) and ilIp(x) are asymptotically equal and approach 
¢$(m)x/m In x. The tables give values of these step functions together with the difference 


A(x) = Ia(x) — ilp(x) 


for x < L. 

Since A(x) is a step function it suffices to tabulate it only at the values of x where it 
changes value, that is at primes belonging to (mod m). These primes are denoted by NV 
and form the arguments of the tables. This makes il[p(N) merely a list of consecutive 
multiples of 7, and this column might well have been omitted. As x varies from one value 
of N to the next, A(x) increases because Ilg(x) increases. The value of A(x) just before the 
next value of N is denoted by A*(JV)) and is tabulated also. 

The modulus 262 is chosen because 131 is the least prime having 3, 5, 7, 11, 13 as 
quadratic residues, 17 and 259 are the least and greatest primitive roots of 131. The last 
three large tables, especially table 24 might some day prove useful as a list of primes of 
these forms. Cunningham's observation that the form km + 1 contains fewer primes < x 
than km + 1, 1 + 1, (1, m coprime) does not seem to hold for m = 262. In fact A(J) in 
table 24 changes sign very often. 


D. H. L. 


370[F].—I. M. VinoGrapov, Osnovy Teorii Chisel [Fundamentals of the 
Theory of Numbers], Moscow-Leningrad, (a) third ed., 10 000 copies, 
1940, 111 p. + an errata sheet. 12.7 X 18.7 cm. Bound, 3 roubles. 
(b) Fourth ed., 3 000 copies, 1944, 142 p. + an errata sheet. 13.8 20 
cm. Paper bound, 4 roubles. 
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(a) This interesting little volume contains two kinds of tables: 


(1) Tables of indices and powers of a primitive root modulo p for » < 100 (p. 104-109). 
These are based on least primitive roots and so are identical with tables of WERTHEIM,! 
and Uspensky & HEASLET.? A comparison with the latter table reveals no discrepancy. 

(2) Table of least primitive roots of primes p < 3000 (p. 110-111). 

Three errata may be noted: p = 1013, for 2, read 3; p = 2593, for 10, read 7; p = 2999, 
for 7, read 17. 

(b) In this edition the first group of tables (p. 135-140) is the same as in the third 
edition but the table of least primitive roots of primes < 3000, in the third edition, has been 
corrected and enlarged to primes < 4000 (p. 141-142). 

D. H. L. 


1G. WERTHEIM, Aufgangsgriinde der Zahlenlehre. Brunswick, 1902, p. 412-417. 
2 J. V. Uspensxy & M. A. HEASLET, Elementary Number Theory. New York and London, 
1939, p. 477-480. ; 


371[G, L].—A. CoLomBanl, “La théorie’ des filtres électriques et les poly- 

' nomes de Tchebichef,”’ Jn. de Physique et de Radium, s. 8, v. 7, Aug. 1946, 

p. 231-243. 21.3 X 26.6 cm. Compare MTAC, v. 1, p. 125, 149f, 385, 
RMT 381, 383. 


; There are two tables for so-called Chebyshev polynomials, p. 236-237. T. I gives S,(x) 

for x = — 2(.1)0, m = [2(1)10; Dj]. Also zeros of S, to 5D. T. II gives values of Xn(x) 
= S,(x) — Sn_i(x), for x1/x2 = x — 2 = — 4(.1)0, x = — 2(.1)2, m = [1(1)10; nD]. Also 
zeros of X1o(x) to 5D. Figs. 2-4, p. 234-235, are graphs of X,(x) for m = 1(1)10. 


Sa(x) = sin (n + 1)0/sin 0, x = 2 cos @. 


372(I].—H. E. Sauzer, ‘Coefficients for facilitating the use of the Gaussian 
quadrature formula,” Jn. Math. Physics, v. 25, 1946, p. 244-246. 
17.5 X 25.5 cm. 


In the Gaussian quadrature formula 


Si feedde =¥ aisles) + Re 
i=1 


the sum extends over the roots x; of the Legendre polynomial P,(x). As these roots are not 
equally spaced, it is not possible to test the smoothness of a set of computed ordinates f(x;) 
by straight-forward differencing, as one would do in the case of a Cotes type formula. 

To examine the (m — 1)st difference of this set of ordinates one must resort to divided 
differences. Much of the cumbersome calculation attending the general divided difference 
process can, in this case, be avoided. In fact this difference can be written 


PH Ci f (xi) 
i-1 
where the coefficients are simply 
Ci) = 2-"(2")/Pa’ (xi). 


This paper contains a small table of these coefficients for n = 3(1)10. The accuracy is 8D 
for n = 3,4, 5; 7D for n = 6, 7, 8; and 6D for n = 9, 10. The author fails to indicate that 
in his notation the roots x; are so ordered that 


%1 <X%2< ees <x, 


a fact which the user of the table will need to know. 
D. H. L. 
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373{[I, L].—H. A. RADEMACHER & I. J. SCHOENBERG, “An iteration method 
for calculation with Laurent series,” Quart. Appl. Math., v. 4, July 1946, 
p. 142-159. 17.5 X 25.4 cm. 


In the authors’ words, “the purpose of this paper is to describe a method whereby ra- 
tional or algebraic operations with Laurent series may be performed with high accuracy at 
the expense of a reasonable amount of labor.” 

The main problem considered is that of solving numerically 


(1) S(w, 2) = ao(z)w™ + ai(z)w™"! + --- + an(z) = 0 


for the coefficients of the Laurent series of a particular branch of w(z), where the a;(z) are 
regular and uniform functions of z in the ring 


R: 1 <|z|< re 


and where neither ao(z) nor the discriminant D(z) of (1) is zero in R. 

The procedure suggested by the authors is first to find an initial approximation and 
then to use an iteration scheme based on a modification of Newton’s algorithm. In this 
modification only one division is needed and that is a preliminary one. A method is given 
of obtaining by trigonometric interpolation a Laurent polynomial, 


F(z) = D0 tas, 
-n 
as a first approximation. It is proven that F,(z) approaches the solution as nm > ©. One of 
the points of this paper is that it is preferable to start with a small value of m and then to 
iterate rather than to use F,(z) for a large value of n. 

The iteration scheme is as follows: Since the discriminant D(z) is not zero there are 
polynomials ¢(w) and ¥(w), with coefficients which are polynomials in the a;(z) divided 
by D(z), such that, 
df(w,z) ~ 1 

dw 


The modified Newton's algorithm can now be expressed by the recurrence formula 


o(w)f(w, z) + yw) 


Wr41 = We — f(we, 2)p (wr). 
This has the usual quadratic convergence of the Newton algorithm. The authors state that 
this method has been used previously by Schwertfeger for the numerical solution of ordinary 
algebraic and transcendental equations. 
The authors show that in the special case of solving 


a(z)w(z) —-1=0 
the recurrence formula reduces to 
Wri = We(2 — attr), 


which is the formula described by Hotelling for inverting matrices. Therefore, in reciproca- 
tion of a Laurent series one can use an inequality of Hotelling and Lonseth to obtain a 
limit for the error due to stopping after any number of steps. 

As an illustration the authors compute the coefficients w, of the Laurent expansion of 
the reciprocal of — Tov 13z) between the first two positive roots of that function. The 
entire computation to 9D is exhibited in tabular form for wa with — 29 < m < 32. The 
remaining coefficients are numerically smaller than 10-*. The paper also contains a descrip- 
tion of how the methods of calculating with Laurent series apply to calculations with 
absolutely convergent Fourier series. 

In conclusion one can say that this article presents in a very convenient form a solution 
to the problem considered, especially for those who will have to do actual computations 
of this sort. 

ABRAHAM HILLMAN 
NBSMTP 











258 RECENT MATHEMATICAL TABLES 


374[K].—La Mont C. Cote, “A simple test of the hypothesis that al- 
ternative events are equally probable,” Ecology, v. 26, 1945, p. 204. 
16.5 X 25.4 cm. 
! 


ke n! 
f P = 2») ———__, for k = 0(1)12, m = [2(1)35; q 
Table III, values o L ie —p! or 0(1)12, = [2(1)35; 5D] 


“Table gives proportion in both tails of (3 + })*. For larger values of m uset = (n — 2E)n74. 
In general, a value is statistically significant (P < 0.05) if E < 3n — ni.” 
Extracts from text 


375(K].— FREDERICK E. Croxton & DUDLEY J. COWDEN, ‘‘Tables to facili- 
tate computation of sampling limits of s, and fiducial limits of sigma,” 
Industrial Quality Control, v. 3, July 1946, p. 18-21. 21.6 X 27.9 cm. 


For samples of size N drawn from a normal distribution with known variance o”, upper 
and lower percentage points of the distribution of s/o are given in Table 1 entitled: “Values 
of s/a at selected probability points for various sample sizes.’’ The sample standard deviation 
is s = [>> (x — £)*/N}}. The probability points of the distribution of s/o are given in pairs 
for probabilities a and 1 — a, witha = .001, .005, .01, .025, .05, .10, for N = [2(1)30; 3D]. 
An approximation is given for N > 30. A table similar to the present one for a = .001, .005 
and for N = 2(1)15 is given in Amer. Standards Assoc., Control Chart Method of Controlling 
Quality Duriug Production, no. ASA Z 1.3—1942, p. 40. 

Table 2, entitled ‘‘Values of o/s for use in computation of selected fiducial limits of « 
for various sample sizes,” may be used to obtain confidence or fiducial limits for ¢. Con- 
fidence levels available are .998, .99, .98, .95, .90, .80 for N = 2(1)30. The entries in this 
table are principally reciprocals of the entries of Table 1. A similar table for confidence 
levels .90, and .98, N = 5(1)30, appears in E. S. Pearson, The Application of Statistical 
Methods to Industrial Standardisation and Quality Control, London, British Standards 
Institution, 1935, p. 69. 

The tabulated values of the tables under review were derived principally from 
CATHERINE M. THompson, “Table of percentage points of the x? distribution,’ Biometrika, 
v. 32, p. 187f, 1941. (See MTAC, v. 1, p. 78). But the .999 points of Table 1 were derived 
from R. A. FisHER & F. YATEs, Statistical Tables for Biological, Agricultural, and Medical 
Research. London, 1938, Table IV, p. 27, while the .001 points of Table 1 were derived from 
tables of F shown in F. E. Croxton & D. J. CowvEn, A pplied General Statistics. New York, 
1939, p. 878-879. 

FREDERICK MOSTELLER 
Harvard University 


EprtoriAL Note: It may be remarked that it was exactly on the pages quoted as 
sources, FISHER & YATES, 1938, p. 27, and Croxton & CowDEN, 1939, p. 878, that we have 
listed errors in the tables in question, namely: MTAC, v. 1, p. 324, and 86. 


376[K].—V. L. Goncuarov, Teoriia Veroiatnostet [Theory of Probabilities}. 
Moscow and Leningrad, 1939, 427 p. + errata slip. 14.4 & 21.7 cm. 
Bound, 11 roubles. An edition of 5000 copies. 


This government ordnance industry publication has three small tables on its last seven 
pages. In the notation of the FMR, Index, these are 


(a) H(x) = 2a e~ dt, for x = [0(.05)2.2; 4D], [2.2(.05)2.75, 3; 6D], [3.5, 4; 9D]. 
(b) H(px), where p = .4769362762:--- is the root of H(x) = 3, for x = [0(.01)3.4(.1)5.4; 5D]. 
(c) xH(px) + p-'x-te-?*2" for x = [0(.05)5.2; 4D]. 


All three tables give first differences. 
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Strange to say, the very well known function H(x) is not tabulated correctly in (a). 
There are two errata: t = 2.4, for .999312, read .999311; t = 2.7, for .999868, read .999866. 
A last-figure error occurs in (b); in the final entry ¢ = 5.40, for 99972, read 99973. 

The function tabulated in (c) is essentially the second iterated integral: 


1 2 fs fe 
mi 8 e*doedt. 
oe} 
This table appears to have several last-figure errors. For example, t = 1.05, for 1466, 
read 1468; t = 1.10, for 1493, read 1494. 


For a discussion of the tables the reader may consult p. 150 and 217. 
DD: Ss 


377(K].—Mrs. CATHERINE M. (THompson) GryLtis & Mrs. MAXINE 
MERRINGTON, “Tables for testing the homogeneity of a set of estimated 
variances,’ Biometrika, v. 33, June 1946, p. 302-304. A preface by H. O. 
Hart.ey & E. S. PEARSON occupies p. 296-301. 19.3 X 27.3 cm. 


These tables are designed to provide 1 per cent. and 5 per cent. points for testing the 
hypothesis that the variances of several normal populations, as estimated from two or more 
independent observations on each, are equal; or hypotheses equivalent thereto. The basis 
for computation was an approximation of HARTLEY! to the distribution of a statistic sug- 
gested by BARTLETT? which differs in weighting factors from the likelihood ratio statistic 
proposed by NEyMAN & Pearson,’ and which has been found more powerful in certain 
cases by BisHop & Narr.‘ 

Hartley’s approximation depends on three parameters: the number of populations &, 
and two functions of the —_— of freedom » of the estimates of the population variances, 


1 1 
a=Z(= )=5, »a-L—- where N = Ln 


7 Me N3’ 
The percentage points vary but ir with c3. The tables are double ony, giving for each 
pair k and c, two values which are approximately the extremes with respect to variation 
in c;. An auxiliary table aids interpolation with cs in the occasional case when these extremes 
cover the computed statistic. 

The percentage points are given to 2D(3 or 4S) for k = 3(1)15 and (the entire range of) 
¢, = 0(.5)5(1)10(2)14. A historical note, several illustrative examples and a discussion of 
the accuracy of the approximation, are provided in the preface. It is found that the approxi- 
mation is ‘‘very good” if the degrees of freedom all exceed 2, and is “‘adequate”’ if some are 
as small as 2. 

J. L. Hopces, Jr. 

Statistical Laboratory, 
Univ. of California, Berkeley. 

1H. O. Hart.ey, “Testing the homogeneity of a set of variances,” Biometrika, v. 31, 
1940, p. 249-255. 

2M. S. BartLett, “‘Properties of sufficiency and statistical tests,” R. Soc. London, 
Proc., v. 160A, 1937, p. 268-282. 

J. NeyMaNn & E. S. PEARSON, “On the problem of & samples,” Akad. umiejetnoéci, 
Bull. Intern., 1931A, p. 460-481. 


‘D. T. Bisnor & U.S. Narr, “A note on certain methods of ating for the homogeneity 
of a set of estimated variance,’ ’R. Statist. Soc., Jn., v. 102, Suppl. v. 6, 1939, p. 89-99. 


378[K, V].—Crecit Hastincs & MARGARET PIEDEM, Miscellaneous Prob- 
ability Tables, calculated and checked under the direction of Dr. H. H. 
GERMOND, 1942-1944. Applied Mathematics Panel, National Defense 
Research Committee, Note no. 14, New York, July 1944. ii, 65 p. 
Offset print. 21.3 X 27.8 cm. These tables are not available for public 
distribution. 
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The main table (p. 6-37 and introduction p. 1-5) is devoted to the two-dimensional 
normal distribution function (or error function) 


V(b, q) = (2x) ™ [" edydx, W = x8 + 9. 


The double entry table to 5D has h and g/h as independent variables, the former ap- 
pearing in rows, the latter in columns; h = 0(.01)4, g/h = .1(.1)1. Differences are tabulated 
in the main rows and columns between the corresponding entries. The differences within 
columns have usually only one digit and never exceed 85 units; the differences in the rows 





keep to three digits. The need for the tables arose in probability problems associated with | 


bombing and fragmentation damage. NBSMTP supplied the key values of V(h, g) from 
which the Table was subtabulated. The computation was carried on for some time before 
the appearance of C. Nicuotson, “The probability integral for two variables,” Biometrika, 
v. 33, part 1, April, 1943, p. 59-72, where a table of V(h, q), to 6D, is given for h = .1(.1)3, 
q = .1(.1)3, «. ; 

Pages 39-45 cover tables for H(x) — xH’(x), where H(x) is defined by 


H(x) = 2(x)-4 if * e-Hat. 


The range is [0(.001)3(.01)4.3; 7D]; for x > 4.3 the entries would equal unity throughout. 


(The entries in the first column on p. 45 are misprinted: the last zero should be deleted | 


everywhere, and the last line on p. 44 should be deleted. This table was computed by using 
NBSMTP, Tables of Probability Functions, v. 1, 1941.) It is stated that the rounding errors 
will occasionally amount to one unit in the last digit. There follow, p. 46-49, tables of the 
inverse of the function H(x) — xH’(x) covering the entire significant range, namely 
(0(.001)1; SD]. 

The next table, p. 51-62, gives values of the function y = 1 — (1 + x)e~*. These tables 
will be useful in particular in connection with the Poisson distribution. The range is 
[0(.001)5(.01)10(.1)15; 5D]. Again, y is practically constant for x > 15. 

Two small tables conclude the collection. On p. 63 we find values of the product xy 
where y is defined by the equation 


xy=1-—e 


The range is x = [0(.01)1; 5D], A. Finally, on p. 65 is a table of 


= xe-2? [| * ot? 
o(x) = xe J e“dt 
for the range [2(.1)7; 5D], A*. 
W. FELLER 
Cornell University 
EprroriaL Note: In our notes on Dawson’s or Poisson’s integral we have listed an 
earlier table of ¢(x), MTAC, v. 1, p. 323, N. Kapzov & S. GwospoweEnr, Z. f. Physik, v. 45, 


1927, p. 133. This table is for the range x = [.1, .5, .8(.2)1.2(.05)2.2; 5D]; A, 1.5-2.2. The 
use for such a table there, arose in discussion of oscillations in electron tubes. 


379[L].—D. CHaLonce & V. KourGaAnorr, “Recherches sur le spectre 
continu du soleil,’ Annales d'Astrophysique, v. 9, 1946, p. 69-96. 
21.5 X 27.4 cm. 


Appendix I contains two tables of “‘la fonction I incompléte d’argument négatif”’ : 


T.(«) = J * yo-tetidt 
T. A, p. 94. x = (0(.01).1;4D]; @ = [.1(.1)1], 
T. B, p. 95-96. a = [.01(.01)1; 4D]; x = [0(.1)1.1]. 
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The authors state that the tables constitute ‘‘un extrait, pour le domaine qui nous 
intéresse ici, d'une table plus étendue qui paraitra prochainement.” No details of the 
calculation are given. 

L. E. CUNNINGHAM 
Astronomy Department 
University of California 
Berkeley 


380[L].—HarvARD UNIVERSITY, Computation Laboratory, Annals, v. 3: 
Tables of the Bessel Functions of the First Kind of Orders Zero and One; 
v. 4: Tables of the Bessel Functions of the First Kind of Orders Two and 
Three, by the Staff of the Laboratory, Professor H. H. Aiken, Technical 
Director, Cambridge, Mass., Harvard Univ. Press, 1947. ii, xxvii, 652 p. 
and viii, 652 p. 19.5 XK 26.7 cm. $10.00 + $10.00. Compare MTAC, v. 2, 
p. 176f, 185f. The offset printing of these volumes is of outstanding 
excellence. 


P. iii, ‘Staff of the Computation Laboratory” [11 members and 13 assistants listed]. 

P. vi, ‘Preface’ by Professor Aiken. Quotations: In Nov. 1944 a conference at the Naval 
Research Laboratory was called to discuss the tabulation of Bessel functions of the first kind 
and of high order, which at that time were needed in connection with various research 
problems of interest to the Navy. As a result the computation project of the Bureau of Ships 
was directed to tabulate the required functions. During discussion it became clear that if 
J100(x) was to be accurate to ten decimal places, an adding, multiplying, and storage capacity 
of more than forty digits would be required of the Automatic Sequence Controlled Calcu- 
lator. Since the multiplying unit of the calculator already supplied forty-six digits and the 
algebraic sign, it was only necessary to link two normal storage registers, each comprising 
twenty-three digits and the algebraic sign, to form a single adding storage register covering 
forty-six digits and the algebraic sign. 

On the first of January 1946, the Computation Project was transferred to the Bureau 
of Ordnance. 

P. ix-xxii, “Introduction” by RicHarp M. BLocu 


Part I. The Bessel Functions, p. ix—xi, 
Part II. The Computation of the Tables, p. xii—xviii, 
Part III. Interpolation in the Tables, p. xix—xxii. 


Part I. Because of the important applications of Bessel Functions to many physical phe- 
nomena, they have been the subject of intensive investigation for many years. Recent 
advances in the theory of frequency modulations, resonance in cavities, waves in various 
media, vibration theory of structures, and other problems of physics and engineering have 
greatly increased the need for extensive tables of J,(x) and Y,(x) covering a large range 
both of the order and of the argument. 

The Staff of the Computation Laboratory is at present engaged in the tabulation of 
Jn(x), 0 < x < 100, m = 0(1)100. The present volumes contain tables for m = 0(1)3, 
x = [0(.001)25(.01)99.99; 18D). 

The Bessel functions satisfy the two relations 


In—1(x) — Jn+i(%) = 2J,'(x) 
2 
(1) Jn-%) + Jnga(x) = = Jn(x). 


For machine computation, the successive application of the recurrence formula (1) provides 
the most feasible method of obtaining the high order functions. Since ten decimals place 
accuracy is to be maintained in the tables of J,(x) for 4 < m < 100, it was necessary to 
investigate the cumulative loss of accuracy which arises in the repetitive use of (1) as the 
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computation proceeds. If Jm_i(x) and Jm(x) (m = 3 for x < 2, m = 2 for x 2 2) are the 
two basic functions upon which the recurrence is constructed, the maximum number of 
decimal places lost is eleven. Consequently the low order functions Jo(x), Ji(x) and J2(x) 
were computed correct to twenty-three places of decimals. Since the figures were available, 
the tables of J,(x) (n = 0, 1, 2, 3) have been printed to 18D, despite the fact that interpola- 
tion within these tables to full accuracy would be extremely difficult with the present 
manual aids to computation. 

The Automatic Sequence Controlled Caiculator is so arranged that after all final results 
are automatically checked, they are printed by the typewriters controlled by the machine 
itself. Certain tables of Bessel functions to eighteen or more D were read against the values 
computed at the Computation Laboratory. These comparisons were made with the values 
of the functions listed in the following three tables: MEISSEL, as in Gray, MATHEWws, & 
MacRosert, A Treatise on Bessel Functions, 1931, p. 286-299, 18D, 1 < x < 24, Ax = 1, 
n = 0(1)3; Hayasui, Tafeln der Besseischen, Theta-, Kugel-, und anderer Funktionen, 1930, 
p. 52-59, m = 0(1)3, .01 < x < 100, selected Ax, 22-103D; Apis, R. Soc. London, Proc., 
v. 66, 1900, p. 40-43, n = 0, 1;0 < x < 6, Ax = .1, 21D. No discrepancies were observed. 
Part II. Values of J,(x) given in other tables were not used. All numerical constants in- 
cluding the coefficients of the ascending power series, the asymptotic series and those re- 
quired for interpolation, were evaluated at the Computation Laboratory, regardless of the 
availability of such material from external sources. 

The ascending power series 


> r x 2r+n 
= 2r+n, we 
J n(x) “Lk ann( =) ’ 
where d;,, = (— 1)"/r!(n + r)!, and k is a normalizing factor, was used to evaluate J,(x) 
for m = 0(1)3 over the range 0 < x <2 with increment Ax = .001, and for m = 0(1)2 over 
the range 2 < x € 25 with increment .01. The values of 10”-a,,, computed to 50D, are 
given for r = 1(1)60,0 < m < 168, m = 0(1)2. 
For the range 0 < x <2,k =1;2£ x < 10,k =5;10 < x < 20,k = 10; 20 < x < 235, 
k = 20. For m = 3 we have |10"-a,,3|, r = 0(1)15, 1 < m & 28, p. xxiii—xxxi. 
There are similar tables (p. xxxii-xxxvii) for the various asymptotic expansions. 
Extracts from introductory text 


381[L].—C. W. Jones, J. C. P. MILter, J. F. C. Conn, & R. C. PANKHURST, 
“Tables of Chebyshev polynomials,” R. Soc. Edinb., Proc., v. 62A, no. 21, 
1946, p. 187-203. 17.5 &K 25.5 cm. 


The main object of the article under review is to present a table of the Chebyshev 
polynomials C,(x) = 2 cos (m arc cos $x) for m = 1(1)12 and x = 0(.02)2. The tabulated 
values are either exact or given to 10D. In addition to the table of C,(x), the article contains 
also short tables of the functions (4 -- x*)4, (2 + x)#, (2 — x)! and arc cos $x required in 
the applications of Chebyshev polynomials discussed in Dr. Miller's article, “Two numerical 
applications of Chebyshev polynomials” (RMT 383). 

The tabular material is preceded by an excellent introduction giving the definition of 
the Chebyshev polynomials C,(x) and S,(x) and of other related functions, the differential 
equations and recurrence relations satisfied by these functions, the explicit power series 
expressions of these functions, the expressions of the twelve powers of x in terms of Cheby- 
shev polynomials C,(x), the orthogonality relations and the generating functions for each of 
the functions under consideration. 

The reviewer agrees’ with the authors’ remarks that the tables will be of particular 
importance to computers. One application particularly worth mentioning is the process of 
interpolation by means of Chebyshev polynomials (RMT 383). The efficacy of this process 
of interpolation is illustrated by the following observation: In a certain region of the Mathieu 
functions ms. in preparation by the NBSMTP, interpolation to the full accuracy of the 
table would require the use of a formula involving differences up to the ninth order; the 











corres} 
polyno 

Fc 
spondi 
NBSN 


(1) 


where 
perio 
perioc 
rise t 
curve 
tions 

bn +1) 


the fe 


In (2 
accef 
every 
infini 


Gm (q) 
Whe 


that 
of th 


(3) 


wher 
regic 
the : 


in w 
(a, ¢ 


upo 
seri 


(4) 


re the 
ber of 

J2(x) 
ilable, 
rpola- 
resent 


esults 
ichine 
values 
values 
ws, & 
= 1, 
1930, 
Proc., 
rved. 
ts in- 
se re- 
of the 


Jn (x) 
! over 
), are 


IN 
i) 
vu 


text 


IRST, 
». 21, 


yshev 
ilated 
itains 
ed in 
erical 


on of 
ential 
series 
heby- 
ich of 


cular 
ess of 
ocess 
thieu 
f the 
; the 








RECENT MATHEMATICAL TABLES 263 


corresponding jnterpolation formula in terms of C,(x) requires only the first five Chebyshev 
polynomials. 

Forty percent of the entries of the table under review were proofread against the corre- 
sponding entries in the more extensive table of Chebyshev polynomials prepared by the 
NBSMTP (see MTAC, v. 1, p. 125); no discrepancies were discovered. 


ARNOLD N. LowAn 


382[L].—N. W. McLacu.an (a) “Computation of the solution of Mathieu’s 
equation,”’ Phil. Mag.,s. 7, v. 36, June 1945 (publ. Jan. 1946), p. 403-414. 
17 X 25.5 cm. (b) “Mathieu functions and their classification,’ Jn. 
Math. Phys., v. 25, Oct., 1946, p. 209-240. 17.3 K 25.4 cm. 


(a) This paper deals with the computation of the solutions of 
@ 
(1) <2 + (a — 2g cos 22)y = 0 
dz? 


where a, q are real parameters. For certain characteristic values, the solutions of (1) are 
periodic, of period + or 2%. Those characteristic values which give rise to even solutions of 
period x and 2x are denoted by dam and d2m41, respectively. The characteristic values giving 
rise to odd solutions of period 2x and x are denoted by bem4: and bem42, respectively. The 
curves @ = a»(q) and a = bm(g) separate the (a — q) plane! into regions in which the solu- 
tions are ‘‘stable’”’ or “‘unstable.” When the parametric point (a, g) lies between a, and 
bmi, the solution is ‘‘stable’’; that is, two independent solutions of (1) may be written in 
the form 


= cos 
(2) y= Lia r+pt ss. 


r= 


In (2), p = O if the subscript m in a, is even, and p = 1 if m is odd. According to currently 
accepted theory, there exists a unique positive value of 8 less than unity corresponding to 
every point (a, q) in this region, such that the solution (2) remains finite as z approaches 
infinity through real values. 

The value of a which, for a fixed g, determines 8, will be denoted by an4g. Between 
Gm(q) and bm+:(g) there lies a family of iso-8 curves, i.e., 8 = constant. When 8 turns out to 
be a proper fraction p/s (in its lowest terms), the solution will be periodic, of period 27s. 
When 8 is irrational, the solutions of (2) will be non-periodic. 

If the parametric point (a, q) lies between 6» and dm, the solutions are “‘unstable”’; 
that is, no solution of the form (2) exists. By Floquet’s theorem, there does exist a solution 
of the form 


) 


(3) y =e Do celtrtpizi, 

r=—@2 
where p = 0 if m is even in bm, dm, and p = 1 if m is odd. When (a, g) lies in an unstable 
region, u is real. It may be readily seen that when yu is a purely imaginary number, (3) yields 
the solutions (2). 

The most important contribution of the paper is to show how 8 may be determined, 
in the stable region, from a knowledge of the characteristic values a», and b,. Let a point 
(a, g) of a stable region be given and let it be desired to determine 8. The author improves 
upon Ince’s method by obtaining some good first approximation to B. By inverting the known 


series for “‘a’’ in terms of (m + 8) and gq, the author obtains the following approximation! to 8: 


= a: a 
2a@-—-1)?—q@ 32(a@ — 1)%(a — 4) 





(4) p= [a 


In (b) the above expansion is extended to include, in the radical, the term — (9a? + 58a 
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+ 29)g*/64(a — 1)5(a — 4)(a — 9). This reviewer believes that, if the first three terms are 
left in their present form, the term involving g® should be? (— 37a? + 319a? — 587a + 17)q8/ 
64(a — 1)5(a — 4)?(a@ — 9). However, the series (4) is useful only when the expression under 
the radical converges rapidly enough with the given terms; and whenever the contribution 
from the corrected term involving q° is small, the uncorrected term will not be much larger 
numerically. 

When (4) cannot be used (a denominator may vanish or g/(a — 1) may be too large), 
then other approximations may be used. Let am and bm4: be the characteristic values for 
the given g, between which a lies. Let \ = (a — am)/(bms1 — Gm). Then X should be an 
approximation to 8. However, usually \ is too crude; a better approximation may be ob- 
tained if it is assumed that, at g = 0, the iso-8 curve intersects the region between a, and 
bm41 in approximately the same ratio, \. Since, at g = 0, am = m?, and bm4: = (m + 1)?, 
it follows from (4) and the above assumption as to \ that 


° pam {isn (EYP 1p 


The author gives still another empirical formula for g. Let 





gem = (at — ant)/(bm414 — On'). 
Replacing \ by ¢m in (5), one obtains 


: pam {fit (MEP -1} 


In an appendix, the author gives a table showing the accuracy of the several approximations 
in a few of the instances in which they were tried. The schedule given below summarizes 
the author’s results; all figures except those in the last column were recalculated by this 
reviewer; in cases of discrepancy, the recalculated figures are given. The results obtained 
were close to those of the author, except on line 2 of the schedule, where the author obtained, 
by formula (4),an amount .57—apparently by neglecting the third term under the radical. 





Formulae Used 
More Accurate 

m a q x (4) (S) (6) Value of 8 

1 3 2 48 _ .56 59 .579 

2 8 4 48 52 53 55 59 

2 6 2 21 3,34 Pe i 27 34 

5 36 16 48 3,62 50 Be | 3,583 

1 2 .0125 33 4142 Al 49 4142 

0 —1.45 2 .52 — 72 .72 52 


The last two examples were supplied by the reviewer. Except in the last example, 
formulae (4), (5), and (6) are better than \, and (4) gives a good approximation in many 
cases. The author recommends the use of this formula wherever possible. It is to be noted 
that the author’s ingenious method of improving on d in (5) and (6) is fruitful in its results, 
especially when (4) cannot be used. 

Once an approximation to 8 has been obtained, the method of improving it by iteration 
or interpolation, in the process of computing the coefficients ¢,, is fairly easy. Thus in the 
recurrence relation 


(7) Ca — (27 + p+ B)* ler — g(cr41 + Crs) = 0, 
one may neglect ¢.(r41) for r sufficiently large, and compute in turn ¢cy_1, «--, Co in terms of 
cr (hence also c,, ---, ¢, in terms of co); and again c_,41, ---, co in terms of c_, (hence also 


C_r41, etc. in terms of co). Then setting r = 0 in (7), the relation between c_1, ci, and Co 
will be satisfied only if 8 is correct, and the divergence of the right-hand side of (7) from zero 
shows how to correct 6 and the coefficients c,. One may of course compute the coefficients in 
terms of any c» rather than in terms of co. Several variations of the computing technique 
are given, with methods of checking the computations. 

The method of approximation may also be used for the unstable regions, once iso-p 
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Soe have been plotted. It is recommended that the coefficients c, be normalized so that 
cf = 1. 

(b) Here are given a great variety of representations for the solutions of Mathieu's 
equation and of Mathieu's “‘modified” equation. Quoting from the author: “The number of 
representations in the guise of series, integral relations, etc., exceeds 300. Of these, about 
200 have not been published hitherto. . . . No attempt is made to show the derivation of 
the new formulae, as this paper would then be much too long.” 

The first part of the paper deals with solutions of the first and second kind for Mathieu 
functions of integral order, and includes the very useful Bessel-function-products solutions, 
previously given in a paper by W. G. Bicktey & N. W. McLacuarian, MTAC, v. 2, Jan. 
1946, p. 1-11. It may be worth pointing out that Bessel-function-products solutions of 
Mathieu's differential equation were given by BRUNo S1EGER‘; and although Sieger’s work 
is again mentioned by Strutt,’ the importance of such solutions seems to have been little 
understood until it was emphasized by Bickley & McLachlan in their paper. (This reviewer 
learned of Sieger’s work from Professor Bickley.) All forms given in the January paper are 
included in this larger one by McLachlan, now under review. In addition, a great many 
variations of the Bessel-function-products are given, which may prove useful from a 
computational standpoint. 

Functions of the third kind (analogous to the well-known Hankel functions) are defined 
in section 11. Except for the normalization factor, these solutions are the same as the ones 
previously defined by L. J. Cau and J. A. Srratron (Jn. Math. Phys., Aug., 1941), see 
MTAC, v. 1, p. 157. The relations between the various solutions, when the parameter g is 
either positive or negative, are also given. 

In addition to formulae relating to Mathieu functions of fractional order (stable solu- 
tions) the author devotes considerable space to the “‘unstable” solutions. It is shown that, 
when the solution is put into the form 


yi = ef” pa Corp pe rte: 
ro 

(with p = 0 if (a, q) lies between a2», and bem and p = 1 if (a, q) lies between ben; and 
G2m+i) then C2, and c_s are conjugate complex numbers, if expressed in terms of Co, real; 
furthermore if p = 1, then Cor, and (¢:/c_1)c_2r_1 are conjugate, if c, is taken to be real. 
It is shown that there exists a real solution of the differential equation which tends to 
zero as z>> — © through real values, if a, g, and u(> 0) are real. Solutions denoted by 
C€msy(+ 5, q) are defined, analogous to the solutions cém(+ z, q) for integral m; similarly 
for S@m4y(+ 3, g). Such solutions are neither even nor odd. It is shown how to construct even 
and odd solutions, but they are less useful than cém4,(+ 2, g) and Sém4y(+ 2, @). 

For the same g and uy, there are two values of a. One is such that the solution approaches 
Cém(z, q) as u4—> 0; the second a corresponds to the solution which approaches se»(z, q) as 
u— 0. Solutions of fractional order corresponding to Mathieu’s modified equation are 
also given. 

The paper contains a number of asymptotic expansions for Mathieu functions of in- 
tegral order, both for large z and large g. These expansions involve certain multipliers (re- 
sulting from the normalization adopted) which cannot readily be expressed asymptotically. 
To this extent the solutions for large g, in a practical case, are really never obtainable by 
the given asymptotic formulae—they are known except for those multipliers which are 
functions of q. 

It is this unfortunate property of the normalization adopted by McLachlan (and the 
English school generally) which is at the crux of the divergence of opinion, regarding the 
normalization scheme, between the English and American schools. 

The author concludes with a section on the zeros of the functions and another on the 
classification of the various solutions of integral order, fractional order, and the unstable 
solutions. A useful iso-u chart, the data for which are credited to Dr. L. J.§Comrre, is 
also given. 

The paper is concisely written and represents a prodigious effort, both as to span 
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covered and the variety of formulae given. It forms a priceless compendium of known results 
(with a very considerable portion of them due to the author himself). The summary given 
above by no means covers all topics treated. It is hoped that the book promised by the 
author, enlarging on the theory covered in this paper, may soon be forthcoming. 


GERTRUDE BLANCH 
NBSMTP 


1 In the author’s om, the horizontal lines in the a-q plane are parallel to the ‘‘g” 
axis. Hence it might have been better to refer to the ‘‘gq-a’”’ plane, and to the taal 
point as (g, a). We shall not, however, depart from the author’ s notation. 

2 Mrs. Ipa Ruopes of the NBSMTP checked this reviewer’s inversion. Dr. McLachlan 
pointed out that the second term of (4), given incorrectly in (a), is correct in (b). 


. 3 These three entries are corrections of the printed values, furnished by Dr. McLachlan. 
DITOR. 


4B. Srecer, “Die Beugung einer ebenen elektrischen Welle an einem Schirm von 
eT uerschnitt,” Ann. d. Phys., s. 4, v. 27, 1908, p. 626-664. 
Strutt, Lamésche- Mathieusche und, verwandte Funktionen der Physik u. 
Techni, diirgebnasse a. Math., v. 1, no. 3), Berlin, 1932, p. 46-48. 


383(L].—J. C. P. MiLier, “Two numerical applications of Chebyshev 
polynomials,” R.S. Edinb., Proc., v. 62, no. 22, 1946, p. 204-210. 
17.5 X 25.7 cm. 

1. The strong convergence of an expansion in Chebyshev polynomials renders them 
useful for interpolation. Let f(a + ¢) be expanded into a series of C,(4t/h). The coefficients 
a, of this expansion are expressible in terms of the derivatives f(a). They are also ex- 
pressible in terms of central differences. The corresponding formulae are derived by opera- 
tional methods and numerical tables given. Even if high central differences are needed, the 
smaller number of a», terms makes interpolation more convenient, particularly in conjunc- 
tion with a table for the C,(x). A numerical example illustrates the advantage of the method. 

2. Since an expansion in Chebyshev polynomials is merely a modified form of a Fourier 
series, a table of the Chebyshev polynomials becomes useful for harmonic synthesis whenever 
the sum of a Fourier series is required for arguments which are convenient numbers in 
x = 2 cos @ rather than in @ itself. An application is given, showing how the evaluation of 
the Mathieu functions can be facilitated by this procedure. 

CorNELIUS LANCzOS 

12012 Renton Avenue 

Seattle 88, Wash. 


384(L].—F RANK H. SLAYMAKER, WILLARD F. MEEKER & LYNN L. MERRILL, 
“The directional characteristics of a free-edge disk mounted in a flat 
baffle or in a parabolic horn,” Acoustical Soc. Amer., Jn., v. 18, Oct. 
1946, p. 363-368. 19.4 X 26.6 cm. 


There are two tables on p. 367. T. I of 

Tiana) 
Ji(\na) 
r/a = 0(.1)1; Ana ~ 3.01, 6.21, 9.37 for m = 1(1)3; m = 1, 2 to 3S, and n = 3 mostly to 4S. 


b(Anr) = 





Jo(An?) + To(An?) 


G = J,(ka sin $)/ka sin ¢; 
Py = — [1(Ana)DnaJ1(Ana)Jo(ka sin @) — ka sin ¢Jo(Ana)J1(ka sin ) J 
+ [Ji(Ana)(An2a? — Ra? sin? @)], if Wn + R sing 
= — 1,(na)[Jo2(Ana) + J12(Ana) J/[2J1(Ana)], if An = R sin ¢; 
Qn = [rnali(Ana)Jo(ka sin ¢) + ka sin @Jo(Ana)J1(ka sin ) ] 


+ (An2a? + ka? sin? ¢). 
T. II is of G, Pa + Qa{n = 1(1)3], mostly to 3S, for ka sin @ = 0(.5)10(1)12. There are 
graphs of P: + Q:, P2 + Qo, Ps + Qs on p. 368. 
R. C. A. 
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385[L].—FRraAnceEsco Tricot, ‘‘Generalizzazione di una formula asintotica 
sui polinomi di Laguerre e sue applicazioni,’’ Accad. delle Scienze di 
Torino, Cl. d. sci. fis., mat., e nat., Atti, v. 76, 1941, p. 288-316. 16.6 XK 25 
cm. 


The Tricomi polynomials of Laguerre! are defined by 
L(t) = (n!)“etd"(e—t") /di" = M(— n, 1, t) 


2 — 1) 
(- 1)"(n!) [ _ oe + ee — eee t+ (- 1yent |- 


(Daeg te og 


This is the case where a = 0 in the more general formula 


Liza(t) = (n!)ett-2d"(e—*t"**) /di" = M(— n,a + 1,2). 


In Tricomi’s paper are the following: 
p. 292, values of e~#Zio(t), for ¢ = [.5(.5)3(1)8(2)34; 5D], A; 
p. 302, a 4D table of the roots of the equations x + sin x = a, for a = 0(.1)3(.02)3.18; 
p. 303, graph of zeros of L(t) for m = 1(1)10; 
p. 315-316, table of e“#Z,(t), m = 1(1)10, ¢ = [.1(.1)1(.25)3(.5)6(1)14(2)34; 4D]. 
The polynomials L,(¢) and their properties were given in E. N. LAGUERRE, “Sur l’intégral 
: e-*dx/x," Bull. Sci. Math., v. 7, 1879, p. 72-81; Oeuvres de Laguerre, Paris, v. 1, 1898, 


p. 428f (L,(t), p. 430). L%4a(t) seems to have been discussed simultaneously with L,(¢), by 
N. SoNIN, in a memoir dated Aug. 1879, Math. Annalen, v. 16, 1880, p. 41 (function T,,*). 
The first one to refer to L%,.(¢t) as generalized Laguerre polynomials appears to have been 
another Russian, WERA MyLLER-LEBEDEFF, Math. Annalen, v. 64, 1927, p. 410. 

The so-called polynomials of Laguerre were introduced into mathematical analysis by 
LAGRANGE,? more than 130 years earlier than Laguerre, in his solution of a dynamical 
problem in which the oscillations of a vertical chain are represented approximately by those 
of a set of similar weights equally spaced on a light string. (See H. BATEMAN, “Lagrange’s 
compound pendulum,” Amer. Math. Mo., v. 38, 1931, p. 1-8.) The polynomials were also 
considered by ABEL, in 1826, “Sur une espéce particuliére de fonctions entiéres nées du 
développement de la fonction (1 — v)~e~*"/@-») suivant les puissances de v.”’ * This function 
is equal to >, Ly(x)v*/k!. In H. BeTHE, “Quantenmechanik der Ein- und Zwei-Elektronen- 

Kk 


probleme,” Handbuch der Physik, second ed., v. 241, Berlin, 1933, p. 289, this result is 
attributed to E. ScHRGDINGER,‘ just 100 years later. 

Laguerre polynomials are also of use in (a) the theory of hydrogen-like atoms; (b) the 
problem of numerical integration over the range 0 to + « [Gauss’s method with Legendre 
polynomials for ranges —1 to +1, or 0 to 1; Hermite’s polynomials for the range — ~ to 
+ 1]; (c) the discussion of the mathematical foundations of the electromagnetic theory of 
the paraboloidal reflector.’ In Bateman’s Bibliography® there are 47 references for Ln(x), 
and 73 for “L,*(x).” 

In preparing this RMT I have been indebted for some assistance from Dr. J. C. P. 
MILLER, and from Dr. ALAN FLETCHER. 


R. ¢. A. 


; 1 *< MTAC, v. 1, p. 361, 425; and v. 2, p. 31 [where L,,(x) is defined without the factor 
(n!)-*], 89. 

2]. L. LaGrancE, “Solution de différents problémes de calcul intégral,” Miscellanea 
Taurinensia, v. 3, 1762-65; Oeuvres, v. 1, Paris, 1867, ‘‘Des oscillations d'un fil fixe par une 
de ses extrémités, et chargé d’un nombre quelconque de poids,” p. 534-536; there are four 
of the polynomials on p. 536. 

3.N. H. ABEL, Oeuvres Completes, Christiania, 1881, v. 2, p. 284. 

4 E. ScHRODINGER, Annalen d. Physik, v. 385, 1926, . 485. 

5 PInNEY, Jn. Math. Phys., v. 25, 1946, p. 49f. Harry Bateman’s Bibliography, p. 77-79. 











268 RECENT MATHEMATICAL TABLES 


386[L, M].—S. A. Kuristranovicn, S. G. MIKHLIN & B. B. Davison, 
Nekotorye novye voprosy mekhaniki sploshnot sredy [Some new questions 
in mechanics of a continuous medium]. Moscow and Leningrad, Akad. 
N., Matematicheskii Institut imeni V. A. Steklova, 1938. 407 p. 16.7 x 
25.3 cm. 

We shall list certain tables in this volume p. 274-336 of a section written by Davison, 
and p. 392-395 of an appendix to the work, presumably written by the joint authors. 
On p. 274 a table is given of 
kax/q = tan-?Ve™le — 1 — (1/p)Vetla — 1 
for kry/q = [0(.2)3; 5S], 1/u = 1(1)4. 


ir do K(1 — #) 
K(t) = —————._, y(t) = tan~! ————-. 
' JS 1 — ¢tsin? @ : , K(t) 
On p. 334 there is a table of $(¢), to 3D, for £ = 0(.00001 ).0001 (.0001 ).0005(.0002).0015, 
.002(.001).01(.003).016, .02(.005).05(.01).1(.02).2(.1).5. 
On p. 335-336 are the following five tables: 


(a) z = fe, (b) e-8#/*, to. 3-45, 


for —dA = .01, .1, .2, .5(.1).7, 1(1)3, 5, 8, 15, 20, and A = 1.05, 1.1, 2(1)8, 12, 15, 20. 


; sal 1 v(t)dt —3w (25 

©, @w=R[ PZ], jewer, 
mostly to 3-4S, for \ = .00001, .005, .01, .05(.05).3, .4(.05).6, .7(.05).95, .99, .995, .9999, 
.99999. In (a) and (d) approximations are given to true results with possible maximum devia- 
tions. Before discussing the remaining tables in the volume we may quote some results from 
B. A. BAKHMETEV, Hydraulics of Open Channels, New York, 1932, p. xv—xvi, 308-311: 
If y is the depth or stage of flow, yo the normal depth of flow or the depth of flow in uniform 





movement, and 7 = y/yo, then the varied flow function B(n) = — _ 1’ 
7” — 


the hydraulic exponent. There are tables of B(n) (f) for » > 1, (g) for 7 < 1, m = 2.8(.2) 
4.2(.4)5.4 
(f) » = 1.001, 1.005(.005)1.02(.01)1.2(.02)1.5(.05)2(.1)3(.5)5(1)10, 20 
(g) » = 0(.02).6(.01).97(.005).995, .999. 
Now on p. 392-395 are tables of 


n being 








ni @-) nk-4 


(h) A(n) = fT an (k) Cla) = fe ao, 


k = 2n + 1, fork = 3(m = 1), 4(m = 3/2), » = [0(.05).6(.01).9(.005)1.05(.01)1.5(.05)2(.1) 
3(.5)5, 6(2)10; 4D], A. 


By trial we found that in the table of A(m), when k = 3, no = 0, for 7 < 1, and 9 = © 


for 7 > 1. We were unable to determine the values of 7 leading to the table of C(n), or to 
that of A(n), when k = 4. 


S.A J.& RCA. 


387[L, M].—NBSMTP, ‘‘Table of the Struve functions L,(x) and H,(x),” 
Jn. Math. Phys., v. 25, Oct. 1946, p. 252-259. 17.3 & 25.3 cm. This is 
Applied Mathematics Panel, Report 59.1, referred to in MTAC, v. 2, 
p. 39. 

The functions in question are, in Watson’s notation, in effect, 


La(z) = 3 an(z); Hale) = ¥ (— 1)"¢a(2), 
m=0 


m=0 





where 


Also 
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where 
Om(x) = ($x) /P(m + 3)0(m + n + 9). 
Also 


L(x) = N Sf ¥ sinh(x cos @)sin*0d0,  N = 2(4x)*/T(n + 4)r(4), 


H,(x) = N a ¥ .in(x cos @)sin6d0, provided that R(x) > — }. 
But for all x 
Lyx) = An(x) + Balx),  Ha(x) = An(x) — Ba(x), 


where 


An(x) = 2) dam(x), Bale) = 2 amsa(z). 
m=0 m=0 

Tables are given of L,(x) and H,(x), m = 0, — 1, — 2, x = [0(.1)10; 7-10S]. Lo(x), L_s(x), 

are also with & and *5‘ modified throughout the range; L_2(x) with & and *8* for x =2.1(.1)10 

For x = 0(.1)2, L_2(x), xL_2(x), &(xL_2), *6*(xL_2) are given. In the case of Ho(x) and 

H_i(x), and H_2(x), for x = 2(.1)10, *8? is given; H_2(x) for x = 0(.1)2, also xH_2(x), has *&. 


Hi(x) = 2/x 7- H_i(x), Li(x) = = 2/x + L_i(x). 


Watson! has tabulated Ho(x) and Hi(x) x = [0(.02)16; 7D], no A; hence linear in- 
terpolation is here correct to about 5D. H_:(x) is readily obtained from Hi(x), and has 
been tabulated before by Arrey,? for x = [0(.02)16;6D], and by JAHNKE & Empe,? for 
x = [0(.01)14.99; 4D]. H_2(x), L_:(x), and L_2(x) seem to be here independently tabulated 
for the first time. 

Kar_ HERMANN STRUVE investigated‘ only the special functions Ho(x) and Hi(x), but 
properties of the general function were later extensively developed by StEMoN® and WALKER.* . 
The function L,(x) bears the same relation to Struve’s function H,(x), as I,(x) bears to 
In(x); Ln(x) = i-*H,, (ix), Lo(x) = — iHolix); Lo'(x) = 2/e — Hilix). Tables of Lo(x) 
and Lo’ (x) for x = [.02, .1, .5, 1(1)12; 6-7S] were given by OWEN.” H,(x) = (— 1)"*"L,(éx), 
and 


H, (xi!) = iL, (xi) = ster,ax + 7 steinx, 


a notation due to McLacuLtan & Meyers® (see MTAC, v. 1, p. 252, 460). For tables of 
y = 4n[Io(x) — Lo(x)], and — y’ = 1 — 4x[Ji(x) — Li(x)] by R. Zurmtw and R. 
Miter, see MTAC, v. 2, p. 59; and on p. 39 a table by Great Britain, Nautical Almanac 
Office, of f(x) = mwe~*(2x)7[Ki(x) + Lu(x) — 2/r]. 

When 1 is half an odd positive integer H,,(x) is expressible in terms of elementary func- 
tions. For example, Hy(x) = B(1 — cosx) = B — J_4(x), where B = [2/(xx)]}. For vari- 
ous tables of J;(x) = H_4(x), see MTAC, v. 1, p. 233. 

mn. © 


1G. N. Watson, A Treatise on the Theory of Bessel Functions, p. 328-329; tables, p 
97. 


2J. R. Arrey, BAAS, Report, 1924, p. 280f, — H_:(x) is also tabulated here for the 
same range. 

3 JAHNKE & EmpE, Tables of Functions, fourth ed., New York, 1945, p. 219f. Ho(x) is 
also tabulated here for the same range, p. 212f, and 218f. There are also Ges other tables 
of Ho(x) and H,(x), to 4D, by S. P. GLazenap, Matematicheskie i Astronomicheskie Tablitsy, 
Leningrad, 1932, p. 110f, x = 0(.02)16, an abridgment of Watson; and by N. W. Mc- 
LACHLAN, Bessel Functions for . Engineers, 1934, p. 176, x = 0(.1)15. 9. 

‘See MTAC, v. 1, p. 305. 

oP. Pans Ueber die Integrale einer nicht homogenen Differentialgleichung zweiter 
Ordnung. Progr. ‘Luisenschule. Berlin, 1890; see Jahrb. Fort. d. Math., 1890, p. 340f. 

6 J. WALKER, The Analytical Theory of ‘Light, Cambridge, 1904, p . 392f. 

'S. P. Owen, “Table of values of the integral o* Ko(¢)dt,” Phil. Mag., s. 6, v. 47, 
1924, p. 736; see also MTAC, v. 1, p. 245, 247, 301. 

; also N. W. McLacuian & A. L. MEYERS, (a) “The ster and stei functions” 
(b) “Integrals involving Bessel and Struve functions,” Phil. Mag., s. 7, v. 21, 1936, p. 
425-436, 437-448. 
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388[L, M].—S. Sko em, ‘‘En del bestemte integraler av formen f f(x)cosakdx 
0 





os f f(x)sin(ax)dx,’’ Norsk Matem. Tids., v. 27, 1945, p. 65-75; 
0 f 
tables, p. 70-71. 15.5 < 23.2 cm. 


S(a, 1) iy er dx = }[e-*Ki(a) — eEi(— a)] 


cosh a shia — sinh a chia 
1/a + 2!/a? + 4!/a5 + --- + (2n)!/a* 
S’'(a, 1) = — $[e*Ei(a) + e*Ei(— a)] = sinha shia — cosha chia i 
= — [1/a? + 3!/a' + 5!/a* + ---]. 
T. I, S(@, 1), for a = [0(.01).1(.1)1(1)10(10)100; 5D]; maximum value at a ~ .8791 is | 
approximately .64996. 
T. II, S’(a, 1), for the same range of a as in T. I; zero value at a ~ .8791, and minimum value 
at a ~ 1.8594 is approximately — .15583. 





389[L, M].—Epmunp C. Stoner, ‘“The demagnetizing factors for ellipsoids,”’ 
Phil. Mag., s. 7, v. 36, Dec. 1945 (publ. Sept. 1946; note added in proof 
28 May 1946), p. 803-821. 17 & 25.4 cm. 
J. A. OsBorn, ‘“‘Demagnetizing factors of the general ellipsoid,”’ Phys. 
Rev., v. 67, 1945, p. 351-357. 19.2 & 26 cm. 


The formulae for the demagnetizing factors, in terms of F and E, are 





abc 
Da = L/4er = F(k, — E(k, 
[Ae = a a Lie) — Ble, 99] 
cos 6 cos } 


=~ —————- [PB 6) — E(t, 6)] 


sin’ @ sin? a 











abc 
= =- F(k, ¢) — E(k, 
Dy = M/4x @-eai@ — we L (k, ¢) — Elk, ¢)] 
+ —— ze, 6) - 
(a? — c2)4(b? — 2) @ aw 6 
sei = cos } [ E(h, 4) — cov oF lh, ¢) — sin? @ sin ¢ cos ] 
sin? @ sin? a cos? a cos 6 
abc 
D. = N/4x = — E(k, 
[he = — age a BO 
b cos @ cos i C7) 
+4 a one [meee — E(k, ¢) |, 
B—c sin’ ¢cos?a cos } 


where k? = (a? — b?)/(a? — c*) = sin? a, cos 6 = b/a, cos = c/a; a, b, cla 2 b 2 c) are 
the semi-axes of the ellipsoid. D, + D, + D. = 1. 

Consider first, ellipsoids of revolution: a polar semi-axis, b equatorial semi-axis, m = a/b, 
uw = b/a; m <1 and yw > 1 an oblate spheroid; m > 1 and uw < 1 a prolate spheroid. Then 


co ds 1 m 
o> .é - cosh! ba ’ 1, 
D. a? f° @+5\% +5 @F—-) [ (nm? — 1) cosh m — 1 ] m > 





1 m - 
“a Sal!- aoa m|, m <1; 
D, = 4(1 = Da). 
Stoner gives tables (p. 816-817) of Da, m or pw = [0(.1)5(.5)10(1)25(5)50(10)150(50)400 
(100)1300; 6D]. In general the m-table will be appropriate for prolate spheroids and the 
u-table for oblate spheroids. To ensure accuracy to the sixth place, the calculations were 


carried out so as to give unit accuracy in the seventh place, and rounded six-place values 
are presented in the tables. 
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Osborn gives two tables (p. 353-354) of demagnetizing factors of the general ellipsoid, 
L/4x, M/4x, N/4x, for (T. 1) cos ¢(= c/a), @ = 10°(10°)70°(5°)85°, 88°, 89°, and for 
cos @ (= b/a), @ = [0(10°)90°; 5D]. Also (T. II) cos @ = .1(.1)1, cos @ = various values. 
There are three large-scale graphs of L/4x, M/4x, N/4x, for 0 < c/a < 1, b/a = 0(.1)1. 
In T. II the values are accurate to 3D and are probably in error several units in the fourth 


| place. 


Stoner has a single graph of these same functions for b/a = .2(.2)1. 


nt. A. 


| 390[L, P].—N. W. McLacuian, Bessel Functions for Engineers, Oxford 





Univ. Press, London, Geoffrey Cumberleye, 1946. xii, 192 p. Reprinted 
photographically. 15.3 23.3 cm. 18 shillings. 


This very useful volume of the Oxford Engineering Science Series first appeared in 1934. 
In the corrected photographic reprint published in London in 1941, two pages were added; 
the new material included an introductory ‘‘Note,” a page of “‘Additional formulae,” and 
20 (instead of 6) “Additional references.” In the Note it is remarked that “The omission 
of contour integral representation of Bessel functions and its technical applications has 
been rectified through publication [by the author] in 1939 of Complex Variable and Opera- 
tional Calculus with Technical Applications,” and a correction of an error on p. 300 of this 
volume is noted. 

We have already referred to various tables in the volume under review (see MTAC, v. 1, 
p. 212, 216, 220, 246, 247, 254, 255, 257, 258, 297). In the right-hand member of formula 147, 
p. 167 one error:still persists; the sign — should be changed to +. The 1946 edition is an 
exact reprint of that of 1941, except for the correction of four signs, two in each of the lines 
— 3and — 5, p. xi, ‘Additional Formulae.” 


R. C. A. 


391[L, S].—C. Srracney & P. J. WALLIs, ‘“Hahn’s functions S,,(a@) and 
Un(a),”” Phil. Mag., s. 7, v. 37, Feb. 1946 [publ. Nov. 1946], p. 87-94. 
16.8 X 15.1 cm. 


“In a paper! on the calculation of fields in certain resonators, Hahn introduced two 
new functions: 
> m? sin? nra 


~ Sale) = elm — wat) 


and 


> a?mn? sin? nra 


Un(a) = with 0 <a <1, 


nai (m? — n°o?)? “i 
and used these functions to shorten his calculations. Since this time, Hahn’s method has 
been used for certain similarly-shaped resonators and Hahn’s two functions usually help 
to shorten the solution considerably. Hahn himself only gave a small table of S,,(a@) and a 
few values of Um(a).3 

“In this report closed expressions are derived for the case of a rational, and are used 
to produce a much more comprehensive table of S,(a) and a slightly smaller table of 
Um(a)/m. In a concluding section integral expressions, power series in a, and asymptotic 
series in m are given which together facilitate the calculation for values of a not given in 
the tables.” 


Tables: — Sn(a), for m = 1(1)10, a2 = [0(.1)1, .25, .75, 3, ¥; 5D]; 
Um(a)/m, for m = 1(1)10, a = [0(.25)1, 4, 3; SD]. 
Extracts from text 


1W. C. Haun, “A new method for the ap of cavity resonators,” Jn. Appl. Phys., 
v. 12, 1941, p. 62-68. There are 2D values of — Sn(a) for m = 1(1)9,a = oiceeay 1: 
also — So(a) for a = 0(4)4, x; and of DD, m = 1(1)4. See MTAC, RMT 208 and 
MTE 69, v. 1, p. 425, 451.—Epirors. 
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392[M].—NATIONAL RESEARCH CouNcIL oF CANADA, Division of Atomic 


Energy. Report no. MT-1 dated Chalk River, Ontario, December 2, | ned 
1946, The Functions E,(x) = f e~=“u-"du, 39 leaves mimeographed on | difference 
1 | 


one side, with covers. Introductory material, p. 1-7 by G. PLAczEk; 
Appendix A, an ‘asymptolic expansion for E,(x), by Dr. GERTRUDE | 395[{N]- 
BLANCH, p. 8; Tables, by NBSMTP, p. 9-39. 20.3 & 27.4 cm. This | 
edition contains corrections of one which appeared in July-August 1946. 


navi 
Witt 
The functions E,(x) play an important role in diffusion theory. The discussion of | PatmovI 
certain integral equations can be simplified by their use; expansions in terms of these func- | and give 
tions are also often found convenient for the numerical evaluation of integrals occurring in | 2), for i 
connection with transport problems. The functions have been defined by SCHLOMILCH,! and | 1 
have been extensively used by SCHWARZSCHILD,? Epp1ncTon,* Hopr,‘ and others. In spite ay ~— 
of this no systematic effort for their tabulation seems to have been made up to the present. | p 
An attempt by Mran & Cuapman' to approximate the functions by “index sums” was not 


| ap~1/ 
| P 
accurate enough for our purposes. 
E,(x) is here tabulated, for nm = 0(1)20, x = [0(.01)2; 7D], [2(.1)10; 7-10D]]. On p. 39 | 
are tables of E2(x) — x In x, for x = [0(.01).5; 7D], and of E3(x) + 3x? In x for x = [0(.01).1; 396(Q] 
7D], for use in interpolation. Since E:(x) = — Ei(— x) = JS° e~“u-du, there are extensive las 
tables of this function in NBSMTP, Tables of Sine, Cosine, and Exponential Integrals, v. 1-2, tig: 
1940, for x = [0(.0001)2; 9D], [0(.001)10; 9S]; [10(.1)15; 14D]. a 
Extracts from introductory text —s 
EpitoriaL Notes: In FMR, Index. p. 207, are given details of 8 tables, E_n(x) | 
= fi? etundu = x) £,° e“u"du, five of them including negative values of n. Eo(x) Cor 


= e*/x, of which values for x = [.1(.001)1(.01)2; 9D] are given by W. L. Miter & T. length « 
R. RosesruGu, “Numerical values of certain functions involving e~*,” R. Soc. Canada, | a unit 
Trans., s. 2, v. 9, 1903, sect. III, p. 102-107. See also TAKEO AKAHIRA, “Tables of e*/x hn 
and /-° e“du/u, from x = 20 tox = 50,” Inst. Phys. Chem. Research, Tokyo, Sci. Papers, en os 
Table no. 3, 1929, p. 180-215; the interval of the table is .02, to 5-6S, A*. 0 to 36( 


10. Scuiémiicu, ‘Ueber Facultatreihen,” Z. Math. u. Phys., v. 4, 1859, p. 390f. 

2 K. SCHWARZSCHILD, (a) ‘“‘Ueber das Gleichgewicht der Sonnenatmosphare,”’ Gesell. d. 
Wissen., Géttingen, Nach., ren Kl., 1906: p. 41f; (b) “Uber Diffusion und Absorption 
in der Sonnenatmosphiare,” Akad. d. Wissen., Berlin, Sitzb., 1914, p. 1183f. 397[U 

*A. S, Eppincton, The Internal Constitution of the Stars, Cambridge, 1926, p. 333. 











*E. Hopr, Mathematical Problems of Radiative Equilibrium (Cambridge Tracts .. ., _ 
no. 31), 1934, p. 21, etc. | Ja 
5 A. M. Mian & S. Cuapman, “Approximate formulae for functions expressed as definite | th 
integrals,” Phil. Mag., s. 7, v. 33, 1942, p. 115f. It is noted that E,(x) arises in the theory M 

of absorption of radiation in an exponential atmosphere. There is a table on p. 119 of ap- { 
proximate values of E,(x), for m = 2(1)8, for x = 0(.5)3(1)6, also .01, .05, .1, .25. TI 
: aroun 
393[M].—W. Soxo.ovsky, “Plastic plane stressed states according to | eight 
Mises,”’ Akad. N., USSR, Leningrad, (Dok.), C. R., n.s. 1946, v. 51, 23, 18 
p. 177. 16.8 K 26 cm. | age of 
There is here a table of | . ¥ 
) Rio d 

1 R(t)dt 

—Q(x) = = - | in 191 
2/%e sint of the 
= 4 — sin“ (2 cos x/34) + } tan [(4 cos x + 3)/R(x)] + } tan [(4 cos x — 3)/R(x)], wines 
where R(x) = (3-4 cos? x), for x = 3x to $x, mostly at interval gsx, to 3D. the in 

| 
. . to be 
394[M].—A. J. C. Witson, “‘The integral breadths of Debye-Scherrer lines ot ad 
produced by divergent X rays,”” Phys. Soc., London, Proc., v. 58, July 7 
1946, p. 407. 18 X 26 cm. | table: 
There is given here a table of D(u) = 4u-*Jo"[C2(u) + S?(u) udu, for u = [0(.1)5;4D], ident 


5D for u <1. For u < 2 the values were calculated from the series for D(u), those for of th 
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tomic 


er2| 4> 2 by numerical integration of four-place tables of C(u) and S(u). In the range .5 to 2 


the greatest difference between the values calculated by the two methods is .0003 ; the mean 
~d on | difference is about .0001. 





Extracts from text 
“ZEK; | 

RUDE | 395[N].—Ericu Micnavup, “Beitrag zur Amortisationsrechnung,” Skandi- 
This | navisk Aktuarietidskrift, 1946, p. 80-84. 15.5 K 23.5 cm. 








ions: With references to earlier discussions by E. LinpeL6r, K. A. Pouxka, A. Bercer, R. 
ion of | Parmevist, H. Hotme, and E. Franckx, the author considers the following five formulae 
fune- | and gives tables for each of them to 7D for half-yearly, quarterly, monthly rates (p = 2, 4, 
ing 7 | 12), fori = 1%(1%)9%: 
and | 
spite | o~z(1-2=*1), op ~3 (1 - P= 4 CAN 9), 
esent. | p 2p p 2p 6p? 
is not | = p-—1., x eee ‘ ae. 6p + 1(1 + p) ] 
ap (9+ 2 ‘), ap = [(1 + 2)¥? — 17/1, ap + 25 | 
p. 39 | 
sorte | 396[(Q].—ENnRIQUE VIDAL ABASCAL, El Problema de la Orbita Aparente en 
rey las Estrellas Dobles Visuales: Diss. Spain, Consejo Superior de Inves- 
isi tigaciones Cientificas, Instituto Nacional de Geofisica, no. 6, Observa- 
- torio de Santiago, Publicaciones, 11, Santiago de Compostela, 1944. xvi, 
ste) 62 p. 21.2 K 27 cm. 
Eo(x) Consider ellipses with common major semi-axis, OA = 1, and eccentricities e = the 
& T. length of OF; = .1(.1).9; then the foci Fi, i = 1(1)9, divide OA into tenths. Suppose that 
a ' aunit circle, with center at O, has been drawn, and P is any point of the circumference, 
pers, then F;P and F;A are the sides of circular sectors, F;A P Fi, whose angle a may increase from 
0 to 360°. A table, p. 53-62, gives the area of such sectors, to4D, e = .1(.1).9, a = 0(1°)360°. 
I. d. R. C. A. 
ytion 
397[U].—FRANcIscoO RADLER DE AQuINo, “Universal’’ Nautical and Aero- 
nautical Tables. Uniform and Universal Solutions Ultra-simplified. Rio de 
i a Janeiro, Imprensa Naval, 1943, 18, 247 p. 17.5 X 24.5 cm. Copies of 
nite this volume may be had from Weems System of Navigation, Annapolis, 
on} Md. at $9.00. 
The author of these tables, a captain in the Brazilian Navy, is well known to navigators 
| around the world, having published more than fifty papers on navigation in the past forty- 
to eight years. Not so well known is the fact that he was born in New York City on January 
51, 23, 1878; his mother was an American, his father a Brazilian. He moved to Brazil at the 
age of 13 and entered the Brazilian Naval Academy at 15. 
This volume is the second Brazilian edition of a book which was first published in 
| Rio de Janeiro in 1903. Editions were published in London in 1910, 1912 with reprintings 
| in 1917 and 1918, and in 1924; and in Annapolis, Md. in 1927 and 1938. The title and content 
«)] of the tables have changed slightly from edition to edition. For those familiar with the 
: earlier editions, it may be said that the principal change in this edition is the reduction of 
| the interval of the argument, latitude, from 1° to 10’. The method and the tables continue 
_ to be universal in that they allow the determination of the altitude and azimuth whatever 
the values of latitude, hour angle, declination and altitude. 
aly The first eighteen pages in this volume include the title page and explanation of the 
tables in English; the next sixteen pages (numbered 1 to 16 also) present similar but not 
Dj, identical material in Portuguese. The principal table was designed to be used in a solution 


for | of the astronomical triangle in which a perpendicular is dropped from the zenith upon the 
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hour circle through the celestial body. The length of the perpendicular and the declination 
of its foot, are called a and b respectively. The angle at the zenith between the perpendicular 
and the meridian (toward the elevated pole) is called a; that between the perpendicular and 
the great circle from the zenith to the celestial body is called 8. L, t, and d denote the latitude 
of the observer, the local hour angle and the declination of the celestial body respectively. 
A is the angle of the astronomical triangle at the celestial body. C is the angular distance 
from the celestial body to the foot of the perpendicular. 

The basic equations for the solution of the two right triangles are obtained by Napier’s 
rules; they are: 


csc a = sec Lesci?; tan b = tan L sect 
tana = csc L ctné; csc h = seca sec C 
tan A = tanacsc C; tan 8 = tan Cesca 


The rules necessary for the use of the equations.and the tables are given on the pages 
numbered 8 in the explanations in English and Portuguese. At first sight, they appear to 
be of the same order of complexity as those in H.O. 208, Dreisonstok (RMT 103), but use 
proves them to be somewhat simpler. 

The principal table (p. 36-215) has as vertical argument the local hour angle, t, 0(10)90°, 
and as horizontal argument the latitude, L, 0(10’)89° 50’. The tabulated quantities are a 
and 6 as defined above, each to the nearest minute of arc, and a. The values of a are given 
in heavy type to distinguish them from those of 6. On the right-hand side of the page, a is 
given to the nearest tenth of a degree for each degree of local hour angle and for the middle 
of the degree of latitude. 

One enters the table with the dead-reckoning latitude rounded off to the nearest ten 
minutes of arc and with the local hour angle of the body to the nearest degree; a and b are 
copied out for these arguments and a is taken for the nearest half degree of latitude. C is 
formed using the equation, C =|b — d|. One re-enters the table looking for a (rounded 
off to the nearest 10’) at the bottom of the page and C (to the nearest degree) along the right- 
hand side of the page. With these arguments, the dark-faced column yields the altitude, h, 
and £ is found in the right-hand column opposite it. The azimuth angle is found by adding 
a and £. 

In the explanation, the author indicates that the determination of 8 by this method 
is weak and proceeds to give three other methods of finding it. The first two involve the 
substitution of arguments; one can look for 90° — C in the 8 column and interpolate the 
value of 8 in the left-hand column above B, or one can interchange the values of a and C 
and interpolate the value of 8 in the column footed A. The third method is to use the last 
equation above with a table of log tangents and log secants. Actually the table given is one 
of values of 10° log tan x and 10° log sec x to the nearest integer for argument x, 0(1’)89° 59’. 

To allow for the minutes discarded in a and C, Aquino suggests the use of: 


Ah, = Aa cos B, Ahe = AC cos A. 


He provides a ‘‘difference of latitude and departure” table to simplify their use; the table 
will also be useful in dead reckoning. To allow for the minutes of latitude discarded, another 
table is provided on the inside of the back cover and the page facing it. This same table is 
offered on the two sides of a separate sheet of cardboard and again on one side of a separate 
folded sheet of heavy paper. The corrections, Ah; and Ahz, may be avoided by the use of 
the equation for csc h, and the log tangent — log secant table. 

The author states that the tables were computed by means of Callet’s and Bagay’s 
seven-place logarithms, with many values determined by Vlacq’s ten-place logarithms; 
Vega’s ten-place table, based on Vlacq, is much more accurate. He further states that if the 
declination be taken to the nearest tenth of a minute, and Ah; and Ah used, the maximum 
possible error in h will be 1.6’ but that the actual error in practise will hardly ever be over 
0.5’. If only Ahz is used, the altitude obtained is always within 5’ of the true calculated alti- 
tude. Although the volume under review is dated 1943 there is in it a yellow sheet dated 
9 April 1946, listing 32 corrigenda. 
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Other tables contained in the volume are four-place logarithms and antilogarithms with 
convenient proportional parts tables, distance to the horizon and dip of the horizon for 
different elevations, combined corrections for refraction, dip of the horizon and, where they 
are significant, semidiameter and parallax for planets and stars, upper and lower limb of 
the sun and lower limb of the moon, for altitudes 8° to 90° and for elevations 0 to 15 meters. 
A similar table for the upper limb of the moon would make a worth-while addition. A small 
auxiliary table allowing one to correct altitudes less than 8° for refraction is a valuable item, 
not often found in navigation tables and especially needed in the polar latitudes. 

A table that remains in use almost half a century while other tables come and go, can 
reasonably be said to have a strong appeal to the average navigator. To appreciate Aquino’s 
great contribution to navigation, one needs only to compare the first edition of this table 
which appeared in 1903 with other tables and methods then in use. It is unlikely that a 
person who has been trained in the use of H.O. 214 or H.O. 218 will change to this table, 
but a person who learned Aquino’s method first might continue to prefer it because of its 
beautiful simplicity, its universality and the small bulk of the tables. 


CHARLES H. SMILEY 
Brown University 


398[U].—W. MyerscoucH & W. HamiLton, Rapid Navigation Tables. 
London, Pitman, 1939. ii, 109 p. 16.4 X 26.6 cm. 10s. 6d. 


These tables are designed for the solution of the astronomical triangle for altitude, h, 
and azimuth, Z, when latitude, Z, hour angle, ¢, and declination, d, are known, and for 
similar problems. In the procedure for the determination of h, Myerscough & Hamilton 
follow such tables as those by SouILLAGOUVET, OGURA, WEEMS (RMT 315), DREISONSTOK 
(RMT 103), and HuGuHes-Comriz (RMT 115), each of which divides the astronomical 
triangle at the zenith into two right triangles. The first triangle is solved for tabular values 
of L and ¢ by a table which gives the remaining parts, as angles or as logarithmic functions 
of angles, without interpolation or other calculation. The second triangle is solved for h by 
logarithmic processes. The several tables differ only in notation and in that one of the 
auxiliary angles used in some of the tables is the complement of that used in the others. 

In the determination of Z the several tables cited show a pleasing variety in method. 
Dreisonstok and Comrie follow BERTIN in deriving the two component parts of Z from the 
same auxiliary triangles as are used in the determination of h. Souillagouet utilizes another 
division of the astronomical triangle in order to get Z in one piece. Weems uses the graphical 
“Rust diagram,” and in his New Line of Position Tables (RMT 315) provides also an in- 
teresting variation on the Bertin procedure. Myerscough and Hamilton, however, follow 
Ogura in using the equation, 


cot Z = cos L (tan dcsct — tan L cot t) 


The most interesting and original feature of Myerscough & Hamilton is the inclusion 
of all data in one table of 91 pages (0 to 90°). At the first entry the page is selected for i, 
and for the left-hand argument LZ or d, as the case may be, the following quantities are 
extracted: 


P = length of side of first auxiliary triangle opposité zenith, deg. & min., 
Q = 10° log sec (side of same triangle opposite pole), 
X = 10tanLcott, Y = 10 tandcsct. 


For the second entry P is combined with d according to rules typical of such tables to give 
a side of the second auxiliary triangle. The page being for the degrees of this argument, 
and the entry by the minutes (using the same figures as were previously used for L and d), 
the following quantity is taken from the R-column of the table: R = 105 log csc (P ~ d). 
For the third entry the R-column is searched for Q + R = 105 log csc h, and h is obtained 
by reading degrees at the top of the selected page and minutes in the left-hand column. The 
equivalence of this procedure to those of the other tables cited is easily recognized. For 











276 RECENT MATHEMATICAL TABLES 


the fourth entry into the tables the page is selected for the latitude, and the Z-column is 
searched for Y — X = 10 cot Z csc L, and the azimuth is read opposite the nearest value. 
The four necessary openings equal those required by the other tables cited, so that the 
prospective user must seek grounds for preference in the arrangement of the tables, which 
is entitled at least to study by other table makers. 

The table would be easier to use if it gave the four values of ¢ for which a given page 
is used and not merely the one in the first quadrant. It would be improved also if the L 
and d argument went from 0 to 90° instead of stopping at 70°. (With these changes and 
two other minor ones the Y-column might be used in the fourth step, in order to eliminate 
the Z-column.) Since the tables are entered with L in the fourth step, there would be nothing 
gained by rearranging the tables for entry with L in the first, as there is in Hughes-Comrie 
and the new Weems. No data on the accuracy of the table are available. 

The tables of Myerscough & Hamilton and of the other authors cited above seek to 
avoid interpolation in the second auxiliary triangle by the use of logarithmic trigonometric 
functions. While there may be some historical justification for such a treatment, it should 
be pointed out that H.O. 214 has accustomed navigators to interpolation. There is, accord- 
ingly, good reason to reexamine the possibilities of such methods as that of Bertin, in which 
both triangles are solved by a single table. The four openings of the various tables cited 
above are reduced to two, with two interpolations of about the same magnitude as those 
customary with H.O. 214 (RMT 399). The Bertin method is possible with the well-known 
Sea and Air Navigation Tables of Captain RADLER DE AQuINO (RMT 397), and with a new 
Spheric Tabulations of R. C. Dove, R. F. D. No. 1, Collegeville, Penna. 


SAMUEL HERRICK 
Department of Astronomy, 
Univ. of California 
Los Angeles 


399[U].—U. S. HyproGrapuic Orrice, Publication No. 214, Tables of Com- 
puted Altitude and Azimuth, Latitudes 80° to 89°, inclusive, Vol. 9. Wash- 
ington, D. C., U. S. Government Printing Office, 1946, 3, xxiv, 263 p. 
22.6 X 29.1 cm. This is the last of nine uniform v. of H.O. 214, each v. 
devoted to 10° of latitude. For sale by the Hydrographic Office and by 
the Superintendent of Documents, Washington, D. C., $2.25 per v.; 
foreign price, postage extra. 


This review will be limited primarily to a discussion of the differences between v. 9 
and the other 8 v. of H.O. 214 which were reviewed earlier (RMT 105). This volume, like 
_ the others, was prepared by the Work Projects Administration, (Philadelphia Project No. 
24831), and presumably is of the same order of accuracy (see v. 2, p. 182f). The interval 
of argument for hour angle is 1° as in the other v.; it might well be 2° or perhaps even 5°, 
save for the loss of uniformity, since the tabulated altitude and azimuth change slowly and 
in a relatively linear fashion. 

The description of the tables and their use is almost entirely new and occupies some 
ten pages more than that in the other v. The use of the pole as an assumed position is 
explained as well as the use of gnomonic, stereographic, azimuthal equidistant and inverse 
Mercator projections. A brief description of grid navigation is given. 

Two ways in which this volume could be improved may be mentioned. The computed 
altitudes might be carried down to 0° or at least to 2°, since in the polar regions, the sun, 
moon, and planets spend a considerable fraction of the time at altitudes less than 5°; there 
are blank spaces available for these data. The second change would be to replace the refrac- 
tion tables given in the front by others especially prepared for the conditions of temperature 
and barometric pressure commonly found in the polar regions. 


CHARLES H. SMILEY 
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MATHEMATICAL TABLES—ERRATA 


References have been made to Errata in ‘“‘A new approximation to r”’; 
RMT 360 (Ser), 367 (Ostrogradskii), 368 (Patz), 369 (Cunningham), 370 
(Vinogradov), 375 (Fisher & Yates, Croxton & Cowden), 376 (Goncharov), 
378 (Hastings & Piedem), 382 (McLachlan), 389 (Stoner), 390 (McLachlan), 
392 (N. R. C. Canada), 397 (Aquino); UMT 54 (Br. Standards Inst.), 
55 (DeMorgan); QR 28 (Gauss). 


99. H. BRANDENBURG, Sechsstellige trigonometrische Tafel, Leipzig, 1932, and 
Ann Arbor, 1945. 


This volume was described in MTAC, v. 1, p. 387f and errors are listed on p. 388, and 
in MTE 28, v. 1, p. 162. Two further errors have been brought to my notice, both by 
correspondents in Palestine. 

(1) Mr. B. Goussinsky, Superintendent of Surveys, points out that the asterisks in 
the entries for the tangents of 0° 27’ 40” and 0° 27’ 50”, p. 32, should be deleted. In the 
7-figure tables tan 0° 27’ 30” is given as 0.007 9996; this has been rounded off (correctly) 
to 0.008 000 when forming the 6-figure tables, thus changing a leading figure and removing 
the necessity for asterisks. I have verified that there are no other similar instances of 
this fault. 

(2) Mr. ALEXANDER Katz, of ,-rusalem, finds that on p. 23, cot 2° 46’ 32”, for 62694 
we must read 62684. This error has taken “evasive action” on several occasions. It occurs 
in both editions of Brandenburg’s 7-figure tables. It is marked in my 1923 edition, and was 
presumably communicated to Brandenburg but does not appear in an extensive list of 
corrections that he prepared in 1927. It still occurs in the 1931 edition, and is marked in 
my copy, but not included in the list that I gave in a review in The Observatory, v. 54, 1931, 
p. 301-302. It was, however, communicated to Brandenburg, and is included in a list of 
errors published by him in 1932. Nevertheless Brandenburg repeats the error in his 6-figure 
table in 1932, and my proof readers, who examined the table, failed to find it! 


Lz. 


EpriroriaL Note: Fora reference to still other Brandenburg errors, see MTAC, v. 2, 
p. 46f. In the 1923 edition of Brandenburg’s 7-figure table are a two-page “Druckfehler- 
Verzeichnis,” and a leaf headed “Berichtigungsbogen zum Uberkleben der Druckfehler.”’ 
In addition to these there i is bound with the Brown University copy of these tables pur- 
chased in 1929, a 4-page ‘“‘Nachtrag zum Druckfehler-Verzeichnis,” on differently colored 
paper, dated 15 March 1927, and listing a very large number of errors, and the names of 
nine individuals assisting in its compilation. 

In 1932 the publishers distributed a one-page “Verzeichnis der Druck-,und Form- 
fehler,” 18 errors, in the 1931 edition of Brandenburg’s 7-figure table. In this sheet, no. 17, 
for 510 lies 519, read 310 lies 319. 

The five errors in this edition of Brandenburg’s 7- figure table, listed by L. J. C. in 
The Observatory, are as follows: p. 95, diff. following sin 4° 59’ 50’, for 485, read 483; p. 111, 
cot 7° 39’ 30”, for 8870, read 88/1; p. 157, cot 15° 13’ 30”, for 2744, read 2743; p. 177, diff. 
following sin 18° 32’ 10”, for 461, read 459; p. 335, cot 44°53’ 0”, for 0808, read 0807. The 
first is not included in the 1932 “Verzeichnis.” 

To conclude the references in MTA C to all known errors in the 1931 and 1932 Branden- 
burg tables we may note that the value of e in each of these tables is given to 30D as, 
“e = 2.718 281 828 459 045 235 339 784 490 662 . . .”” which is correct to only 19D. The 
value is obtained as the sum to 30D of terms 1/n! plus unity, m = 1(1)20. Of the values of 
these 20 terms 4 are not rounded off, and 5 are so rounded. Underneath Brandenburg’s 
incorrect value of e is given Euler's correct value to 23D, taken from his Introductio in 
Analysin Infinitorum, v. 1, § 122. Quite astonishingly the author seems to suggest that he, 
rather than Euler, is correct. 


100. FMR, Index. See MTAC, v. 2, p. 13-18, 136, 178-181, 219-220. 


The notes made by S. A. J. and myself in MTE 89, p. 178-181, were mainly a budget 
of suggestions on minor questions for consideration in a possible new edition. Doctors 
Fletcher and Miller have communicated to us certain statements which should have equal 
publicity. 
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P. 35, 3.14 we suggested “27 dec. Thoman’”’ instead of ‘‘20 dec. Thoman.” Dr. F. writes: 
“As its title implies, Thoman’s book mostly gives 27 dec. But log factorials are to 20 
dec.—I have the book in front of me. One or two other tables in it are also to less than 
27 dec.” 

P. 51, 4.41, for 4.412, read 4.4121; for 4.413, read 4.4132. Dr. M. writes “Delete this item. 
Such an alteration would destroy the effect intended. A, and S, are essentially the 
same function, and the special relations for functions of these types make 2x + 1 and 
x — 4 or 2x — 1 essentially the same, or closely related, arguments. Thus, I regard 
4.4122 as being correctly included in the general heading S,. Likewise 4.4131 and 4.4132, 
since x? + x = (x + 1)? — (x + 1), and S, and V, are related.” 

P. 144. I had written, ‘One wonders at the omission in 8.4 of a reference to Legendre’s table 
of log tan (45° + 4$x).” But there is no such table by Legendre; the table is of In tan 
(45° + 4x), which FMR carefully list on p. 184, as Dr. F. has pointed out. Here is a 
case where an index of names in Part I of the Index would have saved me from this 
slip. While my admiration for the Index has increased rather than diminished since I 
wrote RMT 233, I have more than once found difficulty in determining what printed 
or ms. tables by a given author are referred to in the volume. Hence I have started 
the preparation of a complete card catalogue of names mentioned up to p. 373 of the 
Index. It will only be after its completion that I can feel that all the resources of the 
volume are at my command. Perhaps these remarks may suggest to the authors some 
amplification in a second edition. 

P. 200, We wrote “‘7 dec. Brownlee 1923 (Russell, which one?),” since two Russells are 
mentioned in Part II. Dr. M. comments: ‘‘We do not pretend to give all computers of 
tables, listed under other authorship in Part II. The answer to ‘which one?’ here is, 
the one mentioned in Brownlee 1923.” This remark concerning FMR editorial policy 
would also make unnecessary the reference to L. Brockway under Sherman 1933, p. 433. 

P. 377. I raised the question of a possible edition of Bertrand’s Calcul d. Prob. in 1888. 
Dr. F. reports that such an edition is in the library of the University of Liverpool. 

P. 420, I noted, ‘The title of Newton’s work is Trigonometria Britanica (not Britannica).” 
Dr. F. writes: Spelling ‘Britannica’ was intentional as (i) it is the correct spelling of 
the Latin adjective in question, (ii) Newton (as Mr. Cosens pointed out to me) used 
it in headings, etc. There seems to be one m in the title (as you say) because there was 
no room for two. I dislike departing from the title page, but on this occasion decided, 
rightly or wrongly, to take a broad view.” 

P. 431. I had noted that of Schrén 1860, there was an Italian edition in 1867 and a French 
edition in 1891. Dr. F. remarks that at the University of Liverpool there is a Swedish 
edition, ed. by F. W. Hultman, Braunschweig, Vieweg, and Stockholm, Bonnier, 1868. 

P. 434. I raised the question of listing the 1819 edition of WILLIAM SPENCE, Mathematical 
Essays, either in addition to, or in place of this the listed 1820 edition. Dr. F. writes: 
“We have seen only Spence 1809. Very interested in your 1819 copy. I had four refer- 
ences to 1820: (i) English Catalogue of Books, ‘c. 1820’; (ii) Edinburgh Univ. Lib. Cat.; 
(iii) Univ. Coll.; London, Lib. Cat.; (iv) De Morgan, The Differential and Integral 
Calculus, London, 1842, p. 658-659. (ii)-(iv) give 1820 without ‘circa.’”’ I find that 
the Cat. of the Astor Library in New York also lists a copy of the 1820 edition. Thus 
four libraries having the 1820 edition at the times their catalogues were published, are 
listed. On the other hand the Boston Public Library and the Harvard University 
Library have copies dated 1819, ‘London, Printed for Thomas and George Underwood.” 
The Essays were edited by Sir John F. W. HERSCHEL! as a young man, and his “‘Preface 
to the Essays” is dated ‘Slough Dec. 10, 1818’; hence one would expect this volume 
to have been published in 1819. Since the above was written Dr. F. reported, ‘‘We have 
now seen the Royal Astronomical Society copy of that date.” 

nm. ©. A. 
1 HERSCHEL, CHARLES BABBAGE, and GEORGE Pzacock while still undergraduates 
formed in 1812 an Analytical Society, which through varied publications did important 


work in bringing about reforms in mathematical notation; see W. W. R. BALL, A History of 
the Study of Mathematics at Cambridge, Cambridge. 1889. ‘‘The Analytical School,” p. 117f. 
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101. NBSMTP, Table of Circular and Hyperbolic Tangents and Cotangents 
for Radian Argument, 1943. See MTAC, v. 1, p. 178f. 


For the value of tan 1.5708 read — 272241.80841, not — 27224.18084. 


J. A. MACDONALD 
53 Alexandra Ave., 
Canterbury, E. 7 
Victoria, Australia 


102. D. H. Lenmer, “On the converse of Fermat’s theorem,” Amer. Math. 
Mo., v. 43, 1936, p. 347-354. 
This paper contains a table (p. 349-351) of composite solutions m of the congruence 
2" = 2 (mod nm) having prime factors exceeding 313. A recent recomputation of this 
table by PouLet reveals the following complete list of errata: 


Delete, 68462551 5851 76839733 1019. 

Insert, 44070841 2113 74874869 3533 
70541099 4643 92438581 3331 
71079661 3187 96135601 881. 
74705401 3529 


A relatively unimportant error on p. 351 may be cited. Line 13 purports to give the 
product of all primes p where 17 < p < 101. This number contains the factor 79* and so 
should be replaced by 

775432448746244958042 1688873809769. 
D. H. L. 


103. J. T. Peters, Sechsstellige Tafel der trigonometrischen Funktionen . . . 
von zehn zu zehn Bogensekunden . . ., Berlin, 1929; there was a second 
edition in 1939. See MTAC, v. 1, p. 121, 162. 

The following additional errors have come to light when preparing copy for the new 

Chambers’ six-figure tables. All are on p. 8 of the 1929 edition. 

Cot 0° 27’ 3” for 127.086, read .086 


4 for 126.008, read 127.008 
5 for .930, read 126.930. 


In other words the integers 127 and 126 each need to be lowered one line. 


‘.# 


104. J. T. Peters, Zehnstellige Logarithmentafel. Erster Band: Zehnstellige 
Logarithmen der Zahlen von 1 bis 100000. Berlin, 1922. See MTAC, 
v. 1, p. 57-59. 
The proofs of a new six-figure table now being prepared for Messrs. Chambers were 
compared with this table, revealing, to our great surprise, two errors: 
P. 406 log 69731 for 843 4358 934, read 843 4258 934; 
p. 566 log 93748 for 974 9620 114, read 971 9620 114. 


This comparison would not reveal errors in decimals beyond the sixth. 

About 1924 I noted that Peters, in spite of the great pains he took to ensure complete 
accuracy in the tenth decimal (see his Einleitung, p. vii), and his list of cases in which 15 
or 16 decimals were necessary for this purpose, had missed one such case, 

P.16 log 11275 for 506, read 505; 


the differences are also affected. This error was pointed out as long ago as 1872 by GLAISHER 
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(R. A. S., Mo. Not., v. 32, p. 258—misprinted 358), who gives the mantissa as .05211 65505- 
49998 14---. In the following year Glaisher quoted a letter from the then owner of MICHAEL 
Tay.or’s copy of VLacg (ibid., v. 33, p. 452) saying that some previous owner (GARDINER 
is suggested) had corrected log 11275 by hand. 

Glaisher’s remarks on end-figure errors are quoted in N 72. They were prompted by 
this particular “error.” 


ie 


UNPUBLISHED MATHEMATICAL TABLES 


53[A, B]—NBSMTP, Tables of Circumferences and Areas of Circles. Tables 
prepared for the U. S. Bureau of Ordnance, Navy Department. Compare, 
MTAC, v. 2, p. 86-87. 


These tables are for circles with diameters ranging [.001(.001)10; 6D]. The computa- 
tions were made with IBM equipment, and a manuscript was prepared on the tabulator. 


NBSMTP 


54[L, M].—Carit Hammer, 7: able of selected values of Li(x) = f dt/Int, and 
0 


Zz 


dt/Int, mss. in possession of the author at 304 West 105th St., 
2 
New York City; and in the Library at Brown University. 


These tables are for x = [2(1)10(10)100(100)1000(1000)10 000(10 000)100 000; 8S]. 
The values for Li(x) were previously given by J. voN SOLDNER, Théorie et Tables d'une 
Nouvelle Fonction Transcendante, Munich, 1809, p. 43-49, for x = [0(.01)1(.1)2(.5)3(1)20; 
7D], [22(2)40(5)80(10) 160(20)320(40)640(80)1280; 8S], A® to .8. This table was reprinted 
(without A’, and with a misprint of 1220 for 1280) in A. DEMorGan, The Differential and 
Integral Calculus . . ., London, 1842, p. 662-663. Thus of 90 values given in the ms., 22 
were Soldner’s values. Three other values of Li(x), for x = 10%, 104, 105 were taken from F. 
W. BEssEL, “Untersuchung der durch das Integral / dx/In x ausgedriickten transcendenten 
Function,” Kénigsberger Archiv f. Naturw. u. Math., v. 1, 1811, p. 31, and F. W. BEssEL, 
Abhandlungen, Leipzig, v. 2, 1876, p. 339. The other 65 values were computed by means of 
NBSMTP, (a) Table of Natural Logarithms, v. 2, 1941; (b) Table of Sine and Cosine Inte- 
grals . . ., 1942; (c) Table of Sine, Cosine and Exponential Integrals, 2 v., 1940; (d) Tables 
of Lagrangian Interpolation Coefficients, 1944, using five points for x = 500 to 9000, and 
seven points for x = 20 000 to 90 000. 


C. HAMMER 


55[P].—SIDNEY JOHNSTON, Roller Chain Transmission Basic Exact Centre 
Distance Tables. Ms., iii + 10 sheets typed on one side. 20.3 X 32 cm. 
Original in possession of the author at 81 Fountain St., Manchester 2, 
England; carbon copy in the Library of Brown University. Among the 
“References” in the ms. are the following: (a) K. B. Jacos, “Driving 
chains and theory application to power transmission,’ Engineering and 
Shipbuilding Draughtsmen’s Assoc., Trans., 1928-29 (also as a pamphlet, 
T. 5, centre distance tables, p. 47-54); (b) Machinery’s Handbook, New 
York, Industrial Press, twelfth ed., 1943, p. 861-2. 


These tables are intended to serve the mechanical engineer in solving the bothersome 
problem of the design of roller chain transmission. Suppose that a roller chain of N links 
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and pitch P connects two chain gears, the number of whose teeth are T and ¢. The tables 
enable one to find the distance S between the centers of the gears. The formula on which 
the tables are based is: 
S = 4P(D — d)cscx 
where 
D = csc xT and d = csc at, 


and x is the least positive solution of 


cotx’+cx =a, where 
c=nr(T—t)/(D-—d) and a=[(N— }(T+H]/(D -—d) 


The table is one of double entry giving the factor csc x to 6, 5 and 4D in terms of ¢ and a. 
The value of c in a typical case is slightly greater than unity. The present tables are for 
c = 1, 1.01 and 1.02 only. The second variable a has the following range 


a = [1.75(.002)1.8(.0025)1.835(.005)1.93(.01)2.07, 2.075(.0125)2.325(.025)2.95(.05) 
4.15, 4.125(.125)6.5(.25)10; 6D], [11.25(1.25)22.5(2.5)35(5)105; 5D], 
[100(50)300(100) 1000; 4D]. 


The partial differences of csc x with respect to c and a are given together with the Bessel 
coefficient A(A — 1)/4. There is an auxiliary table of csc x7 for T = [10(1)160; 8D] for 
finding the pitch diameters d and D. Comparing this table with a similar smaller table to 
4D in British STANDARDS INSTITUTION, Specifications for Steel Roller Chains and Chain 
Wheels, revised April 1934, no. 228-1934, p. 18-19, one finds in the latter table, 10 last-figure 
unit errors: in excess for T = 44, 46, 63, 76, and in defect for T = 11, 22, 91, 95, 107, 133; 
also at T = 127, for 49.4295, read 40.4295. 

The present tables, if made available to machine design people, should do much to 
replace the crude approximations usually resorted to in dealing with this comparatively 
precise problem. A few handbooks give the exact formula 

N=T+2r(T — #)(tanA — A) 
where 
A = arccos P(T — t)/2xS 
or an equivalent formula but the reviewer has not found any other tables for obtaining S 


directly. 
D. H. L. 


MECHANICAL AIDS TO COMPUTATION 


The reader is referred to the first two articles of this issue, dealing with 
MAC topics. 


28[Z].—Akad. N., SSSR, Moscow, IJzvestiia, Otdelenie Tekhnicheskikh Nauk, 
1946, no. 8, September, p. 1065-1200. 16.7 X 25.8 cm. Entirely in 
Russian. 


This issue is almost wholly devoted to material dealing with computation and computing 
mechanisms. The contents are as follows: 


N. G. Bruevicn, “The present state and problems of the theory of mechanism precision,” 
p. 1065-1079. 

I. fA. Akususxtl, “An outline of punched cards machines,” p. 1081-1120. 

L. I. GUTENMAKHER, “Electrical models of physical phenomena and their applications in 
technology and physics,” p. 1121-1146. In the literature list there are several references 
to American authors including two to BusH. 

L. Livsterntk, “Problems of computational mathematics,” p. 1147-1156. 
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L. fa. NeIsHuLER, “Tabulation of functions,” p. 1157-1176. There are here two references 
(p. 1157, 1175) to MTAC and to U.S. A. as “the country with the greatest development 
in the industry of calculating machines.” 

CuronicLE: M. L. Bykuovskxil, ‘The new differential analyzer of Bush,” p. 1177-1198, 
An illustrated description based on the long article of V. Bus & S. H. CaLpwe t, 
“A new type of differential analyzer,” Franklin Inst., Jn., v. 240, 1945, p. 255-326; 
see MTAC, v. 2, p. 89-91. 


R. C. A. 


NOTES 


68. Doctor ComriE’s ADDRESS.—We regret that we omitted to state 
in connection with the L. J. C. article, published v. 2, p. 149-159, which 
has been much in demand, that it was the address which he delivered 31 
October 1945 at the Conference on Advanced Computation Techniques 
(MTAC, v. 2, p. 65-68), as chairman of subcommittee Z of the Committee 
on Mathematical Tables and Other Aids to Computation. 


69. GIBBS’ PHENOMENON.—I feel that the Note on the sine integral in 
MTAC, v. 2, p. 195, will give the impression that any description of the 
Gibbs’ phenomenon in which the number 1.08949 (approx.) occurs is wrong, 
and that this number should be replaced by 1.17898 (approx.) as the result 
of a new and careful evaluation of K = (2/2)Siz. This is not so, as the 
appropriate number depends upon the way the phenomenon is described. 

The series F(t) represented by 3(4 — t) foro < t < x,and by — 3(4 +2) 
for —  <t <0 [not by }(¢ — =) as stated in the Note if we interpret ¢ 
algebraically] is an odd function with a discontinuity or jump of 7 from 
— 3x to $x at ¢ = 0. The Fourier series representing this function, i.e., 


> sin nt/n, 

n= 
exhibits the Gibbs phenomenon as an overshoot at each end of amount, 
say, 5, and, measured from the origin, the function jumps each way by an 
amount 32 + 6 which is given by 


Six 


f sin ¢ dt/t = $4 + siz, 
0 
where sir = — f sin ¢ dt/t ~ .28114, sothat 6 = siz. 


The phenomenon may be defined as in MTAC by the ratio K of 4x + 6 
to $m, ie., K = (2/x)[$e + si x] = 1 + (2/x) six ~ 1.17898. This ratio is 
also that of the jump including both the overshoots to the jump itself, i.e., 
a + 26 to x. 

We can, however, define the phenomenon by the ratio K’ of the jump 
+ either overshoot to the jump itself, i.e., * + 6 to 7, so that 


K’ = (1/r)(e + 6) = 1 + (1/x) si r ~ 1.08949. 


When the jump is not necessarily at the origin nor of amount 7 sym- 
metrically disposed about the ¢ axis, it is usual to describe the phenomenon 
as an overshoot at each end by an amount which is about 9% of the jump 
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itself,| a definition which follows naturally from the ratio K’ above. (See, 
for instance, T. v. KARMAN & M. A. Biot, Mathematical Methods in Engineer- 
ing, New York and London, 1940, p. 335). The ratio K is equivalent to 
saying that the sum of the two overshoots is about 18% of the jump. 

K. Knopp in his Theory and Application of Infinite Series (London, 1928) 
stated on p. 380 [p. 392 of the 1931 German ed.] that the first maximum has 
a value $7(1.08949) i.e. our K’, whereas his diagram on p. 379 clearly sug- 
gests that the value is really K. Obviously somewhere in the literature there 
has been a confusion between the definition of the overshoots as a per- 
centage of the jump, and Knopp has used the correct diagram with a wrong 
description. 

It is useful therefore to draw attention to the matter, but the Note in 
MTAC has obscured the issue by relating it to the accurate computation 
of the sine integral. Actually, to explain the difference between K and K’ 
there was no need to use the Harvard Automatic Sequence Controlled 
Calculator and prove K = 1.17897975. . . . If we eliminate si between K 
and K’, we have K’ = 3(K +1), so that we can obtain the interesting 
relation referred to in the footnote to the article in MTAC, without having 
to know the value of si z at all. 

G. MILLINGTON 


Marconi’s Wireless Telegraph Co., Ltd., 
Great Baddow, Chelmsford, Essex, England 


1 This result was stated by M. Bécher, Annals Math., s. 2, v. 7, 1906, p. 131.—Eprtor. 

Note sy R. C. A.: We are glad to have Mr. Millington’s communication which -will 
doubtless interest many readers. It seems desirable, however, to make clear that justification 
for his sweeping first paragraph (to which he returns in the last) is doubtful. I now refer to 
the following five places (which are the only ones) of my Note where the 1.08949 is involved: 
(i) Zygmund states that K = 1.089490, which is, of course, an error; (ii-iii) the two Sz4sz 
papers before 1944 (Amer. Math. Soc., Trams., v. 53, 1943, p. 440, and v. 54, 1943, p. 
497) where Zygmund’s incorrect value is copied. (iv) Harpy & RoGostNski state that 
k= fo" sintdt/t = 1.71 . . . (whenit should have been 1.85 . . .) which is really equiva- 
lent to Zygmund’s erroneous statement, as I pointed out. There remains, then, one and 
only one statement involving 1.08949, namely: (v) that of ZALCWAssER, who obtains the 
limit $4(1.089) which, as we remarked in the Note, isexactly the error of HARDY & ROGOSINSKI. 
Thus in every one of the “‘five places’, indubitable errors were listed. In not one of these 
errors enters the question of misinterpretation of ‘‘the way the phenomenon is described”’. 


70. THE GRAEFFE PRocEss.—In view of our previous articles on this 
subject by D. H. L. “The Graeffe process as applied to power series,” 
MTAC, v. 1, p. 377f, and by Mr. MiTcHELL, “‘The Graeffe process,” MTAC, 
v. 2, p. 57f, the following references to three articles published elsewhere may 
be noted: E. Bopewic, “‘On Graeffe’s method for solving algebraic equa- 
tions,” Quart. Appl. Math., v. 4, 1946, p. 177-190; Jose L. Massera, “El 
método de Graffe para resolver ecuaciones algebraicas,””’ Montevideo, Uni- 
versidad, Facultad de Ingenieria, Boletin, afio 10, v. 3, Dec. 1945, p. 1-20; 
R. San Juan, ‘“‘Complementos al método de Graffe para la resolucién de 
ecuaciones algébricas,”’ Revista Matem. Hispano-Americana, s. 3, v. 1, 1939, 
p. 1-14. 


71. WAS THERE AN ITALIAN REPRINT OF VEGA’S Thesaurus AFTER 1896?. 
—We may begin by quoting the following five authorities which state that 
there was: (a) In Jahrb. ai.d. Fortschritte d. Math., 1910, p. 1054, is the entry, 
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“G. Vega, Thesaurus, logarithmorum completus. Vollstandige Sammlung 
grésser logarithmisch-trigonometrischer Tafeln. Neudruck. Mailand. 684 
S. 4°.” (b) H. ANDovER, Nouvelles Tables Trigonométriques Fondamentales, 
Paris, 1911, p. vi; in listing the photozincographic reproductions of the 
Thesaurus, Istituto Geografico Militare of Florence, is the statement “un 
troisiéme tirage vient d’étre effectué (1910).”’ (c) Modern Instruments and 
Methods of Calculation, ed. by E. M. HorsspurGu, London, Bell, and Royal 
Soc. of Edinburgh, [1914], p. 50, 52, “reprinted Milan, 1909.’’ (d) F. J. 
DuartTE, Nouvelles Tables Logarithmiques . .., Paris, 1933, p. xxiv, ap- 
parently quotes (b) as his authority for a third reprint by the Istituto in 
1910. (e) FMR, Index, 1946, p. 440, lists, as in (b) and (d) a third Istituto 
reprint of 1910. A possible explanation of the 1909 date in (c) is that there 
was confusion with the sixteenth edition: of Cremona’s Italian translation 
in that year, of Vega’s Manuel logarithmique et trigonométrique (see Intern. 
Cat. Sci. Lit., v. 9A, p. 125. 

In spite of such an array of authorities I was puzzled that I could find 
(i) no mathematical bibliography except (a), and no Italian bibliography, 
which listed an Italian edition after 1896; (ii) not a single library which had 
a copy, and (iii) no review in any periodical. Hence I wrote to the Director 
of the Istituto Geografico Militare offering to purchase a copy of the 1909 
or 1910 edition. In his reply dated 2 Oct. 1946 occurs the following para- 
graph, in translation: ‘‘As a matter of fact besides not possessing any copy 
of the 1910 edition of Vega’s Thesaurus we cannot assure you that such a 
copy was an exact reprint of the 1896 edition, since the documents relative 
to it in the Library were destroyed by the Germans. From oral testimony 
gathered from clerks who were more or less directly concerned with the 
reprinting of the Thesaurus it would seem that the 1910 edition would have 
been in complete conformity with the 1896 edition.’’ As a result of this 
somewhat inconclusive paragraph, I published the rather vague statement 
about edition 4, MTAC, v. 2, p. 163. 

Shortly after reading this L. J. C. sent me a copy of a translation of a 
statement which he had received from the Istituto during 1922-24 (when 
he was on the staff of Swarthmore College), in reply to his varied queries 
including one about a 1909 or 1910 edition. The following sentences there 
occur: ‘“‘No other reprint of these tables has been issued by this Institute 
since the issue of the 1896 edition, and there is no record of any edition 
published in Milan in 1908 or subsequent to that date.” ‘“The zincographic 
plates prepared for the reprint are still preserved in the Institute.’’ L. J. C. 
also drew my attention to his published statement in the Observatory, v. 52, 
1929, p. 325 about these plates not being used for any edition subsequent 
to 1896. Hence I am now inclined to subscribe to the following statement 
by L. J. C. in a personal letter: “I will not believe in the 1909 or 1910 
edition until I have an affidavit from somebody who has seen one.” 

R. C. A. 


72. WHAT Is AN Error?—Those who delight in pointing out trivial end- 
figure errors may like to be reminded of the following words of wisdom from 
Glaisher’s pen (R. A. S., Mo. Not., v. 32, 1872, p. 261): ‘‘The increase of the 
last figure in tables, when the succeeding figures are greater than 500..., 
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seems to deserve more attention than it has received. Errata, such as some 
noticed in this communication, where the succeeding figures are 499. .., are 
by no means uncommon; and it appears that the discoverers of them imagine 
they are doing some service by noting them. Take, for example, one of the 
cases in this note: the figures starting from the tenth are 5 49998. ..; if we 
take 5 as the tenth figure, the error is 49998..., if 6, the error is 50002, 
differing by 00004. Now, as our table only professes to give 10 places cor- 
rectly (regard being paid to the magnitude of the figure in the eleventh place) 
a difference in the fifteenth place does not come in question at all: 5 and 6 
are both equally correct; they only differ by quantities, which throughout all 
the rest of the table we agree to neglect. It is a matter of regret that all such 
valueless refinements are not avoided by the author always explaining the 
exact convention on which the last figure is increased. A very convenient 
arrangement would be to understand that when x figures of a number were 
tabulated, the error was less than 6 in the next figure; or, if the calculator 
wished to be more accurate, 56 in the next two figures. To obtain a table of 
x figures, it is usual to calculate x + 1, or x + 2 figures, and the inconveni- 
ence of extending the calculation further in the particular case when the 
next figure is 5, or the next two 50, is, in many cases, excessive, and as the 
result is of no additional value when obtained, a figure “‘wrong’’ under these 
circumstances ought not to be styled an error. Probably a good many of 
Lefort’s errata are of this class. Babbage, in the introduction to his well- 
known table of seven-figure logarithms, states, that in ninety-three instances 
the next three figures in Vega were 500, and that in all these cases the log- 
arithms were carried to more than ten places to determine whether the 
figures were really 500... or 499. .., and decide whether the least figure was 
to be increased or not. This appears to me to have been quite needless. It 
sets up an unnecessary and artificial standard of accuracy for the numbers 
whose seventh, eighth, and ninth figures happen to be 4, 9, 9 or 5,0, 0. To 
the user of a table of seven-figure logarithms it is a matter of really no 
importance whether his error is 499 or 501; he is content to make an error 
of 5, and an additional error of + 0.01 is of no consequence.” 

With this I thoroughly agree—so much so that on more than one occa- 
sion I have written to our beloved editor saying ‘‘I have found . . . errors 
of less than one unit in . . . tables, but am not sending them to you, lest 
you should be tempted to publish them.” In the Introduction to the 
BAASMTC, Mathematical Tables, v. 6, Bessel Functions, Part I, which I 
edited, I wrote: ‘‘In general all values have been computed to two decimals 
more than are given in these tables; the error of any tabulated value should 
not exceed + 0.52 units of the last decimal.” Incidentally Professor H. H. 
AIKEN informs me that he has not found any error exceeding my assigned 
limits in the functions Jy and J;. It would be futile to go to the trouble of 
altering any end figures where the error lies between + 0.50 and + 0.52, 
since an interpolate is liable, in any case, to be in error by at least a 
whole unit. 

Errors, including those in the last decimal, often enable the table de- 
tective to ascertain how a table has been computed. One has only to instance 
the case of Buckingham, who refused to give any information about the 
compilation of his eight-place tables. But he was hopelessly betrayed by his 
errors (see MTAC, v. 1, 1943, p. 88f). The 2000 errors (exclusive of those 
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of a unit in the last decimal) in HayAsut, Siebenstellige . . . Tafeln (1926) 
give ample proof (a) that he did not check by differencing and (b) that he 
used a building-up process for intermediate values. Curiously enough, his 
end figures are fairly reliable. A run of errors in an early volume of Davis! 
showed (a) that he had made independent subtabulations in each interval, 
(b) that he had relied on repetition—the poorest possible check, and (c) that 
he had not checked by differencing. Other errors showed (d) that he had 
neglected second differences when interpolating 10-figure logarithms, and 
(e) that he had taken 10-figure logarithms of rounded-off quantities con- 
taining only five or six significant figures. But to his credit be it said that 
he was an apt “pupil” and can be trusted not to fall into any of these 
traps again! 

Mrs. GirFrorp’s end figures,” especially in the tangents, show the neglect 
of higher order differences; at one point there is a perfect wave in each 10” 
interval, with an amplitude of 3 units. The observation that her sines near 
90° were often in error by 99, 100 or 101 units led to a confession (in the 
true Sherlock Holmes style) that she “pre-fabricated” the first six decimals, 


and later added the seventh and eighth, with the not unnatural result that 


the sixth is often one out! 

DUuFFIELD’s claim to have computed his logarithms to 12 decimals, in- 
creasing the tenth when the last two were 50 or more, is immediately shown 
to be false by the fact that his end-figure errors are (with a few exceptions, 
which can be accounted for) the same as those of Vega!® 

The fact that BENson had copied from BRANDENBURG was revealed by 
his end-figure errors. He, too, was forced into a confession (MTAC, v. 1, 
1943, p. 9) that shows he had not been honest either in his compilation, or 
in his preface. Ives, who also wrote a deceitful preface, provided at least a 
part of the clue to his plagiarisms by his errors.‘ 

I have seen a 5-figure navigational table which contained just five per 
cent of end-figure errors, because it had been prepared from a six-figure 
table, but rounding off all 5’s in the same direction. 

| Ee ed 

EpiroriaL Notes: In MTAC, v. 1, p. 144 and 58, accuracy of half a unit in the last 
decimal place by PETERS, and PETERS & grein has been noted; and also on P. 145 accuracy 
less than .502 in the last decimal place. See further ‘ ‘Cayley and tabulation,” p. 98. In the 
quotation of a passage from Glaisher’s pen “‘Lefort’s errata” are those referred to in MTAC, 


v. 2, p. 164-165. The well-known 7D table of CHARLES BABBAGE is Table of the Logarithms 
of the Natural Numbers from 1 to 108 000, London 1827, and various later editions. 


1H. T. Davis, Tables of the Higher Mathematical Functions, v. 1. Bloomington, Ind., 
1933. 

2 E. Girrorp, see MTAC, v. 1, p. 11, 24f, 64f. 

%See MTAC, v. 2, p. 164 

4H. C. Ives, see MTAC, v. 1, p. Of. 
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21. PORTRAITS AND BIOGRAPHIES OF BRITISH MATHEMATICAL TABLE 
MAKERS.—Where may portraits be seen, or copies possibly be procured, of 
any of the following individuals: PETER BARLOW (1776-1862), HENRY 
BricGs (1561-1630), OLIVER ByRNE (publs. of 1838-77), ALLAN JOSEPH 
CHAMPNEYs CUNNINGHAM (1842-1928), James Dopson (d. 1757), RICHARD 
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FarLeyY (publs. of 1840-56), HERSCHELL E. Firtrpowsk1 (publs. of 1849-57), 
Emma GIFFORD (1861-1936), PETER Gray (1807?-1887), Epwarp LINDSAY 
IncE (1891-1941), HENRY SHERWIN (publs. of 1705-1741), RoBert SHor- 
TREDE (1800-1868), JoHN SPEIDELL (publ. of 1619), MicHAEL TAYLOR 
(1756-1789)? None of them are listed in Ball’s Collection of portraits (E. M. 
HorsBuRGH, Modern Instruments and Methods of Calculation, London, 1914). 
Karl Pearson tells us (Logarithmetica Britannica, part IX, 1924) that he 
vainly sought a portrait of Briggs. The Dict. Nat. Biog. contains sketches of 
Briggs and Gray, and there is biographical material about Cunningham in 
(a) Poggendorff, (b) London Math. Soc., Jn., v. 3, 1928 (by A. E. Western), 
and about Ince in (a) Who’s Who 1941, (b) London Mathem. Soc., Jn., v. 16, 
1941 (by E. T. Whittaker), (c) Edinb. Math. Soc., Proc., s. 2, v. 6, 1941 
(by A. W. Young), (d) Nature, v. 148, 1941 (by A. C. Aiken). Where may 
one find biographical data concerning the other persons listed? 
me Se 
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28. TABLES OF TAN~(m/n) (Q14, v. 1, p. 431; QR18, v. 1, p. 460; 20, 
v. 2, p. 62; 24, p. 147).—In MTAC, v. 2, p. 63, footnote 2, J. C. P. MILLER 
suggested the following problem: Determine all positive integers N which 
are fundamental in the sense that there is no relation of the form 


(1) tan-! N = }); tan n; 
(where the \; are integers and the n; are positive integers less than N). 

This problem can be completely soived using methods of elementary 
number theory (C. F. Gauss, Werke, v. 2, 1863, and 1876, p. 477 and 523, 
and the papers by C. ST@RMER and others quoted in MTAC, v. 2, p. 28). 
It can be shown that N is fundamental if and only if N* + 1 has a prime 
factor which is not a factor of a number n? + 1 with m < N. An effective 
construction for the relation of the form (1) in the case when WN is not 
fundamental can be given; a table giving all such relations with N satisfy- 
ing N? + 1 < 100 000 has been prepared. The positive integers which are 
fundamental form a minimal “‘integral’’ basis for the set {tan-!#} and 
therefore also for the set! {tan-!(m/n)}. The connection between this 
“integral’’ basis and the “rational” basis which was apparently known to 
Gauss? can be made clear by the following remark. Corresponding to a 
prime p = 2 or 4n + 1, Gauss has tan! (a/b) where a and 6 are determined 
uniquely by the conditions a 2 b > 0 and a? + b? = p while we have 
tan-! N where N is the least positive integer such that N? + 1 is a multiple 
of p. A table has been prepared, covering all such primes p < 500, expressing 
tan-! (a/b) in terms of tan! N and tan-! m with m < N. These results will 
be submitted to the London Math. So., Jn. 

The reduction of tan-' 1 40333 78718 was attempted in order to 
test the reduction algorithm. This is reducible since the prime factors of 
1 40333 78718? + 1 all occur as factors of m? + 1 with m < 1 40333 78718. 
Use of the factorisation given by Gauss (Werke, v. 2, p. 481) led to a contra- 
diction which was found to be due to an extra factor 13 in the decomposition 
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given by Gauss. The entry in his table corresponding to 1 40333 78718 
should read 
5-5-13-17-17-61-61-61-61-73-73-157-181 


and the required reduction is tan~ 1 40333 78718 = — tan 28 — 2 tan 27 
+ tan 19 — 4 tan“ 11 — tan“ 5 — 2 tan! 4 — tan“ 2 + 20 tan™ 1. 
Joun Topp 


Univ. of London King’s College 


1D. H. LexnMer, Duke Math. Jn., v. 4, 1938, p. 323-340. 
2C. F. Gauss, l.c., p. 523. See C. StgRMER, Archiv f. Math. og Naturv., v. 19, 1896, 
no. 3, p. 1-96, especially p. 77. 


CORRIGENDA 


A . 360, 1. — 4, for x = 1(2)6(1)10, read x = 1(1)6(2)10. 
A . 68, 275, for PEDERSON, read PEDERSEN. 


1,p 
2,p 
2 


. 2, p. 195, 1. 6, for $(t — x), read — $(t + x); p. 228, last line, delete ‘‘266, read 288,” and 
replace by ‘“‘Lewin, read Levin.” 
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CLASSIFICATION OF TABLES, AND SUBCOMMITTEES 


A. Arithmetical Tables. Mathematical Constants 
B. Powers 
C. Logarithms 
D. Circular Functions 
E. Hyperbolic and Exponential Functions 
Professor Davis, chairman, Professor ELDER, Professor KetcHuM, Doctor Lowan 
Theory of Numbers 
Professor LEHMER 
G. Higher Algebra 
Professor LEHMER 
H. Numerical Solution of Equations 
J. Summation of Series 





i Finite Differences. Interpolation 

K. Statistics 

Professor W1LKs, chairman, Professor COCHRAN, Professor EISENHART, Professor ‘ 
FELLER, Professor HOEL 


L. Higher Mathematical Functions 
M. Integrals 





N. Interest and Investment 
0. Actuarial Science 
Mister Eston, chairman, Mister Thompson, Mister WiLLIAMSON 


P. Engineering 





Q. Astronomy 
Doctor ECKERT, chairman, Doctor GOLDBERG, Miss KRAMPE 
R. Geodesy 
§. Physics, Geophysics, Crystallography 
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